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THE STELLAR TEMPERATURE SCALE’ 
A. PANNEKOEK 


ABSTRACT 


The methods of determining stellar temperatures from color and energy distribution, 
from ionization, and from the Boltzmann function are discussed and are shown to 
depend upon complicated conditions in the atmospheres of the stars which are not 
sufficiently known at the present time to permit the derivation of reliable temperatures. 
Direct methods of determining the effective temperatures of stars give 3300° for class 
cMo and 10,500° for class Ao. These values are regarded as well established, their un- 
certainty being estimated at no more than 5 per cent. It is suggested that from the 
yes. F degree of accordance between these temperatures and those obtained by other indirect 
ee methods information may be obtained concerning the conditions prevailing in the at- 

ee) mosphere of the stars. 


a 
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1. To find the temperature of a star, the most direct and the most 
promising method is to make use of the energy distribution in the ) 
spectrum. According to Planck’s law for black-body radiation, the 
logarithmic intensity difference for two light sources is a linear func- 
tion of 1/X, and the gradient is equal to 0.43 X 14,600/T, after a cor- 
rection factor (1 — e~*/*")-" has been applied, which is only impor- 
tant for high temperatures and long wave-lengths. Measurements of 
these gradients—relative to the sun, to a standard star, or to a light 
source of known temperature—have been made by Wilsing,? Rosen- 
berg, H. H. Plaskett,4 Sampson,’ Jensen,°® and by Greaves with his 

«Part of a symposium on astrophysics, presented at the Harvard Tercentenary 
Conference on Arts and Sciences, September 3, 1936. 

2 Publ. Ap. Obs. Potsdam, 56, 1909; 74, 1919. 

3 Nova Acta K. Leop. Akad., 101, No. 2, 1914. 

4 Pub. Dom. Ap. Obs. Victoria, 2, No. 12, 1923. 

5 M.N., 83, 174, 1923; 85, 212, 1925. © A.N., 248, 242, 1933- 
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fellow-observers at Greenwich.’ The first results showed serious dis- 
crepancies in the temperature scales. Careful comparisons and re- 
ductions of these data by A. Brill’ have shown that their chief cause 
was a change of the gradient with wave-length. That the logarith- 
mic intensity difference is not a linear function of 1/X is confirmed 
by the measures of Jensen. So it appears that the radiation of a star 
cannot be black-body radiation. 

Of course, there was no other reason to suppose that the stars ra- 
diate as black bodies than that in a first coarse approximation they 
seemed to do so. This hypothesis was introduced for the sound rea- 
son that it was the only means of deriving a temperature from the 
spectral intensity gradient. This temperature was called “color tem- 
perature’; it is of great practical importance because the most easily 
determinable data on color, such as color index and effective wave- 
length, depend on the same gradient and, therefore, are strongly 
correlated to this color temperature. 

The ionization theory afforded another method of finding stellar 
temperatures, by making use of the intensities of absorption lines. 
The state of ionization of an element depends on temperature and 
pressure. To eliminate the uncertainties of first appearance, of 
abundance and of absolute values, Fowler and Milne? introduced the 
method of maxima of absorption lines. Observations show the spec- 
tral class for which a line has its maximum intensity. By theory the 
temperature is computed for which the atoms producing this line 
have a maximum concentration in the state of ionization and excita- 
tion required to produce it. Three assumptions are involved in this 
deduction: the abundance of the element is constant over the spec- 
tral sequence; the absorption coefficient is constant for all spectral 
classes; and the pressure is the same in all these atmospheres. Be- 
cause this pressure is unknown, an absolute determination is not pos- 
sible. By the assumption that the Balmer series has its maximum at 
10,000° (Ao), the pressure was fixed; and it determined the tempera- 
ture of other maxima. The application of this method for lower tem- 
peratures is not easy, because there are few pronounced maxima. 
For high temperatures, however, just where the gradient ceases to 

7 Greenwich Obs., 1932; M.N., 94, 488, 1934. 

8 A.N., 218, 209, 1923; 219, 27, 356, 1923. °M.N., 83, 403, 1923; 84, 499, 1924. 
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give reliable results, it proved most useful. Especially through the 
work of Mrs. Cecilia H. Payne-Gaposchkin” the temperature scale 
of the B and O stars could be determined by the maxima of Sz, C, O, 
and N in different stages of ionization. 

These ionization temperatures may be considered as an extrapola- 
tion to high temperatures of a scale based on spectral energy distri- 
bution. Hence, the assumption of the stars radiating as black bodies 

which we know they do not—is underlying the entire temperature 
scale. 

Now it seems that such an assumption cannot be avoided. Even 
in the case of the sun we make use of it; the effective temperature is 
the temperature of a black body that emits the same quantity of 
energy. For this reason the color temperature is sometimes called 
another kind of effective temperature. There is, however, an impor- 
tant difference: the effective temperature measures the total stream 
of energy passing through a unit surface area of the star. This stream 
is determined by the interior conditions; it determines (for given 
conditions of composition of the atmosphere and of gravitation at 
the surface) the state of the stellar atmosphere, the temperature and 
the ionization in each layer, the spectral distribution of the emitted 
light, and the line intensities. The amount of this energy stream is 
the independent datum, and all the other quantities are functions of 
it. The effective temperature is a number indicating this energy 
stream. Therefore, the effective temperature is not an uncertain 
value founded on an incorrect assumption, but it is the exact expres- 
sion of a fundamental datum of the star. When we speak of the tem- 
perature of a star, this effective temperature is meant. Color tem- 
perature and ionization temperature have the meaning, that they 
may be considered as more or less correct substitutes for the effec- 
tive temperature. 

Table V in Brill’s discussion of the stellar temperatures in the 
Handbuch der Astrophysik™ shows the large discrepancies between 
the different determinations. The differences between the tempera- 
ture scales resulting from different comprehensive discussions, 
though smaller, still show their uncertainty, as may be seen from a 

t0 Harvard Circ., Nos. 252 and 256, 1924. 

1" Handbuch d. Astrophysik, 5, Part I, 150-151, 1932. 
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comparison of the scale derived by Mrs. Payne-Gaposchkin and 
quoted by Eddington,” the scale given by Russell, Dugan and Stew- 
art,’ and the scale given by Brill as the result of his discussions."4 
That these figures do not even mark the limits of uncertainty is 
shown in the result of the careful determination at Greenwich of the 


























Bo Ao Fo Go Ko Mo 
C. H. Payne-Gaposch- 
a ee eee 20,000 | 10,000 7500 5600 4000 3000 
Russell, Dugan, and 
Stewart... .........] 23,000 | 11,000 7400 5600 4200 (3100, M2) 
re | as eee 22,000 | 13,500 8550 5800 4370 3240 





color temperature by comparison with laboratory standards. Such 
absolute measures are difficult, and the elimination of the instrument 
introduced the uncertain element of atmospheric coloring. Still, the 
main result cannot be much in error, and it gave 18,000’ for the color 
temperature of Ao stars. 

2. The theoretical determination of the absorption coefficients of 
atoms under the conditions prevailing in stellar atmospheres, by the 
work of Kramers,'’ Sugiura,’ Gaunt,'? and Menzel and Pekeris,"™ 
provides a new foundation for the determination of stellar tempera- 
tures. It must be possible now, at least in the first approximation, to 
derive, first, the energy distribution in the stellar spectrum as a 
function of the temperature and, second, the variation of the line 
intensities with temperature for the real atmosphere, without the 
assumptions that were necessary ten years ago. Of course, theoreti- 
cal discussions on the structure of the atmospheres are needed, which 
can be given as yet only to a certain degree of approximation. For a 
wave-length in the continuous spectrum, the intensity is found by 
solving the equations for the integrated light: 

dJ dH 
dé es dé 

12 The Internal Constitution of the Stars, p. 141, 1926. "3 Astronomy, 2, 753, 1927. 

14 Handbuch d. Astrophysik, 5, Part I, 193, 1932. 

5S Phil. Mag., 46, 836, 1923. 

6 Scientific Papers, Inst. Phys. Chem. Res., Tokyo, 17, No. 339, 89, 1931. 

17 Phil. Trans., A, 229, 163, 1930. 1% M.N., 96, 77, 1935. 


= k(J — E); J = 2H for =o. 
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where £ = — f{pdh is the homogeneous depth, 47J = the total ra- 
diation, 47H = the net stream, 47E = the black-body radiation, 
k = the absorption coefficient, all holding for a special X. The so- 
lution for the energy stream emitted at the surface is: 


I " 
H = E()e“#V3 kV 3 dé . 
2+ Al | 


Denoting by & the mean absorption coefficient (Rosseland mean), 
we have 


T+ = Ti(t + $kE) (T, = the surface temperature), 
E= Cr (e@/AT —;)('= Cr (e@ ATo(t+3RE)— 2 __ 1)", 
ds \ ds : 


In my former computation, a linear function for E was assumed,” 
namely, the first two terms of a Taylor development of this expres- 
sion. Since the errors of this approximation may be rather large in 
some cases, we will use here the exact expression, reducing the inte- 
gral to the form used by Lindblad” and Milne." Putting 
bs V3 4 
kV 3dt = dt, T,  %? ni tol 


we have 


+ @ tf edt TC; ol 
4rH = 4 7 | reer eo Taam Aa 4 ( :) f(a, p) , 
2+V3N J, ealtts r 2+vV34 


sil ase~'dt 
f(a, p) = eT a 
For this function f(a, p) tables are given by Milne (using Lindblad’s 
data); since they, however, extend only from p = o to p = 2, they 
had to be extended to special large values of p needed here. The 
results for the net stream 47H are compared with the black-body 
radiation for the effective temperature 
1Cy 


rE(T,) = 5 


(ew—1), 


19 Tbid., 95, 520, 1935. 
20 Upsala U. Arrskrift, P. 33, 1920. 41 Phil. Trans., A, 223, 201, 1922. 
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4 
where a, = c,/AT; = aWV1/2. In the tabulated values the factor 
mc, has been omitted in 47H as well asin tE. The values for k/k 
were taken from the data of Amsterdam Publications No. 4 (for log 


g = 4.4). 

In Table I the results are given for effective temperatures T, = 
5040/0.2, etc., corresponding to T, = T,W1/2. They afford the 
differences in the gradient g = c,/T(1 — e~**")-" between the 
stellar and the black-body radiation for the intervals 6000-5000 A 
and 5000-4000 A. The differences from the former computation are 
not large, and the theoretical color temperatures show the same char- 
acter. The high color temperature of 18,000° found at Greenwich for 
the Ao stars, for which usually an effective temperature of about 


TABLE I 








| 
log p log 4H /c: log E/c: | Ag | 7) 














5040 
| | | | 
6000 | 5000 | 4000 6000 | 5000 | 4000 | 6000 | 5000 | 4000 | 5500 | 4500 | 5500 | 4500 
| | | | | | 

0.2.. —I orl—o 88|)—0.65| o 837] I.103} 1.416) o gro} 1.179| I 492|+0 02] 0.00] 0.90 | 0.79 
ee —0.76|/— .56)/— Fs 0.528, 0.773] 1.062] 0.611] 0.851| 1.118) — 03| — IO} 1.05 | 0.90 
0.4..../-0.37/— .16/+ .10] 0.296] 0.539] 0.838] 0.356] 0.564] 0.780/— .24/— .38] 1.06 | 0.86 
0.5....;—-0.14/-+ .07/+ .34] 0 ror! ©.342| 0.637] 0.122] 0.295| 0.458)— .47|/— .61| 1.07 | 0.87 
0.6.. —0.06/-+ .I7/+ .41] 9 900) ° 130| ° 389) 9.899] 0.037] 0.143}— .63}— .70] 1.16 | 1.05 
0.7....|/+0.08/-+ .20/+ 24| 9.733] 9.904] 0.019] 9.684] 9.783] 9.831)/— .50/— .31] 1.55 | 1.71 
0.8..../+0.14]/+ .08]— .09] 9.550] 9 606] 9.503] 9.472] 9.533] 9.520/+ .08/+ .41] 2.39 | 2.70 
0.9..../+0.08)— .05}— .31] 9.330] 9.290) 9.052] 9.262) 9.284) 9.211;/+ 43| + .76} 3.0% | 3.33 
16....|4+0 ——* 15|—0©.47| 9 103) 8 983) 8.635] 9.054] 9.035] 8 an tiie 79) +o 7 3.56 | 3.83 

















10,000 is assumed, is not only confirmed but is even surpassed by 
the theoretical values. Computation shows constancy and even tem- 
porary decrease of color temperature, in going from Ao to Bo, at 
variance with observation, which, after allowance has been made for 
space reddening, indicates a slow but steady increase. Hence it is 
not possible to derive effective temperatures from the observed gra- 
dients in the case of hot stars, so long as the qualitative discrepancy 
between theofy.and observation is not removed. 

In order to investigate it more closely, the energy distribution over 
the entire spectrum was derived for the case of an effective tem- 
perature of 10,080° (5040/T = 0.5). The data for the computation 
and the results are given in Table II. 

The variations in surface radiation depend on p = 0.87k/k; the 
mean absorption k determines the increase of temperature with 
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depth, and the special absorption k(A) determines from what depth 
we find radiation of this \ reaching the surface. For these high- 
temperature stars where the absorption coefficient is determined by 
hydrogen, it jumps at each band edge to a higher value on the blue 
side but then decreases continually with \% down to the next edge. 
The emitted radiation is strong on the longer wave-length side of 
each edge and falls there abruptly to a smaller value. (We disregard 
here the smoothing of the highest intensity peak at each edge in 
consequence of the line-crowding.) This small emission at the short 

















TABLE II 
| by as 
r p | a log 4H/cr | log E/cx 
BIO srt cn carne spe 0. 33 1.41 g.20-10 | 9.25—10 
gooo a alas ON 0.69 1.88 9.54—10 9.64—10 
& — 
8206 SnaWeenahat oe 2.06 els eS 
0.30 9.69—10 
7000 ©.49 2.42 9.91—10 9.95—I0 
4000 2.19 4.22 0.64 °.46 
3040 whl 4.63 oil Per Ore Ieee aT 
0.15 . —f2 
x i 0.27 a 
3000 0.27 5.63 pie 0.56 
2500 0.47 6.70 oe 0.54 
2000 °.gI 8.45 ‘ 39 0.41 
1500 sci Menten 2.17 11.27 45 ©.00 
: 0.62 > 9 
SOO a so Par Be Re eas 7.30 16.90 aa 8.82—10 
9.0 ) ” 
Ono tas sha Bee pt 18.54 1.09 8.42—10 
0.6 .—~ 10 
wave-length side of each edge is often spoken of as an “absorption 


band,”’ for which a correction should be applied to get the clear spec- 
trum. Figure 1 shows that this is not a correct statement of the 
phenomena. The whole spectrum is distorted; and corrections should 
also be applied, even larger ones, to the bright parts on the long 
wave-length side of the edges. 

The energy-curve of the spectrum emitted by the stellar atmo- 
spheres in Figure 1 shows, first, that in the Paschen region (between 
8206 and A 3646, where the third hydrogen-level determines the 
absorption) the gradient is stronger than in the black-body radia- 
tion; this is the phenomenon already treated, which appears in the 
color temperature. In the Balmer region (between A 3646 and X g11, 
where the second hydrogen-level dominates the absorption) the de- 
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viation is still more striking; here the emitted radiation increases in 
a steep ascent for the wave-lengths just above the Lyman edge, to 
dwindle to nothing below g11 A. These regions are, it is true, not 
observable; but this strong distortion of the spectrum must produce 
considerable effects in other phenomena. The radiations of wave- 
lengths 1215 A, 1024 A, 971 A (La, LB, Ly), producing the transi- 























4 1 1 4 4 4 1 | 1 1 
9/11, idea 
(1000 A); j 10 


FIG. 1 


tions from the lowest to the second, the third, and the fourth level, 
have here an intensity 22 times, 120 times, and 215 times larger than 
in the black-body radiation. Hence, in the upper layers of the at- 
mosphere the higher levels must be more strongly occupied than in 
the case of thermodynamic equilibrium, and the more so the higher 
the level considered. This must increase the general absorption; 
moreover, it offers an explanation of the remaining discrepancies in 
the color temperature of the A stars. 
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Since the transitions 1-2 and 1-3 chiefly dominate the occupa- 
tion of these higher levels, we may expect that the number of atoms 
in the third level is increased 5 to 6 times more than the number of 
atoms in the second level. Hence, the absorption due to ionizations 
from the third level is, in the same ratio, stronger than from the sec- 
ond level, relative to the former computations based on the Boltz- 
mann ratios. Since the mean absorption for this temperature” is 
chiefly a combination of second- and third-level absorption, this 
means that k/k and p must be smaller in the Paschen region and 
larger in the Balmer region. This will produce a smaller difference 
in intensity gradient between the stellar and the black-body spec- 
trum for the Paschen region than was found in our first computation 

a change in the right direction. It will not be easy, however, to 
determine the exact amount. What we computed here was the spec- 
tral-energy distribution in the radiation emitted at the surface of 
the star. For the deeper layers the composition of the radiation ap- 
proaches black-body radiation at an increasing rate; and all these 
deeper layers, with their different occupation of higher levels and 
their, consequently, different absorption coefficients, determine the 
real observed spectrum. Accordingly, in a second approximation a 
more refined treatment will be necessary in which the conditions in 
each atmospheric layer and its contributions to the total radiation 
are taken into account. 

There is still another point to be considered. We have assumed the 
surface temperature to be 7, = T, 1/2 (in this case 8470°), as 
follows from the Schwarzschild relation for black-body radiation. 
(The deviation from the linear gradient in the surface layers, giving 
the exact ratio 7,/T, = 1.23 as computed by Hopf, does not con- 
cern us here, because for us 7, means only a constant in the linear 
function we use at greater depth.) For the distorted energy-curve of 
the real stellar radiation, the same ratio will not hold. Now that it 
has been found in the first approximation, based on the surface tem- 
perature 7, = 8470°, we are able to integrate the radiation over the 
entire spectrum to find the total amount of emitted energy and to 


22 Cf. Table 8, Pub. Astro. Inst. Amsterdam, 4, 1935. 
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see to what effective temperature it corresponds. By reading from 
the curve and by numerical integration, we find 


Hd» 
og (4 = 0.290, 
Cr 
E wt (T,\4 
log f! dX = log - ( = 0.20 


whereas 


| 
No 


un 


when computed with 7, = 10,080°. The total amount of energy 
comes out somewhat larger for the stellar atmosphere and corre- 
sponds to an effective temperature somewhat higher, 10,400’. This 
means that for a star with effective temperature 10,080° the surface 
temperature is somewhat lower than that computed by the Schwarz- 
schild relation. The difference is not large enough to upset the fore- 
going results; it has to be taken into account in a second approxi- 
mation. 

3. For the stars of low temperature there are difficulties of another 
kind. The absorption coefficient is due to the ionization of metal 
atoms from different levels. Because we do not observe any band 
edge corresponding to a metal-atom level, we assume that there is 
some smoothing process at work; and for the absorption coefficient 
we take the value computed by a smoothed formula: 


k = Const. T~?y—3(e#/*T — 1). 


The computations were made for the temperatures 5040 and 3150 
(5040/T = 1.0 and 1.6). The result is contained in the values of 
Table III. The absorption coefficient is a minimum for a large in- 
frared wave-length (11,000 A for 5040, 17,000 A for 3150°), and 
increases strongly with decreasing wave-length in the visual part of 
the spectrum (because at low temperatures the lowest levels are 
strongly occupied relative to the higher levels). The consequence is 
that the short wave-length part of the spectrum, including the visible 
region, is much weakened (as is shown in Figure 2 for 7, = 3150), 
whereas the maximum in the near infrared is increased and dis- 
placed toward the longer wave-length side and the far infrared again 
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is weakened. The slope in the visible part of the spectrum deviates 
strongly from the corresponding black-body radiation and indicates 


TABLE III 




















Tr = 5040 T, = 3150 
| | | | 

a | p | log 4H /cy | log E/cy a p | log 4H /c; | log E/c, 

50,000 0.68 0.30 | Le a | oe Pe 
25,000 1.35 | 0.85 | 7.69 7.68 2.16 | t.45. |. 746 7.30 
20,000 1.69 | 1.12] 8.05 | 8.00 2.70 | 1.66 7.04 a ee 
15,000 2.25] 1.48] 8.43 | 8.37 3.61 | 1.74 | 7.94 7.82 
12,000 2.62 | 2.70 S.7% | 8.62 4.§8 1 1.85 8.08 7.97 
gooo 3.76 | 1.70 8.99 | 8.88 6.01 | 1.02 8.08 8.04 
7000 4.83 | 1.45 Q.12 Q.02 y ay ee) eo 7.80 7.95 
6000 5.63 | 53 Q.II 9.05 g.O1 | 0.28 7.48 7.81 
5000 6.76 | 0.76] 9.00 | 9.04 | 10.82 | 0.10 | 6.95 7.55 
4000 8.45 | 0.36 8.64 | 8.90 13.52 | 0.02 | 6.18 7.05 
3000 11.27 | 0.08 7.84 | 8.49 18.03 | 0.001 | 4.85 6.03 

| | 
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far too low a temperature. All these deviations are strongest for the 
lowest temperatures. Integrating the energy values of Table III 
over the whole spectrum, we find here again a total stream of energy 
somewhat larger than belongs to the corresponding 7, (log = 9.032 
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instead of 9.001 for T, = 4230, 8.221 instead of 8.185 for T, = 
2650) and corresponding to a black-body radiation of 5130° and 
322c°, only slightly larger than the values of 7, V2. 

We can make use of the measurements of Jensen to derive tem- 
peratures from the spectral gradients of low-temperature stars. Jen- 
sen measured for a range of wave-lengths the relative intensities of 
a number of stars relative to a Ursae Minoris. By means of seven 
stars, which he has in common with the Greenwich observers 
(n UMa, a Leo, a Lyr, a CrB, a Agl, a UMi, a Aur), we can use the 
absolute Greenwich gradients to reduce his results for a Orionis to 
absolute values (i.e., intensities relative to a black body of T = ©) 
and compare them with H and E multiplied by A‘. The result of 
the comparison is shown in Figure 3. (The shorter wave-lengths, from 
4000 A downward, are not reliable because here in the class A com- 
parison stars the hydrogen depression makes itself felt.) Compared 
with the black-body-curve, the slope of the observed intensity-curve 
indicates a temperature below 3150°. On the contrary, the theoreti- 
cal curve computed for the stellar atmosphere shows a steeper slope 
than the star. Hence, 7, for the star must be above 3150°—or, more 
exactly, above 3220 —and the theoretical curve belonging to T, = 
3600° seems to fit better. An exact result is not possible because 
Jensen’s measurements show a variation of the gradient with wave- 
length which is not present in the theoretical curves. 

There are other sources of uncertainty here. It is doubtful whether 
for these low temperatures we can rely upon the smoothed formula 
for the absorption coefficient where it is due only to the ionization 
of a small amount of sodium and potassium atoms. Even if we were 
able to observe and measure the continuous background of the spec- 
trum, its gradient would give information about the variation of the 
absorption coefficient rather than about the temperature. But it is 
well known that, owing to the crowding of absorption lines, chiefly 
those constituting the molecular absorption bands, the sifting-out of 
the continuous spectrum is well nigh out of the question. Thus, for 
such low-temperature stars the gradient of the visible spectrum, ob- 
served with present means, cannot give very reliable results for the 
temperature. 
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4. The theoretical knowledge of the coefficient of continuous ab- 
sorption must provide the foundation for the use of line intensities 
for the derivation of stellar temperature. The intensity of an absorp- 
tion line depends on the ratio s/k of two quantities: s, the diffusion 
coefficient in the line (depending on abundance, state of ionization, 
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level, and transition probability of the atoms involved); and k, the 
continuous absorption coefficient. Because k is variable with tem- 
perature in a known way, the maximum intensity of an absorption 
line is not simply identical with the maximum concentration of the 
atoms producing it, but one can be derived from the other. More- 
over, it is not necessary to assume a pressure; all the layers in an at- 
mosphere with different pressures take part in the formation of the 
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line, and the theoretical treatment of the atmosphere gives the line 
intensity, where the intensity of the total stream of energy (i.e., the 
effective temperature 7) and the surface gravity g appear as the 
determining parameters. The only assumption we cannot avoid mak- 
ing use of is that of constancy of abundance over the range of spectral 
classes. 

The most conspicuous of the line maxima used in determining 
stellar temperatures is the maximum of the Balmer series of hydro- 
gen in class Ao. Fowler and Milne were able to explain it easily be- 
cause the hydrogen atoms in the second level of energy have a clearly 
depicted maximum of concentration, for which 10,000° was assumed. 
If, however, the variation of k is taken into account, we get an ex- 
tremely flat and ill-defined maximum at a lower temperature (7000° 
to 8000” for log g = 5 to 2; cf. Figure 6 in Amsterdam Publications 
No. 4). Unsold has shown? that in a pure hydrogen atmosphere the 
intensity of the Balmer lines must steadily increase for decreasing 
temperature. The reason is simply that the diffusion coefficient s for 
these lines depends on the concentration of H in the second level, 
whereas they are situated in the Paschen region, where the absorp- 
tion coefficient k depends on the concentration of H in the third level, 
and that the relative number of the atoms in the second state com- 
pared to those in the third state increases with decreasing tempera- 
ture. Going down below 10,000° the intensity of the Balmer lines 
must go on increasing because, though JN, and s are strongly de- 
creasing, NV, and & are decreasing still more strongly. As soon, how- 
ever, as k for hydrogen falls below the absorption coefficient due to 
the metal atoms, the & of the stellar atmosphere begins to increase; 
and now the intensity of the Balmer lines decreases rapidly. Their 
maximum, therefore, is not due to themselves but to the transition 
of the continuous hydrogen to the metal absorption. If the relative 
abundance of the hydrogen to the metal atoms is assumed to be 50 
instead of 1000, the maximum is raised 1000°. It is not likely that in 
this way we can escape the difficulties. Moreover, the intensity of 
the Balmer lines is an expression chiefly of their width, which is due 
to the Stark effect and depends, in the first place, on the surface 
gravity. Hence, in comparisons with observational data a careful 


S25. f ADp.,'8, 32, 1034. 
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arrangement according to the gravity g is first necessary before the 
variation with 7, can come out in the right way. We may expect 
that further approximations in the theoretical treatment of the A 
stars may clear up the discrepancies; but, for the moment, the maxi- 
mum of the Balmer lines is not a good indicator of temperature. 

This hydrogen case is an example of a more general phenomenon. 
Between 10,000° and 6000° the continuous absorption coefficient 
varies irregularly with temperature, because here is the transition 
from hydrogen to metal absorption. The variations of line intensi- 
ties in this realm are determined more by these irregularities than 
by the variations in the element producing the line. The exact 
course of these irregularities depends on the assumed composition 
of the atmosphere, especially on the relative abundance of hydrogen 
and metal atoms. In the maxima of absorption lines within these 
temperatures, composition and temperature are interrelated. 

Below 6000° and above 10,000° the absorption coefficient runs 
smoothly and slowly with temperature, so that the maxima are re- 
lated in a simple way to temperature. Lines absorbed by the lowest 
level of ionized atoms (as H and K of Ca*) present very flat and un- 
certain maxima; hence they cannot give exact results for the tem- 
perature. The matter is different for lines absorbed by higher levels, 
where the temperature of the maximum depends chiefly on the ex- 
citation potential. Here, however, the crowding and blending of lines 
make determinations difficult, and exact photometric measures of 
line intensities are lacking on the whole. Moreover, the variation 
of k with temperature depends here on the assumed composition of 
the metal atom mixture. 

The case is more favorable for high-temperature stars above 
10,000°, where the absorption coefficient depends on hydrogen only 
and the absorption lines can be separated easily. Thus for the B 
and O stars the former results of Fowler and Milne and of Mrs. 
Payne-Gaposchkin from metalloid lines can be put upon better 
foundations. After the method indicated and used in Amsterdam 
Publications No. 4, the line intensities have been computed for a 
range of temperature, then plotted, and the temperature of the maxi- 
mum read. The pressure and the atmospheric parameters s, o, 0’ 
were taken for two cases of surface gravity, log g = 4.4 (dwarfs) 
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and log g = 2.4 (giants). In Table IV the results for the tempera- 
tures of the maxima are given, and in the last columns the spectral 
classes of these maxima are indicated, first according to old results 
of H. H. Plaskett and Mrs. Payne-Gaposchkin, and then according 
to the latest results of E. G. Williams.** The resulting scale of tem- 
peratures is not very different from the former investigations; for 
Bo the corresponding temperature is nearly 25,000”. 





TABLE IV 


: H. H. Pi. on ; 

Ion. P. Exc. P. log g = 4.4] logg = 2.4 pe ee E. G. W. 
Si 16.3 9.8 12,900 II,200° Ao 
( 24.3 18.0 21,000 16,800 Sea ohat os B3 
Nt. 29.5 18.4 21,000 18,700 B2 31-2 
Set 22 2 18.9 22,000 19,400 Bi-2 Bi-2 
He. 24.5 20.9 18,000 14,400 B2-3 
O* 34.8 23.3 24,000 21,000 Bi Bi 
S72 44.9 22.6 31,500 27,000 Oo Oo 
\ 47.4 27.3 35,000 31,500 O7-8 Og 
( 47.7 290.4 30,000 28,000 Og? Og-Bo 
O 54-9 33-7 37,400 33,000 Og-Bo O8 
N® Ye 30.7 39 ,000 33,6000 ae tane 
Het. 54.2 50.8 44,000 39,000: }........ 
Gere 64.2 55-5 50,000: | 44,000 























5. Now that the intensity of a line expressed, for example, as the 
half-width for 50 percent residualintensity, or as the equivalent width, 
can be measured as well as computed, it is not necessary to make use 
of maxima. The intensities themselves may be used, and so a larger 
quantity of data may be utilized. The number of good intensity 
measures, it is true, is only small as yet. We can, therefore, make 
use of a few cases only as instances of the method. 

For the line K of Ca* the half-width in the solar spectrum is 
6.3 A; this corresponds (according to formulas 28 and 46 in Amster- 
dam Publications No. 4) to log abundance = —4.71, that is, the 
number of Ca atoms is 20, if all the metals together are taken as 
toco. By means of this abundance we compute log AX’ and the 
half-width Ad itself (in the case of disappearing atoms): 


5040/T = 0.4 0.5 0.55 0.6 0.65 
log AA? = —2.58 —1.04 —0.13 +0.63 +1.28 
Ar = 0.051 0.30 0.86 2.1 4.4A 


24 Ap. J., 83, 319, 1930. 
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From some slit-spectra of a Lyrae, taken in August 1935, with a 
two-prism spectrograph attached to the Wyeth reflector at the Har- 
vard College Observatory, I find an equivalent width of 0.97 A. 
From theoretical profiles, we found for the relation between half- 
width and equivalent width: 


log E.W........ 0.40 0.20 0.00 9.80 9.60 9.40 9.20 9.00 8.80 
loe EL.W........ 9.89 9.69 9.48 9.28 9.06 8.95 8.78 8.68 8.54 


Thus, an equivalent width of 0.97 A means a half-width of 0.29 A, 
corresponding to 5040/7 = 0.50. Hence, assuming for a Lyrae the 
same abundance of Ca atoms as for the sun, we find a temperature 
of 10,100”. 

In the same way we can make use of the results of Mrs. E. T. R. 
Williams-Vyssotsky”> for the equivalent width of the line K in a 
number of class A stars. We find (taking log g = 4.4): 


| | 
| | 
| Sirius Mean Ao A2 A3 | As 
E. W. ..| 0.62 | 1.08 2.94 | 4.03 Aue 
i err. ma 0.33 0.91 | 1.25 1.39 
5040/T | 0.47 0.50 0.55 0.57 0.58 
i. | 10,700 | 10,100 200 8800 8700 


It must be remarked that the assumed relative abundance of 
hydrogen and metal atoms goes into this result in its full amount, be- 
cause it determines the relative values of k& for A stars and the sun. 
If hydrogen atoms are 100 times (instead of tooo times) more abun- 
dant than metal atoms, the temperatures deduced are raised 1200”. 
Moreover, the assumption of equal abundance of calcium in the sun 
and the A stars is involved, whereas after the results of W. W.Mor- 
gan” the abundance of different metals must be different in different 
A stars. The deviations, however, from constant composition do not 
seem to amount to large factors; since the variation of the concentra- 
tion of active atoms with temperature is very rapid, the uncertain- 
ties of temperature resulting from these differences will not be large. 

There may be a way to escape these uncertainties and to eliminate 
the differences of abundance and of absorption coefficient at the 
same time, namely, by comparing the line intensities for ionized and 


25 Harvard Circ., No. 348, 1930. 7° Pub. Yerkes Obs., 7, No. 3, 1935. 
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neutral atoms of the same element. The line \ 4227 of Ca and the 
K line \ 3933 of Cat both belong to the lowest level; therefore, no 
other factors are involved than the variations of k with wave-length, 
which are small for this small distance in \ and are neglected in the 
following computation. For a Lyrae the equivalent width of Ca 
4227 was found from the Victoria spectra to be 0.044 A. Accord- 
ing to the curve of growth derived by Struve and Elvey for this 
star,?’ this value is already situated on the wing part of the curve, 
where the equivalent width increases as the square root of the con- 
centration. Hence, the ratio of the concentrations of Ca and Cat 
atoms is the square of the ratio of the equivalent widths of \ 422 

and A 3933, 1.e. (0.044:0.97)?, of which the logarithm is — 2.70. 
For the sun we have the half-width of \ 4227, taken from the pro- 
file derived by Redman,” 0.57 A; compared with the half-width 
6.3 A of the K line, we find the logarithm of the ratio of the con- 


centrations log (0.57:6.3) = —2.08. Hence, we find (x = rate of 
ionization): 
c x 
log a Lyr : ©] =0.62. 
1—x/) ~*~ I-2 
Computing this ratio after the ionization formula «/(1 — x) = 


K/P (taking log g = 4.4), we find: 
T = 10,080° 8400° 7200° 6300° 5800° (©) 
log K/P = 4.61 3.89 *, 38 2.83 2.46 


The observational result would give a temperature for a Lyrae 
nearly 1000° above the solar temperature. In this Ao star both 
\ 4227 and K have strongly decreased relative to the sun, but nearly 
in the same ratio, whereas theory demands that \ 4227 should have 
decreased much more rapidly. Indeed, for a temperature of 10,000° 
we should expect \ 4227 to be nearly 7 times fainter, so that it would 
have entirely disappeared in the spectra of Ao stars. Of course, we 
could assume that the line seen at that wave-length is not the Ca 
line 4 4227 but has some other origin; indeed, a faint 77+ line 
corresponds in wave-length. We shall see, however, that the same 

77 Ap. J., 79, 409, 1934. 

8 M.N., 95, 751, 1935. 
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phenomenon presents itself for other elements; hence, it is not 1ikely 
that we can so easily explain away the discrepancy. 

6. Adams and Russell derived”? in 1928 the scale of stellar tem- 
peratures by comparing the number of atoms in different levels with 
the Boltzmann function. The factors of abundance and transition 
probability for each line were eliminated by using for each star rela- 
tive values compared with the sun. Then the logarithm of the con- 
centration ratio is a linear function of the excitation potential E: 

log Y = Const. — B( 54° — ©) 

1 1 

The relative concentrations (or numbers) of active atoms producing 
the lines were found by comparing the stellar with the solar-line in- 
tensities and calibrating the Rowland scale of solar intensities by 
means of multiplets. The Fe spectrum bears nearly go per cent of 
the weight, because here the levels extend continuously from E = o 
to 4.5 volts. It appeared that for a Orionis and a Scorpii, log Y was 
not a linear function of EZ; for large E the variation in log Y slowed 
down. So it was necessary to replace £ in the formula by a quad- 
ratic expression E — 0.083’. This means that the limiting slope 
of the curve for E = o was used to derive the difference in 5040/7. 
The same quadratic expression was then used for the other stars. 
In this way temperatures of gg0o° for Sirius and 2900” and 3000° 
for a Ori and a Sco were found. 

Here the line intensity, after calibrating, was directly connected 
with the number or the concentration of atoms. We now know that 
line intensity is directly connected with s/k, so that numbers ex- 
pressing line intensities must be multiplied by factors k to find num- 
bers expressing the concentration of atoms. Hence, corrections = 
log (k../kR©) must be applied to log Y. The values of log k were 
taken from Amsterdam Publications No. 4, using 5040/T = 1.6, 
log g = 1.4, for a Ori and a Sco; 5040/T = 0.5, log g = 4.4, for 
Sirius; and 5040/T = 0.88, log g = 4.4, for the sun. The differences 
Star minus Sun vary with wave-length; Dr. Russell kindly provided 
me with the separate wave-lengths of the lines used, so that for each 
group the reduction could be computed. In Table V the & correct- 


29 Ap. J., 68, 9, 1928. 
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ions and the uncorrected and the corrected values of log Y are 
given. It appears that by thus taking account of the variations of 
the absorption coefficient with temperature and wave-length, the 
constant part of log Y is changed considerably, whereas the slope is 





TABLE V 
| a ORIONIS, a SCORPII SIRIUS 
|—_—— — —— 
E 
, log Y : log Y 
log } log Ak | : log } log Ak 

| corr. | corr 

ee eee + een, = | 
0.06 | 2.05|\+2.30/—1.30/—1.40] 0.75 | 0.90 |—2.66/+2.13] 9.47 
0.94 +1.28/-+1.23) 1.54) 1.571 9.74 | 9.66 2.13) 2.20) 0.16 
1.54 I+-0.54/-+0.97| 1.33) 1.42] 9.21 | 9.55 2.30} 2.11] 9.81 
2.21 .|+0.30/+0.47] 1.44| 1.35) 8.86 | 9.12 2.00] 2.25] 0.25 
2.56 =0.00/—0.05| 1.67] 1.70) 6.27 | 3:25 2:00| 2:43) O:43 
2.89 —0.19]/—0.03] 1.40) 1.41] 8.41 | 8.56 E2Os| 2.22) -0.27 
2.2 |-° 23|=0.10) 1.41 1.43} 8.26 | 8.47 EuQ9| .2.231°0:.20 
ue Gy -|==O.54)--0.27] F252] 1-50) 7:04 | 3.23 E791 “289)- 0/85 
4.03 —0.76)/—0.40) 1.57) 1.57) 7.67 | 8.03 1.77} 2.36] 0.59 
4.50 =6,57) 0277! — 1. 50|-—1 60] 7.55 | 7-63 |—1.60)/+2.38] 0.78 

| | | 

| | 








changed by only a small amount. Figure 4, where the corrected 
log Y are represented, shows that the curvature for a Ori and 
a Sco still exists. Retaining Russell’s formula, and hence assuming 
that in all the stars the occupation of the higher levels follows 
E —o.083 E’, we find the expressions 

log Y = 0.88 —o.82 (E — 0.083?) for a Ori and a Sco, 

log Y = —o.50 +0.41 (E — 0.0837?) for Sirius . 
Assuming that the temperature determining the observed distribu- 
tion over the levels coincides with the effective temperature, we 
find the values of 5040/T = 1.70 and 0.47; hence T = 2960° for 
a Ori and a Sco and T = 10,700° for Sirius. We have to consider 
here that the quadratic expression is only an interpolation formula 
which certainly does not hold for higher values of F, so that the 
limiting slope for E = o, on which depends the temperature found, 
has the uncertainty of an extrapolated value. 

The quadratic deviation from the Boltzmann distribution, the 
so-called ““Adams-Russell phenomenon,” has not found an adequate 
explanation as yet. Dr. Russell, in a letter, suggested the possibility 
that for these two Mo stars the value of log Y for the lowest E may 
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be systematically wrong, because the lines belonging to the lowest 
levels are extremely strong in these stars and so cannot be exactly 
compared with the solar lines. If the result for E = 0.06 is excluded, 
the remaining values, as Figure 4 shows, may be represented by a 
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straight line, which, of course, has a smaller slope. For Sirius, log Y 
as a function of E is as well represented by a straight line as by the 
curve of Adams and Russell. The resulting linear expressions are 


a OrI, a SCO — SUN SIRIUS — SUN 
log Y uncorr. = +1.60 —0.57E log Y uncorr. = —2.55 +o.21E 
log k/kO = —1.40 —0.03E log k/kO = +2.15 +0.05F 
log Y corr. = +0.20 —o0.60# log Y corr. = —o0.40 +0.26E 
Then for 5040/T we have 1.48 and 0.62; hence T = 34c0° for a Ori 


and a Sco, and T = 8130° for Sirius. Because the mean slope used 
here is smaller than the limiting slope in the quadratic formula, the 
scale of temperatures is narrower. 

The data of Adams and Russell used here can give still another 
result. The constant in the foregoing formulas represents the rela- 
tive number of atoms for the star and the sun for E = 0, 1.e., for 
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the lowest level; it depends only on abundance and state of ioniza- 
tion. Comparing these values for Fe and Fet, the abundance of the 
element is eliminated in their ratio; and the same holds for the values 
of k, inso far as we may assume that the average wave-length for the 
Fe and the Fet lines was the same. Hence, in the same way as in 
the case of Ca and Cat‘, this ratio is equal to the ratio of x/(1 — x) 
for the star and sun. The values of log Y in the case of Sirius given 
by Adams and Russell for Fe*+, 77, and Ti+ were first corrected for 
k and then reduced to E = o by means of the slope derived from Fe. 


TABLE VI 














| | | 
g , > og } 
I | log } log } —o.26E) log Vo | E log } , 6E| log Yo 
| corr. corr. 
Fe" | Ti 
a -. « | 
2.5... —0o.18|/+2.02}—0.73/-+1 29|| 0 03 —2 00/-+0 20/—0.01|/+0c.19 
3.56 '—9.54/-+1.66]/—0.95|+0.71]|] 0.83 2.11] 0.09 22|— .13 
| | | 1.02 1.07| 1.13 27\+ .86 
a es oy” verre 1.13] 1 07 | 371+ .80 
Tit 1.82 1.82] 0.38} -47|— .09 
ROD. csv 1.00] 1.20] s4i+ .66 
| 6 ; . 2.26 —0.90/+1.30}—0.59/+0.71 
1.16.......{—1.05/+1.15]—0.30|+0 85| _ | | 
bi87 |—0.78] 1.42] 41| 1.01] 
2.00 ...,-0.08] 2.28] 52| 1.76) 
2.93 =o 75|+1.45|—0.76|+0 69} 
| | | 7 
Taking averages for log Y., the ratio log Yo — log Y? = log 


[x/(1 — x) St. : x/(1 — x) O]. We have 
log Y, for Fe —o.40, for Fe+ +1.00 log Yo /Y? =+1.40; 
log Y, for Ti +0.43, for Ti+ +1.08 log Vi/Y"?= +0.65. 


Computing log x/(1 — x) = log K/P for different temperatures 


(log g = 4.4), we find: 


10,080 | 8400° | 7200 | 6300° | s5800°(©) 
Fe 3.68 | 2.78 3 | F327 | 0.87 
: 2.10 | 1.04 
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The corresponding temperature for Fe is 7500°; for Ti, 6600”. 
Hence we meet with the same phenomenon as in the case of Ca: 
the decrease of the lines of the neutral atoms relative to the en- 
hanced lines from the sun to Sirius is far too small for their large 
difference in temperature. Assuming a temperature of 10,000° for 
Sirius, all those arc lines which we see in the Sirius spectrum should 
be invisible after the ionization formula, provided the right curve 
of growth has been used. But these faint lines probably are situated 
on that part of the curve of growth where a large variation of atom 
concentration corresponds to a small change in line intensity. So it 
is probable, here as well as in the case of Ca, that pecularities of the 
curve of growth, depending upon turbulence and collisions in the 
star and upon the Weisskopf-Wigner corrections of the atoms, play 
a rdle and make the exact determination of temperatures by this 
method difficult. These computations of Y, and of their ratio have 
also been made with the quadratic formula; they give so nearly the 
same results that it is unnecessary to give them in detail. For the 
low temperature Mo stars the ratio log Y{/Y? = —o.55 for Fe 
and —1.32 for 77. It is not well possible to use these results for the 
derivation of temperatures because, owing to the chromospheric 
character of their atmospheres, we do not know what pressure to 
use in the computation of the ionization. On the other hand, in so 
far as these relations are trustworthy, they may better be used to 
derive some information about the pressure in the relevant atmos- 
pheric layers, once their temperature is known. 

7. The effective temperature itself is derived for the sun by meas- 
uring the total radiation, the solar constant, and reducing it by 
means of the angular semi-diameter of the sun to radiation per unit 
surface. For the stars this method can be applied in a direct way 
only, where the angular diameter can be determined by interferom- 
eter measures, i.e., for red giants and supergiants. The total radia- 
tion, the bolometric magnitude, is then derived from the visual mag- 
nitude by means of the ‘“‘reduction to bolometric magnitude,” such 
as is given by Eddington:” for 3600°, —o.95 mag.; for 3000°, —1.71 
mag.; for 2540° —2.59 mag. Then for a Sco we find from the angu- 
lar semidiameter of 0%021 and from the visual magnitude 1.2 the 


3° Op. cit., p. 138, 1926. 
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temperature 298c0° (red. to bol. mag. = —1.74; bol. mag.= —o.54). 
Eddington’s table of reduction from visual to bolometric magnitude 
has, however, been computed for black-body radiation. Hence, this 
derivation does not materially differ from the computation made in 
the textbook of Russell, Dugan, and Stewart, where the black-body 
radiation for \ 5290 is compared with the visual magnitude of the 
star. A better approximation will be reached if, for the reduction 
from visual to bolometric magnitude, we make use of the spectral 
intensity-curve computed here for the stellar radiation. The curve 
for 3150 in Figure 3 gives a reduction of — 2.70 mag., i.e., 1.2 mag. 
larger than the black-body-curve gives for the same temperature. 
Making use of the interpolated values: for 3360°, — 2.02 mag.; for 
3600°, —1.46 mag., we find now from the same data on a Sco the 
temperature 3310° (red. to bol. mag. = —2.18). 

Though in this way the systematic error of former derivations is 
avoided, the method remains somewhat unsatisfactory because the 
total radiation is derived from the measurement of a minute part of 
it at the extreme border of the spectrum. To find the effective tem- 
perature directly, bolometric or radiometric measurements of the 
total radiation are needed. Pettit and Nicholson have measured 
stellar radiations with vacuum thermocouples,** and by means of 
absolute standardization of their instrument in various ways they 
could derive stellar temperatures from the interferometer diameters. 
They found 3270° for a Ori and a Sco.? The measures had first to 
be reduced to “‘no atmosphere,” i.e., corrected for the atmospheric 
absorption bands in the red and infrared; this correction was derived 
from measures on solar energy-curves, of course, by supposing the 
stellar radiation to be black-body radiation of a certain low tempera- 
ture. This involves the possibility of a certain error; the radiation of 
the star has a different energy distribution from a black body, and for 
its derivation even the knowledge of k may not be sufficient because 
of the strong molecular absorption bands in the stellar spectrum it- 
self. Still, it is not likely that appreciable errors will arise, because 
the atmospheric bands are well distributed over the most contribut- 
ing parts of the spectrum. Hence this determination deserves a large 
weight. We might suppose that a certain limb-darkening of these 
stars, now theoretically determinable from the knowledge of k, could 


3 Ap. J., 68, 279, 1928. 32 [bid., p. 301. 
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produce an important correction to the angular diameters derived 
from interferometer measures. The figures of Table III, however, 
show that p and k/k for visual rays is very small; hence there is no 
appreciable limb-darkening here. 

The accordance of the two results based upon interferometer 
measures is very satisfactory now, so that 3300° may be assumed as 
the effective temperature of these cMo stars. 

8. The same data may be found in an indirect way from eclipsing 
variables. If all the data about the eclipses, and hence the ratios of 
the stellar diameters to the orbit, are well known, the radial velocity 
gives the diameters in kilometers. Then we want the parallax to 
derive angular diameters, which, when combined with the observed 
radiation, give the temperatures. 

This method has been applied by S. Gaposchkin*s to a large num- 
ber of eclipsing variables. Because most of the parallaxes were highly 
uncertain, only a rough accordance with the scale usually adopted, 
with many large deviations, could be ascertained. Lately Pilowsky*4 
made another statistical application of this method. A number of 
eclipsing variables with known radial velocities gave stellar diam- 
eters, which were plotted as a function of the spectral class. From a 
number of binaries the Russell diagram was derived, giving absolute 
magnitude, i.e., total radiations, as a function of spectral class. 
Combining both curves, we find the temperature as a function of 
spectral class. The uncertainty here arises from the individual devia- 
tions from the statistical curves; in the first-named curve a deviation 
of 0.3 in log R occurs, which means a factor of 1.4 in T. In order to 
strengthen the correlation, the statistical relations of radius, abso- 
lute magnitude, and spectral class to mass were determined; but 
here the deviations from the curve are still larger. Though a general 
concordance with the usual scale of temperatures is arrived at, it is 
not likely that this statistical treatment will allow great accuracy. 
The physical correlation between spectral class and temperature is 
much narrower and more exact than the statistical correlation be- 
tween either of them and such data as mass, diameter, and absolute 
magnitude. So it will be better to look for separate stars for which all 
the data are well known. 

We will consider here only spectral class A, as a well-marked 


33 A.N., 248, 213, 1933. 34 Zs. f. Ap., 11, 265, 1936. 
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point in the scale of stellar temperatures. Among them we have 
only three stars with sufficiently reliable trigonometric parallaxes: 
B Persei, 07031+4; 8 Aurigae, 07037+ 4; and a Coronae, 0’053+ 10. 
For the first star the partial eclipse leaves the ratio of stellar to orbi- 
tal radius rather indeterminate. For the first and the third stars the 
velocity of the companion has not been observed; so the relative or- 
bit can be derived only from uncertain estimates of mass. Only 
6 Aur remains, for which the computation has been given already in 
Russell, Dugan, and Stewart’s textbook. For the two extreme sup- 
positions of uniform disk and of completely darkened limb, Shapley 
found r/a = 0.147 and 0.159. With m for each component, 2.07 + 
0.75 = 2.82, with the reduction to bolometric magnitude —o.24 
(for 10,000°, following Eddington), with the radial velocity of 220 
km/sec, and with the period of 3.960 days we find for 7 the values 
9630° and 10,130°. From the energy distribution of Table II and of 
Figure 2, however, the reduction to bolometric magnitude is found 
to be 0.23 mag. larger than for black-body radiation; taking this 
difference into account, the values for T are 10,300° and 10,740’. 
Hence, 10,500° may be assumed as the temperature of this Ao star. 

9. Thus, two points in the scale of stellar temperatures may be 
considered as well established: 10,500° for Ao and 3300° for cMo; 
we may estimate the limit of uncertainty of these values to be 
nearly 5 per cent of their value. They have been derived by direct 
methods of determining the effective temperature itself. The results 
from other methods which gave the first results for temperatures, 
from color and energy distribution, from ionization, and from the 
Boltzmann function are so strongly dependent on complicated con- 
ditions in stellar atmospheres that it is, as yet, not possible to derive 
reliable temperatures from them. On the other hand, the tempera- 
tures now established may enable us to make a decision between dif- 
ferent assumptions made in the discussion of energy gradients and 
line intensities. With the progress of the theoretical treatment of 
stellar atmospheres and their absorption coefficients, an increasing 
accordance may be expected. But this degree of accordance will 
serve rather as a test of the theories and as a source of knowledge 
concerning stellar atmospheric conditions than as a determination 


of temperature. 








A SPECTROSCOPIC EXAMINATION 
OF METEORITES* 


ARTHUR S. KING 


ABSTRACT 

Spectrograms of the material from thirteen meteorites, belonging to the nickel-iron, 
stony, and graphite types, were examined, and thirty-two of the constituent elements 
identified. Over half of these are not usually detected in chemical and mineralogical 
analyses of meteorites, some being among the rare terrestrial elements. Estimates of the 
abundance of the elements identified in iron and in stony meteorites, based on the degree 
of development of the spectra and the types of lines appearing, showed a fair agreement 
with the corresponding abundances found in chemical analyses by Noddack and Nod- 
dack. The distinctive differences in composition of iron and stony meteorites are dis- 
cussed, also the evidence that elements abundant in the stony type may be disintegra- 
tion products of those composing iron meteorites. The band spectra of meteoritic and 
of terrestrial graphites indicated no difference in the relative abundance of the.carbon 
isotopes. 

The composition of meteorites has been studied for the most part 
by the methods of chemical and mineralogical analysis used for ter- 
restrial minerals. As spectroscopic analysis has seldom been em- 
ployed, the following investigation was carried out to determine its 
usefulness in supplementing the other methods, as well as to study a 
collection of meteorites most of which had not been examined spec- 
troscopically. 

Earlier studies, involving in each case an examination of the spec- 
tra of a number of meteorites, have been made by Lockyer and his 
co-workers,’ also by Hartley and Ramage,’ Jewell,’ and Crookes.‘ 
These investigations showed the leading spectral characteristics of 
iron and stony meteorites (siderites and aerolites) and have been 
confirmed by my own work on the spectra of thirteen meteorites, 
including four irons, eight stones, and a graphite nodule. In addition 
to the identification of elements, I have attempted to indicate their 
relative abundance as shown by the strength of the more persistent 
lines and the degree of development of the spectra. 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 556. 

* The Meteoritic Hypothesis, 1890. 

2 Ap. J., 9, 221, 1899. 

3 [bid., p. 229. + Phil. Trans. R. Soc., A, 217, 411, 1918. 
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Specimens of the following meteorites, named after the locality of 
the fall, have been examined: 
Irons: Canyon Diablo, Ariz. 

Mount Joy, Pa. 
Staunton, Va. 
Xiquipilco, Mex. 

Stones: Pasamonte, N. Mex. 
Sioux County, Neb. 
Allegan, Mich. 
Ulysses, Kan. 
Gruver, Tex. 
Richardton, N.D. 
Weldona, Colo. 
Tryon, Neb. 

Graphite nodule: Elden, Ariz. 


The arc was regularly used as a light-source, the lower carbon elec- 
trode being charged with fragments or powder of the specimen. The 
spectrum of blank carbons, usually those used afterward with me- 
teorite material, was photographed on the same plate. Since an iron 
meteorite will conduct sufficiently to be used as an electrode, two 
pieces of Canyon Diablo iron were employed in this way; and in 
addition, clean drillings from this meteorite were used in the carbon 
arc. The Mount Joy iron was also used as electrode, against a carbon 
cathode. In examining the Elden graphite nodule, the arc was used 
with fragments as electrodes, and also with silver rods containing the 
material. Several spectrograms of the carbon-tube furnace contain- 
ing meteorite material were made for comparison with the arc spec- 
tra, the spectrum of the empty tube at the same temperature being 
taken first. As the lines of interest in this work are strong in all 
sources of moderate excitation, the arc served well for this survey. 

Since the arc spectrum of a meteorite contains thousands of lines 
in the region easily photographed, the spectrum of the meteorite 
was accompanied on the same plate by arc spectra of its most abun- 
dant constituents. From these spectra a great part of the meteoritic 
lines, perhaps gg per cent, could be identified by inspection. For iron 
meteorites, the arc spectra of iron, nickel, and cobalt were used for 
this purpose. For stony meteorites, arc spectra of iron, chromium, 
manganese, calcium, and magnesium were used. The meteoritic lines 
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not identified at once from these arcs either belonged to elements 
having relatively few lines, such as sodium, aluminium, and silicon, or 
were the most persistent lines of elements present in very small quan- 
tity. In either case they were readily identified from tables such as 
Kayser’s Hauptlinien. A final step was to look for traces of the 
ultimate and other highly persistent lines of elements suspected to be 
present but not noted in the first examination. The spectrograms 
were made in the second order of the 15-foot concave grating, except 
that for the yellow and red the first order was used. The exposures 
covered the range from \ 2500 to \ 6750. 
RESULTS 

Table I gives the elements identified, arranged according to atom- 
ic number, with estimates of abundance based on their spectra: “1”’ 
indicates that the most persistent lines were visible, and “5,” that 
the spectrum was developed so fully as to indicate a high content of 
the element. In assembling the data a tabulation of this sort was 
made for each meteorite, but among meteorites of the same type the 
variations for most elements were found to be small. The abundance 
given is that usually prevailing, although the range of differences is 
noted for some elements in the case of individual irons or stones. 
Even meteorites quite different in appearance and hardness, such as 
the Pasamonte and Weldona stones, showed very similar spectra, a 
fact probably to be explained by the presence in the two varieties of 
different compounds of the same elements. 

As is well known from other analyses, the typical iron meteorite is 
an iron-nickel alloy, although a number of the stronger lines are due 
to cobalt, copper, and silicon. The spectra of the irons showed the 
content of chromium and manganese to be very low as compared 
with that of such terrestrial irons as have been examined by the 
writer. The stones gave very intense spectra of magnesium and sodi- 
um, the D lines of the latter sometimes spreading over a range of 
60 A.s Calcium, silicon, chromium, and manganese showed rich 
spectra; while iron, though evidently much less abundant than in 
iron meteorites, was uniformly present in large quantity. In general, 

5 Sodium, according to chemical analysis, is much less abundant than magnesium; 


and this high intensity of sodium lines is probably due to the concentration of spectral 
energy in a few lines. 
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the more friable stones gave the stronger spectra of silicon, alumini- 
um, and calcium. 

The variation of nickel content in different stones is noteworthy. 
Two very friable stones, Pasamonte and Sioux County, showed only 


TABLE I 


ABUNDANCES OF ELEMENTS IN METEORITES 











ESTIMATED FROM SPECTRA | CHEMICAL ANALYsIS (I. AND W. NopDAck) 
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| 
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ELEMENT l 
Irons | Stones Irons Stones 
er eee ord kee | I 2 | — 5.0 X10° 
"Sa et eden eet | * | . | 3.9 X10-4 24. Mio 
Ce er ee I 5 _ | 7.18 X10 3 
Mes oa sancelnes 2 5 312 M10 | i .so- Mia 
Al... | 2-3 I-4 4.0 X105 L265 10 
Si | 3 3-5 8.0 X10 3 2. 542X106 * 
t hae I ? 1.47X10°3 5.00 X10"4 
| I Z _ 2.63 X10 3 
6 2 | 4 5.0 X10 4 I.92 X10? 
OME eo istie ci se - | I +- Et Jero-4 
Prater eu hnorcner I 2-3 4.0 X10°5 2.10 X10°3 
V.. | I | I 6.2 X10°° 3:0 Mia 
ce: I 4 2.4 X104 | 5.0 X103 
| Deer | I 4 3.0 X104 | 2.05 X103 
Fe 5 4 9.02X10% | 1.276X10° 
“(ihe Eee anae aeeer 3 | I-2 5.47X10 3 PSI. Ki10°% 
Ni oe cane 5 | I-3 8.46X107? 2.01 X10 3 
Cu. 3 2 | 3.05X10 4 E-e6- O10 
CE ere ee | I I— | L.15x10-4 a4, 10°" 
Ga | 2 I | L.3 MIC? + 
Ge 2 I 2.36X10 4 2.3% “16-5 
Rb. — I ° 4.5 M10" 
Sr I I + 9.2 Kios 
Mo I— — | 1.66X10°5 2:5 1e°° 
Ru I— _ 2.39X10 5 ro) 
Ag | I I 3.2 X10° fo) 
Sn I -- I.02X10 4 £0 “Xto0° 
ee ; I— | - | 2.0 X10° no: ie? 
S55 area ee ~ | I— | — to Mio 
Ba I I | a 2.0 X105 
| I— — ta SOR” ° 
ea rl I I— | 5.3 XI05 | 5-0 xX10°° 





* Usually masked, owing to the use of carbon electrodes. Exceptions noted in text. 


faint traces of the strongest nickel lines. The other stones, which, 
excepting Allegan, were of much harder structure, gave nickel lines 
of considerable strength. The iron spectrum was well developed 
for all. Crystallographic analysis would probably account for this 
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difference by the presence of non-nickel iron compounds, such as 
troilite, daubreelite, or lawrencite in the crumbly stones; while inclu- 
sions of some of the numerous nickel-iron compounds occur in the 
harder ones. If, however, we consider that the friable meteorites 
represent the most advanced stage of the stone-forming process, a 
possible connection with the hypothesis of G. N. Lewis® is suggested. 
According to this hypothesis, the composition of stony meteorites 
results from the atomic disintegration of the constituents of iron 
meteorites. Silicon, aluminium, calcium, and magnesium, known 
disintegration products of iron and nickel, are prominent in the stony 
type. If iron and nickel are equally subject to transmutation, the 
usual proportion in iron meteorites (about 11 to 1) would leave 
enough iron to give a strong spectrum after the nickel content had 
become very low. In fact, as Lewis notes, the splitting-off of an 
alpha particle from nickel should give iron and increase the disparity 
between iron and nickel in the stones. An atomic disintegration the- 
ory should account also for the large increase of manganese (not 
considered by Lewis) over that present in the irons. Evidence on 
this point is supplied through the transmuting, by Fermi and his co- 
workers,’ of cobalt into manganese by neutron bombardment. Co- 
balt, fairly abundant in the irons, is even harder to detect than nickel 
in the Pasamonte and Sioux County stones, while manganese is 
plentiful. 

However, the Allegan meteorite, a very friable stone, showed nickel 
and cobalt lines of fair strength, indicating a content of both, in 
agreement with the chemical analysis of this stone by Merrill and 
Stokes.* The evidence in favor of nuclear changes seems just suffi- 
cient to justify keeping the possibility in mind in future work with a 
large variety of stony meteorites. 

The Elden meteorite, found about 35 miles from the Canyon 
Diablo Meteor Crater in Arizona, is a rare specimen, being com- 
posed largely of carbon in the form of graphite (46.5 per cent’), and 
is the largest known of the type. The graphite forms a matrix 

6 Phys. Rev., 46, 897, 1934. 

7 Proc. Roy. Soc., A, 146, 483, 1934. 

8 Proc. Wash. Acad. Sci., 2, 41, 1900. 


9 From analysis by W. A. Sloane, U.S. Bureau of Mines. 
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crossed by numerous intersecting veins composed chiefly of nickel- 
iron. The meteorite was described by L. F. Brady,’ who suggested 
that it might have been an inclusion in iron belonging to the Canyon 
Diablo fall, the iron shell having fallen away. Evidence in confirma- 
tion of this explanation was given by R. E. S. Heineman," who 
found a very similar, though smaller, nodule imbedded in a section 
of Canyon Diablo iron. 

The spectrum showed strong CN bands in the arc and C 11 lines in 
the spark. Aside from this evidence of high carbon content, the 
Elden spectrum was mainly of iron and nickel and very similar to 
that of the Canyon Diablo iron.” Exceptional features, however, 
were stronger spectra of the alkali metals and earths than is usual 
with irons and evidence of more barium than in any other of the 
meteorites examined. 

The last two columns of Table I give the abundance found in chem- 
ical analyses of iron and stony meteorites by Ida and Walter Nod- 
dack."’ These results are averages for each type, made from mixtures 
of material from 16 irons and 42 stones. Iron particles were extracted 
magnetically as far as possible from the pulverized stone mixture, so 
that the iron content for the stones is low. A plus sign indicates that 
the element was detected, sometimes by spectroscopic means, but 
the amount was not measurable. 

Table I shows a fairly close agreement between the spectroscopic 
and chemical abundance for both iron and stony meteorites. Ele- 
ments found by Noddack and Noddack to be present in their mix- 
tures in amounts greater than a few parts in a million, but not listed 
in my table, are gases, such as oxygen and chlorine, or metals whose 

10 Amer. Jour. Sci., 21, 173, 1931. 


1 Tbid., 23, 417, 1932. 

2 The writer is indebted to Professor F. A. Jenkins of the University of California for 
a comparison of the relative abundance of the two carbon isotopes, C” and C3, in the 
Elden meteorite and in terrestrial graphite. The bands at \ 4737 and A 4744.5, due 
to C?C and CC, respectively, were produced in a discharge tube of the form pre- 
viously used in the study of isotope bands (Wooldridge and Jenkins, Phys. Rev., 49, 
404, 1936). Microphotometer measures of the bands given by the two graphites showed 
definitely that there is no difference in the relative intensity greater than the error of 
measurement. For details of this experiment see Jenkins and King, Pub. A.S.P., 48, 
December, 1936. 


3 Naturwissenschaften, 18, 757, 1930. 
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persistent lines are in an inaccessible spectral region. A fundamental 
limitation to the detection of very small percentages by the spectro- 
scopic method must be recognized. The most persistent lines are 
those produced in transitions from the ground level in the atom. If 
the energy involved in such a transition is concentrated in one line, 
or at most a very few lines, the spectroscopic sensitivity may be 
very great. But when this energy is divided among a considerable 
number of lines, perhaps several multiplets, as is usual in the com- 
plex spectra, the possibility of detecting very smal]l amounts be- 
comes correspondingly less. The few elements not detected in my 
spectra, but for which Noddack and Noddack obtained measure- 
ments, are of this kind. Examples are the rare earths cerium, neodym- 
ium, and samarium, measured as three to four parts in a million, 
and the platinum metals, also of very low content. The latter have 
been the special object of chemical analyses of the Canyon Diablo 
meteorite," and evidence of the whole platinum group has been 
found. Ruthenium only was identified in my spectra, by means of a 
very faint line. Noddack and Noddack give it the highest abundance 
in the group, by a small margin; but it may be noted that other 
chemical analyses do not agree on this point. In the case of plati- 
num, the strongest line, \ 3064.71, would be effectively concealed by 
a strong nickel line. 

The conditions of my experiments were unfavorable for the detec- 
tion of either oxygen or carbon. Oxygen is one of the chief constitu- 
ents of stony meteorites; but oxide bands in the open arc, as well as 
the infrared oxygen lines, could as well have come from the air. Car- 
bon, given an abundance of about 4 X 10 4 by Noddack and 
Noddack in both irons and stones, and found in meteorites as dia- 
mond and graphite, would regularly have been masked through the 
use of carbon electrodes. The Canyon Diablo meteorite, when used 
as electrodes, showed the CN bands faintly. The strength of the car- 
bon spectrum from the Elden nodule has been mentioned. Sulphur, 
another prominent meteoritic substance, is difficult in arc spectra on 
account of the high excitation potential of lines in the accessible 
region. 


4 See summary by J. L. Howe, Science, 66, 220, 1927. 








514 ARTHUR S§. KING 


The grounds for identification of some of the elements present in 
small amounts may now be discussed. 

The presence of /ithium in meteorites has often been left uncertain 
by analysis, though in stones it would be expected to accompany the 
large content of sodium. The ultimate line \ 6708 could be seen in 
the irons and showed considerable strength in the stones, as did 
A 6103.6. 

Very faint lines of phosphorus were measured at AX 2533.99, 

553-27, and 2554.92 in the Canyon Diablo iron. Spectra of the 
stones in this region may not have been sufficiently strong. 

Scandium was identified in stones by means of faint lines at 
d 3911.81 and A 4023.69. 

The strongest lines of t2tanium and vanadium could be seen in both 
irons and stones, while in the stones a large number of titanium lines 
were identified. 

Zinc gave \ 4722 and J 4810 distinctly in the Canyon Diablo and 
Xiquipilco irons. 

Gallium, found in very small amounts in terrestrial minerals, 
though rather widely distributed in igneous rocks, showed \ 4033.01 
and 4172.05. The first, being only 0.06 A from an ultimate line 
of manganese, was masked in the stones, but was easily detected in 
the irons on account of the faintness of the manganese line. The sec- 
ond, A 4172, is 0.08 A from an iron line of medium strength but was 
usually resolved in both types. Gallium, though seldom detected in 
meteorites by other methods, was identified in the early spectro- 
scopic investigations; and its presence in meteorites has been taken 
as evidence of its probable presence in the sun," the gallium lines in 
the solar spectrum being at that time difficult of detection on ac- 
count of blends. 

Germanium, another rare terrestrial element, showed its ten 
strongest lines from \ 2592 to A 4686 distinctly in the irons. They 
were weaker in the harder stones, such as Tryon and Weldona, and 
very faint in the soft stones, Pasamonte and Sioux County. Several 
of these lines are found in the solar spectrum. Goldschmidt” ques- 
tions the extreme scarcity of germanium in terrestrial rocks on the 

1s Hartley and Ramage, Ap. J., 9, 214, 1899. 

6 Proc. R. Inst., 26, 73, 1929. 
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ground that it is very generally associated with aluminium and in 
chemical reactions is likely to be masked by the latter. 

Rubidium gave \ 7800 and \ 7947 faintly in the Pasamonte stone 
and still weaker in the Elden graphite. These lines did not appear in 
the Canyon Diablo iron. 

Strontium could, as a rule, be detected by \ 4607. A spark spec- 
trum of the Elden graphite, however, gave \ 4077 and Xd 4215 of 
Sr 1 with considerable strength. 

Molybdenum and ruthenium were identified in the Canyon Diablo 
iron by A 3193.98 and A 3498.95, respectively. The latter was very 
faint but measured in close agreement with the accepted wave- 
length. 

The ultimate lines of silver, \ 3281 and X 3383, could, asa rule, be 
seen in both irons and stones. These lines were sometimes strong in 
Canyon Diablo spectrograms, but varied in strength according to the 
part of the specimen examined. 

Tin lines were very faint and often blended with those of other 
elements; AX 3034.12, 3177.05, and 3263.33 were identified. 

Antimony was recognized by \ 3232.52 only. 

The ultimate line of caesium, \ 8521.15, could be seen very faintly 
in the Pasamonte stone and in the Elden graphite, not in the Canyon 
Diablo iron. Other speciments were not examined in this region. 

The barium line X 5535.53 was of fair strength in both irons and 


x. 
stones. Sometimes other 


Je 
lines appeared, especially in the Elden 
graphite. 

Gold showed X 3675.95 very faintly on two Canyon Diablo spec- 
trograms. 

The lead lines AA 3639.6, 3683.5, and 4057.8 could usually be seen, 
sometimes also \ 2833.1. They were stronger in the irons and quite 
variable on different spectrograms. 


CONCLUDING REMARKS 


Leaving out of account a search for elements present in small quan- 
tity as exhaustive as that of Noddack and Noddack, we find that 
fully half of the elements listed in Table I, as identified by their spec- 
tra, are not, as a rule, detected in chemical and mineralogical analy- 
ses of meteorites. A quantitative spectrum analysis would have as 
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its main object the meteoritic content of these less abundant ele- 
ments. For this purpose a very smal] amount of the material would 
be sufficient; and with standardized samples of various elements at 
hand and a routine established, the analysis could be made in much 
less time than by the other methods. For the production of certain 
spectra special experimental arrangements would be necessary. How 
these would compare in difficulty with the chemical operations re- 
quired for some groups of elements, the writer is not prepared to say. 

Comparisons of abundances of the elements in meteorites with 
those in terrestrial minerals have frequently been made.’ Such a 
comparison by Noddack and Noddack with the results of their analy- 
ses of sixteen hundred minerals showed a relatively higher content 
of magnesium, iron, and nickel in meteorites than in the earth’s 
crust. Apart from these elements, perhaps the most noteworthy dif- 
ference is a considerably larger amount of gallium and germanium in 
meteorites. 

Meteoritic elements which are really abundant are well within the 
first third of the atomic-number series, nickel, No. 28, being the last 
plentiful element. A similar condition holds for elements in the 
earth’s crust, according to the data of Washington and Clarke;® but 
this predominance of the lighter elements is shown in the paper of 
Noddack and Noddack to be more pronounced for meteorites, owing 
chiefly to a relatively higher terrestrial abundance of strontium, 
barium, zirconium, tungsten, and the rare earths. In the case of 
compounds, the differences are more conspicuous, some constituents 
of meteorites not being found among terrestrial minerals. 


For their co-operation in supplying specimens of meteorites, I am 
indebted to Professor H. H. Nininger, of Denver; Professor F. C. 
Leonard, of the University of California at Los Angeles; Dr. C. O. 
Lampland, of the Lowell Observatory; and Mr. L. H. Brady, of the 
Museum of Northern Arizona, Flagstaff. 

CARNEGIE INSTITUTION OF WASHINGTON 


Mount WILSON OBSERVATORY 
August 1936 


17 See, for example, G. P. Merrill, Proc. Amer. Phil. Soc., 65, 119, 1926; Farrington, 
Meteorites (Chicago, 1915), p. 215; I. and W. Noddack, op. cit. 
8 U.S. Geol. Surv. Prof. Paper, No. 127. 

















EFFECTS OF RED SHIFTS ON THE DISTRIBUTION 
OF NEBULAE 


EDWIN HUBBLE* 


ABSTRACT 


Two new surveys of nebulae, to mpg= 18.47 and 21.03, are described which, together 
with three earlier surveys to mpg=109, 19.4, and 20, indicate the apparent distribution 
in depth. The numbers of nebulae per square degree, NV, to limiting magnitudes m, are 
represented by the relation - 


log N=0.6 (m—Am) —9.052+0.005 , 


where 
log Am=o0.2 (m—Am)—4.239+0.008 . 


The corrections, Am, are interpreted as apparent departures from large-scale uniform 
distribution due to the effect of red shifts on apparent luminosities. 

The effects to be expected are calculated on the two assumptions that red shifts are 
(a) velocity shifts and (6) not velocity shifts. With an average T>=6000° assumed for 
integrated nebular radiation, Am=4.0 dd/X (velocity shifts) or 3.0 d\/A (not velocity 
shifts). The empirical expression, derived from the apparent distribution, together 
with the velocity-magnitude relation, is Am= 2.94 dd/X. 

If red shifts are not velocity shifts, the apparent distribution agrees with that in an 
Einstein static model of the universe or an expanding homogeneous model with an in- 
appreciable rate of expansion, provided spatial curvature is negligible. 

If red shifts are velocity shifts which measure the rate of expansion, the expanding 
models are definitely inconsistent with the observations unless a large positive curvature 
(small, closed universe) is postulated. The maximum value of the present radius of 
curvature would be of the order of 4.7 X 108 light-years; and the mean density, of the 
general order of 10°”. The high density suggests that the expanding models are a forced 
interpretation of the observational results. 


Surveys with large reflectors on Mount Wilson and Mount 
Hamilton show that the distribution of nebulae’ over the sky down 
to limiting magnitudes? between 19 and 20 is moderately uniform 
except as affected by local obscuration within the galactic system 
and by the occasional great clusters. The distribution in depth, as 
indicated by the rate of increase of the numbers of nebulae to suc- 
cessive limiting magnitudes, also appears to be approximately uni- 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 557. 

1 Mt. W. Contr., No. 485; Ap. J., 79, 9, 1934; Lick Obs. Bull., 16,177 (No. 458), 1934. 

2 Unless otherwise stated, the terms “nebula” and “magnitude” in the present dis- 
cussion refer to extragalactic nebulae and to photographic magnitudes on the inter- 
national scale, respectively. 
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form when tentative corrections are applied for the effect of red 
shifts on apparent luminosity. The second conclusion is less reliable 
than the first and requires further investigation before it can be used 
with confidence in cosmological theory. The following contribution 
is a re-examination of the subject, carried out in greater detail and 
with the aid of additional material. 

The observational data are five homogeneous groups of nebular 
counts, each giving the number of nebulae per unit area to a particu- 
lar limiting magnitude. The five pairs of values determine an 
empirical relation, log V,=/(m), which diverges systematically from 
the theoretical relation representing uniform distribution. The de- 
partures from uniformity, expressed as corrections to the observed 
magnitudes, are then formulated as a relation Am =f(m), which sums 
up the observational results of the investigation. This relation in- 
cludes the effect of red shifts, in addition to observational errors, 
space absorption, and real departures from uniform distribution; 
hence its interpretation is of critical importance for cosmological 
theory. 

I. OBSERVATIONAL DATA 
PREVIOUS SURVEYS 

Three groups of nebular counts are available from previous 
surveys. These represent the 1-hour exposures on Eastman 4o 
plates with the too-inch and 60-inch reflectors already discussed 
by the writer,’ and the 1-hour exposures on Imperial Eclipse plates 
with the 36-inch Crossley reflector discussed by Mayall.‘ Since the 
details for individual plates are already published, they will not be 
repeated here. The limiting magnitudes were 20.0, 19.4, and 19.0, 
with uncertainties of the estimated order of 0.1 mag. The mean 
values of log N per plate were derived from all available material 
in the polar caps (8> 40°), reduced to standard exposures of 1 hour 
by corrections derived from empirical relations between log N and 
log E. The results were reduced to mean log N per square degree 
with the aid of corrections for coma and for plate area. 

For the present investigation the data were used as follows: The 
1-hour exposures with the 60-inch and 1oo-inch reflectors (214 and 

3 Mt. W. Contr., No. 485; Ap. J., 79, 8, 1934. 

4 Lick Obs. Bull., 16, 177 (No. 458), 1934. 
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228 plates, respectively) were extracted from the Mount Wilson 
surveys and reduced separately. Mayall’s values, representing all 
the 36-inch material, were adopted as published, except that lati- 
tude corrections were applied in order to make the data compa- 
rable with the Mount Wilson counts. These corrections’ increased 
the mean log N by the increment 0.032 without materially affecting 
the dispersion. The results are included in Table IV. 
NEW SURVEYS 

The two additional groups of counts represent 2-hour exposures® 
with the 1oo-inch reflector and 20-minute exposures with the 60- 
inch, both on Eastman 4o plates. The fields, as before, are re- 
stricted to the polar caps; and the counts are reduced to log N per 
square degree in the same manner as in the earlier surveys. Each 
plate was examined three times, the last examination being a con- 
tinuous review of each group as a unit. It seems probable that a 
more critical standard was attained in the identification of threshold 
nebulae. This advantage, combined with the appreciably higher 
sensitivity of recent emulsions, led to relatively fainter limits, and 
hence larger numbers of nebulae, than would be expected for cor- 
responding exposures in the former surveys. Details for individual 
plates are listed in Tables I and II. 

LIMITING MAGNITUDES 

The limiting magnitudes for the earlier surveys were given to the 
nearest 0.1 mag., and the published values are adopted without 
change. Since the new surveys represent extreme cases and play 
important réles in the correlation of log NV,, and m, the limits are 
given to another decimal place. The limiting magnitude of the 
counts on the 20-minute exposures was derived from about fifteen 
plates of six fields in which nebular magnitudes down to m=19.0 

$s Mayall concluded that “some soit of latitude effect is probably involved in the 
Crossley material” but that “the data are insufficient to establish the nature of the 
variations.” The more extensive material in the Mount Wilson surveys appears to 
establish a cosecant law. 

6 Thirteen of the exposures with the 1oo-inch are not precisely 2 hours but range 
from 100 minutes to 150 minutes. The corrections reducing log N to the standard ex- 
posure range from 0.02 to 0.13 and average 0.066. The uncertainties in these small cor- 
rections are believed to be negligible. 
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TABLE I 
SURVEY FIELDS 
100-Inch Reflector; 2-Hour Exposures 
































| | 
B N | E* | Z | Q Ni | log Mt 
+73° 276° | 135 | 20° FG 157 2.33 
63. Oo | 150 40 GE | 263 | 4.54 
60. 280 | 150 36 | GE | 250 | 2.37 
S23. 17 | 120 19 GE | 238 2.47 
52. 14 110 31 G 121 | 2.30 
| | | | 
51 | 22 || 120 28 E 328 «=| = (2.58 
4034.3 30 120 36 G m7 2.44 
eee 15 120 36 GE 216 | 2.47 
Pl, eer 159 135 20 FG 1446 | 2.35 
eS ee 164 105 19 | GE 203 | 2.50 
Pe 329 110 45 G 142 | 2.43 
ee ee, 20 100 34 G 185 2.58 
+40... 143 120 35 GE 23 2.53 
—40....... 50 130 24 G 200 2.45 
ar | 57 120 34 G 157 2.41 
ee 7 130 25 | E | 214 2.38 
OO ee Reernee? 130 120 22 GE 190 2.41 
Be iss ase 31 IIs | 43 | FG 160 2.53 
Pi nl 30 120) | 38 GE 222 2.51 
Pee 32 120 | 34 G 116 | 2.27 
ne eee ‘I 61 120 28 FG ya |  @.et 
eee 34 120 43 | G 154 2.42 
45.----- 45 120 32 3 144 2.35 
ee 160 120 | 41 G 175 2.46 
Sy ee 79 120) | 21 G ro | 2.35 
| 

ee a 104 120 30 | GI 241 a8 
Oy Se 102 120) | 33 GI 182 2.30 
oe 38 115 | 47 G 185 2.52 
BOG aces 50 r0~—C i} 30 GI 192 2.40 
BOs ocx 70 120 32 FG 198 2.56 
Oa gets oa’ 120 120 31 FG 172 2.50 
roe 34 110 4! GE 132 2.30 
ss. 100 120 28 GE 223 2.44 
Se. 140 120 37 G 210 | 2.49 
58. 82 120 34 E 204 «| 2.54 
§8........ 123 135 39 F 145 | 2.45 
Oo aia a s0 3 89 120 37 GE 246 2.49 
Gao cecok 50 120 43 G 171 | 2.41 
ae 110 120 35 GE 200 | 2.52 
eae 95 120 36 FG 231 | 2.60 
ee 40 135 51 GE 301 | 2.56 

















* F =Exposure-time in minutes. 
t log N is reduced to the standard exposure, E =120, by the corrections A log N =1.33 log E. (Mt. W 
Contr., No. 485; Ap. J., '79, 8, 1934). 














EFFECTS OF RED SHIFTS 521 
TABLE II 
SURVEY FIELDS 
60-Inch Reflector; 20-Minute Exposures 
B » Z Q Vr | log N]] B N z Q Nr |logN 
Ee S z acme: ; a nae Cee Siete: 
75 320°| 32 |1 18 | 1.30||—43 61°| 42°) GI 30 | 1.65 
74 ) aye 28 | GE 10 | 1.07|] 43 66 34 | E 16 | 1.31 
73 313 25 | E 18 | 1.20|| 44 4! 34 | GE ¥2.1°323 
63 323 24 | GE 24 | 1.46]| 44 55 38 | E Ig | 1.40 
63 325 24 | E 15 | 1.22) 44 66 40 | I 54) 2237 
} 

62 338 28 | GI 23. | 1.44l| 44 80 20 | I 46 | 1.74 
62.. 340 33 E 16:| .25]} 45 40 30 GE 18 | 3.41 
60.... 20 31 | E 10 | 1.05]} 45 55 39 | E 8 | 1.01 
60.. 135 3216 14 | 1.30]| 45 80 24 | E 28 | 1.53 
60 136 aq) G 17 | 1.30|| 45 102 28 | I 38 | 1.66 
60.. 138 42 I 6 | 0.96]] 45 109 33 | E #7 | 2.92 
60 170 38 G 5 | 0.871] 47 79 29 | I ' 27 | 1.51 
59 24 30 | E 21 | 1.37|| 47 102 25 | GE 14 | 1.27 
59 28 37 | E 21 | 1.38]] 49 it ae 38 | GE 7 a ee 
SB... 34 42 | GE 6 | 0.91]| 49 39 35 | GE 21 | 1.45 
58 36 31 | GI 13, | I.21]| 49 61 31 E 37 | 1.65 
58 40 38 GE 5 | 0.82]| 490 99 32 E 52 | 1.80 
5s 21 38 | E 13,| 1.19|| 49 149 39 | E 8 | 1.00 
54 | 110 | 37 | GE 2 I.50]| 50 60 | 35 | GE| 28] 1.58 
ae, | 111 39 | GE| 23] 1.490]} 5° 75 24 | GE 6 | 0.89 
54.. 193 25 | GE 25 | 1.50]| 50° 81 32 | E 8 | 0.98 
Ba. 195 23 | GE 40 | 1.69|| 50 99 35 | E tz } ¥-20 
53 196 21 | GE 8 | o.g9]] 50 150} 39| 1 18 | 1.36 
53 225 40 | GE 20 | 1.43i| 5 75 24 |G 16 | 1.36 
Bae. 294 | 42 | GE tp 1 2.37 52 81 38 | GE 19 | 1.4! 
er 247 43 | GE 10 | 1.14|| 51 89 28 | E 7 | 0.91 
BOW. 05} 540 31 | G 6 | 0.96]} 51 160 44|E r5 | 1-20 
50. I4I 34} G 8 | 1.09]] 51 162 45 | E 10.1 E.352 
46 133 28 | E 23 | 1.44]] 52 81 42|E ¥4 | 3:95 
462-3) ESA 30 | GE 27 | 1.56]| 52 89 32 | 37 | 3-30 
46. 136 34 | E 26] 1.52]| 52 163 48 | E 14 | 1.28 
42. oc2sh Bae 42 | GE 8] 1.07]] 53 52 32 | E I9 | 1.35 
ae 123 45 | GE 4|0.78]] 5 100 26 | E 19 | 1.34 
4I 121 39 | GE 7 | 1.02|| 53 110 27 | E 29 | 1.52 
AO cvs 2 30 | G 16 | 1.4! 54 36 42 | GE rr | 1.18 
40 ..| 142] 34] FG 0.90] 54 51 34 | E 13 | 1.19 
+40 a) 344 39 | G 26 | 1.65]} 54. 100 29 | E 42 1 2.590 
—ql.... 171 47 | E go | 1.11 54 III 28 | E ka he ee 
Pe) 36 | GE 18 | 1.41 54 129 34 | E 16 | 1.28 
—=A2....5-) 372 47) E + | 2.001—SS5.. 35 45 | E Ky lee es 
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TABLE Il—Continued 

















] | | 
B nN £ Q Nx | log N || B [> 2 ll Q Ny | log N 
an of Sed |e y A ADL, | ase il aes 
—55' 50] 37°| E 26 | 1.49||—6r......| 70° | 42°| GI 10 | 1.12 
55 73 29 | GE OX 13.488 162....-s] OF 1 3? a ae (es i 
55 77 33 | E Of e.008 102-55) cose) “aks 41 | 1.69 
55 100 ay © 22 | 1.40|] 63... 99 | 36 | I 20 | 1.36 
55 III 30 | E 27 | 240 (63.... 118 | 41 | I a3 I t.42 
55 125 49 | I 22 | 1.40|] 64 100 | 36] E 5 ae a ae 
55 126 34 | E 38 | <-06 -64......] 2to.| “aE 1B 18 | 1.34 
55 130 | 36] E 17 | 1.30]| 65 | 100 | 38] E 34 | 1.60 
ee 2] 39| GE 9] 1.07|| 65 | 103 2|/E 25 | 1.48 
56 73 | 321E 55 | 1.80]] 65......] 120 | 40 | E 21 | 1.39 
56 77 31 | E 18 | 1.32|| 66 85 | 40|E 16 | 1.26 
56 127 36 | E 37 | 1.64]] 66 | 104 39 | E 43 | 1.69 
<7 40 40 | E 7 11.34|| 67 | 85 | 38 | E ¥3 1 2.27 
57 41 42 | GI 1] 1.17]| 68 64 43 | E 26 | 1.49 
57 60} 36|G g | 1.12|| 69 | 63 | 45 | E 16 | 1.28 
57 123 33 | GE 31 | I.59]|—70. 62 | 48] E 7 o ee 
60 | 39 | GE 16 | 4.32 | 
_ | | 
58 123 35 | E 26 | 1.49 | 
59..--.-| 51] 40] E 18 | 1.33 
BO cy LES 39 | E 18 | 1.33 | 
59 123 | 38| E 18 | 1.33] | 
59 150 | 44] GE 13] 1.24 
59 152 45 | GE 24) 1.51 | 
60 50 44 | E i} £27 
—60 71 38 | GE 26 | 1.53 



































or 19.2 have been measured from schraffierkassette images. The 
scale was determined from the North Polar Sequence and three 
Selected Areas in which the scale has been extended by Baade 
beyond the reliable ranges in the Mount Wilson Catalogue of Photo- 
graphic Magnitudes.’ The fifteen plates were scattered at intervals 
throughout the period of time required to complete the survey. 
Since the total range in the limits of the individual plates is only 
0.25 mag., the mean value, 18.47+0.03, is believed to be reliable 
except for uncertainties in the scale of the comparison stars and in 
the zero point of the scale of nebular magnitudes. 

The estimation of the limiting magnitude for 2-hour exposures 
necessarily involved considerable extrapolation. A dozen 20-minute 
exposures with the 1o0o-inch telescope, scattered throughout the 


7 Seares and others, Carnegie Institution of Washington Publication, No. 402, 1930. 
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group, were centered on the same fields as were used in calibrating 
the 20-minute exposures with the 60-inch. The mean limiting magni- 
tude of these plates, reduced to the zenith and to excellent quality, 
is 19.11+0.02, very close to the limit of the reliable nebular magni- 
tudes determined with the schraffierkassette. Hence, for 2-hour ex- 
posures, the limiting magnitude should be 


19.11 + 2.5 p log — = 19.11 + 1.94), 
where p is the Schwarzschild exponent in the reciprocity law. 

The value of p for surface images of equal density on fully de- 
veloped Eastman 4o plates increases as the density diminishes and 
approaches unity at the threshold itself. The value unity was used 
for deriving magnitudes in the previous surveys, and observational 
data were presented which justified the procedure. A re-examina- 
tion of the material with special reference to very small surface 
images suggests values of p ranging from o.9g5 to 1.0 for the particu- 
lar conditions under discussion, and indicates that the mean should 
give the desired results to a close approximation. The limiting mag- 
nitude for the 2-hour exposures is then 21.00+0.04, with uncer- 
tainties of much the same order as those for the nebular magnitudes 
around m=1g in the standard fields. 

A second method of calibration, also applied in the previous sur- 
veys, depends upon a known relation between diameters and surface 
brightness of threshold images.’ Several 2-hour exposures on Se- 
lected Areas, in which Baade has extended the scale of stellar magni- 
tudes, were made under conditions of poor seeing. Stars with 
known magnitudes were thus smeared out into surface images simi- 
lar to, although considerably larger than, nebular images at the 
threshold of identification on plates of excellent quality. The same 
sorts of results are obtained with small extra-focal images, although 
the shadow of the Newtonian flat introduces uncertainties. The 
relative diameters of the threshold images of these stars and the 
threshold images of nebulae on good-seeing plates with equal ex- 
posures indicate the relative surface brightness and hence the differ- 
ence in the integrated magnitudes. For instance, a 2-hour exposure 


8 Mt. W. Contr., No. 453; Ap. J., 76, 106, 1932. 
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of S.A. 68, with diameters of faint images about 8"5, registered 
stars down to about 19.9. Diameters of threshold nebulae on the 
sharpest survey plates average about 3‘2. From the relative areas, 
the nebular images should be 2.12 mag. fainter than the stellar 
images, but the correction for surface brightness as a function of 
diameter reduces the increment to 1.15 mag.; hence the threshold 
magnitude for the nebulae is 21.05. 

Four 2-hour exposures are available and, in addition, six 1-hour 
exposures and five schraffierkassette plates (images 16” X 16”) with 


TABLE III 
LIMITING MAGNITUDE FOR 120-MINUTE EXPOSURES 














| | 
Me 
MEAN | No. oF MEAN | RE - 
Exp. PLATES | Diam. | | err 
Observed Reduced* 
| 
Min | | | 
<a AOE ree |} 12 | 3%2 | 19.11t0.02 | 21.00+0.02 | Good seeing 
en 6 | 6 | 19-57 06 | 21.02 os5 | Poor seeing 
<< ee e | ¥ | 20.20 06 | 21.02 06 | Poor seeing 
oa ; 4 5 | 16 | 19.08+0.06 21.08+0.06 | Schraffierkassette 
Mean. | | Oe | 21.03 +0.03 | 





* The reduced values of mg refer to nebulae with diameters of 3"2 at the threshold of identification on 
120-minute exposures at the zenith, made with the 1oo-inch reflector on Eastman 40 plates. Exposures with 


a 
+ 


the schraffierkassette range from 110 to 180 minutes, E=145, log E =2.164. Mean diameters are approxi- 
mate; the reduced limiting magnitudes are the means of those for the individual plates in each group. 
exposures ranging from 2 to 3 hours. Results from the latter two 
groups involve extrapolations in exposure times as well as in di- 
ameters. Limiting magnitudes from the three groups, together 
with the limit from the 20-minute exposures, are listed in Table 
III. The limiting magnitudes, m,, refer to nebulae with diameters of 
3.2 at the threshold of identification on 2-hour exposures at the 
zenith. The simple mean, 21.03+0.03, is adopted as the most prob- 
able value. The small probable error, of course, merely indicates 
the consistency of the data and does not fully eliminate the possi- 
bility of systematic errors. 

_ The series of five limiting magnitudes form a consistent system in 
which the scale appears to approximate closely the conventional 
scale over the range involved, and the accidental errors are less than 
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o.1 mag. This conclusion is drawn from examination of the Mount 
Wilson data alone, but Mayall’s data for the Crossley counts are 
thoroughly consistent and may be included in the general statement. 

The zero point of the scale is uncertain by the amount of the 
errors in the standard nebular magnitudes between 18.5 and 19.0. 
These arise from two sources. One, the uncertainties in the scale of 
the standard sequences of stellar magnitudes, is believed to be small 
since it involves the North Polar Sequence and three independent 
extensions of that scale from about m=17.0 downward. The other 
arises from the method of measuring nebular images and involves 
the question of the relative amount of luminosity contributed by 
extremely faint but relatively large outer regions of nebulae. A 
rough rule is that reliable total magnitudes may be derived by 
smearing the nebulae and comparison stars over areas with di- 
ameters three times the diameters of the nebulae, as estimated by 
simple inspection of focal plates.? These conditions were fulfilled in 
measuring the very faint nebulae in the standard regions, and hence 
the uncertainties arising from outer regions are believed to be less 
than o.1 mag. 

For many purposes, including investigations of distribution and 
associated problems, the zero point is relatively unimportant; the 
significant feature is the scale. The present series of magnitudes can 
be represented to the nearest o.1 mag. by extrapolations from m= 
19.0, based on the reciprocity law with an exponent p between 0.95 
and 1.0. The reliability of the scale thus depends largely on the value 
of p. Since it is improbable that / should be greater than unity, any 
appreciable revision would probably reduce the range in the ob- 
served magnitudes and make the faint magnitudes brighter with 
respect to the bright magnitudes. 

SIGNIFICANCE OF LOG V 

The mean values, log WV and log N, for the five surveys are listed 
in Table IV, together with the dispersions, o, calculated from the 
residuals relative to log NV. The three quantities should satisfy the 
relation 

log VN = log N+ 1.15107, 
9 The rule is based on a comparison of magnitudes from extrafocal and schraffierkas- 


sette images with those derived by Stebbins and Whitford from accurate measures with 
photoelectric cells, made directly at the foci of reflecting telescopes. 
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provided that the frequency distributions of log NV are normal error- 
curves; the small differences, C—O, show that this condition is very 
closely approximated. 

The surveys represent random sampling with average samples of 
five different sizes. The significance of the dispersions is somewhat 
obscured by the fact that the sizes of samples do not correspond to 
the adopted log NV but represent the numbers of nebulae actually 
identified on the plates before the counts were reduced to standard 
conditions. The latter data are included in the table and show that 
the dispersion tends to diminish as the size of the sample increases. 


TABLE IV 
DATA FROM INDIVIDUAL SURVEYS* 


} = log V 


l 
| | | 
Tel. | £ | we | WM m flog | o | log ¥ 
| | | (obs) 
ane | Watch cig ee De a - 
60-in. 20 | I2r | 19 | 18.47 | 1.829 | 0.224 | 1.887 | 1.885 
36-in. ; | 60 | 284 | 35 | 19.0 | 2.103 | .225 2.161 | 
60-in. ; 60 | 214 | 4t i io.4 | 2.312] 156 | 2.340 | 2.342 
100-iN. . | 60 | 228 | 52 | 20.0 | 2.650 | .176 | 2.686 | 2.685 
100-in. |} 120 | 4I | 200 21.03 | 3.154 | 0.084 | 3.162 3.162 
| | | 





* Tel. =telescope; E =exposure time in minutes; VP =number of plates; VW: =average number of nebu- 

lae actually identified per plate; m=limiting magnitude of survey; Jog ¥ =mean logarithm of number of 
nebulae per square degree; o =dispersion in log NV; log W =log V+1.15107; log W (obs) =logarithm of the 
mean number of nebulae per square degree (each plate reduced to standard conditions) 
This tendency, together with the close agreement of the results for 
the two galactic hemispheres and the absence of conspicuous 
systematic variations in latitude or longitude, emphasizes once more 
the random character of the large-scale distribution of nebulae 
throughout the sample of the universe available to inspection. 

While the large-scale distribution appears to be essentially uni- 
form, the small-scale distribution is very appreciably influenced by 
the well-known tendency toward clustering. The phenomena might 
be roughly represented by an originally uniform distribution from 
which the nebulae have tended to gather about various points until 
now they are found in all stages from random scattering, through 
groups of various sizes, up to the occasional great clusters.’? The 

© This representation is purely formal and has no genetic implications. For pur- 
poses of speculation, the reverse development seems preferable, namely, that nebulae 
were formed in great clusters whose gradual evaporation has populated the general 


field. 
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tendency appears to operate on a relatively modest scale—at least 
no clusters are known with more than a few hundred members; 
hence irregularities tend to average out among samples that are 
large compared to a single cluster.’ In small samples irregularities 
are conspicuous, and, as an empirical fact, they lead to a random 
distribution of log N rather than of NV. Accidental errors of observa- 
tion and reduction also tend to introduce random scattering of per- 
centage deviations (since they are generally equivalent to variations 
in limiting magnitude), but the effects are relatively small and are 
presumably constant for all the surveys. 

The significant quantity, for the present purpose, is the arith- 
metical mean, N, related to the geometrical mean by the expression 
already given, in which o represents the combined effects of irregu- 
larities and accidental errors. If the errors are small compared to 
the irregularities, as in the surveys using small samples, the relation 
may be used as it stands, or N may be calculated directly from the 
data. On the other hand, if the samples are so large that the dis- 
persion is dominated by errors rather than by irregularities, the 
term including o may be ignored, and WN derived from the simple 


relation 


log N (distribution) = log NV (observed) . 


The dispersion arising from accidental errors alone is probably 
of the order of 0.06 in log N (0.1 mag) and about equal to that aris- 
ing from irregularities in the survey using the largest samples. Ap- 
propriate corrections for accidental errors would consist of a con- 
stant term, 1.151 (0.06)?=0.004, to be subtracted from each log N 
in Table IV. However, the term is much smaller than the various 
probable errors, and it may be neglected or thrown in with the un- 
certainties of the zero point of the magnitude scale. It appears, 
therefore, that the tabulated log NV may be used directly for com- 
paring the numbers of nebulae to the successive limits of the various 
surveys. 

t Samples in the regions of great clusters were rejected in the present surveys; and, 
for this reason, irregularities tend to average out for samples that are large compared 
to single groups of nebulae. 
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COMPLETENESS OF THE COUNTS 

Completeness of the counts was investigated by comparing 
plates of the same regions made with different apertures and ex- 
posures, the 2-hour exposures with the 1oo-inch being taken as 
standards. In this way it was possible to count all the nebulae 
actually recorded, regardless of whether or not they were identified 
in the course of the survey, on plates of 20-minute and 1-hour ex- 
posures with the 60-inch and of 1-hour exposures with the 1oo-inch. 
The numbers actually recorded, V., were then compared with the 
numbers previously identified, N,, reduced to quality E, corrected 
by the coma factor, and then designated as NV... Since the investiga- 
tions were restricted to the central regions of the plates, the two 
corrections (for quality and coma) concern the identification, and 
not the registration, of nebular images. The differences, V.—N., 
represent the numbers of nebulae recorded below the threshold of 
identification. This threshold is above the threshold of registration 
by the amount 

log V2 — log N, 
i“ 0.6 ; 





Amn 


Additional evidence is furnished by similar comparisons of 20- 
minute and 1-hour exposures with the 60-inch, the 1-hour exposures 
with the 1oo-inch being used as standards; and also by comparisons 
of the 20-minute exposures with the 1-hour exposures with the 60- 
inch as standards. 

The available material, much of which is listed in the surveys, 
includes about fifteen comparisons of 20-minute exposures, about 
ten of the 1-hour exposures with the 60-inch, and six of the 1-hour 
exposures with the too-inch, and leads to values for Am, of 
0.60+ 0.04, 0.55+0.04, and 0.57+0.06, respectively. The last value, 
representing 1-hour exposures with the 1oo-inch, is necessarily the 
least reliable. These direct tests of the completeness of the counts 
are consistent with the conclusion derived from measures of neb- 
ular magnitudes, namely, that the value of Am, is of the general 
order of half a magnitude. 

The results for the three surveys indicate a satisfactory degree of 
homogeneity and suggest no conspicuous systematic variation de- 
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pending upon limiting magnitude. For these reasons it is assumed 
that the Crossley counts and the 2-hour survey with the 1oo-inch 
behave consistently, although no direct controls are available. 
Corrections for incompleteness would probably make the limiting 
magnitudes in Table IV slightly brighter, but the amounts would be 
very small and approximately constant throughout the series. 
Therefore, they may be included among the uncertainties of the 
zero point of the magnitude scale. 

Finally, the effect of systematic variations in plate densities aris- 
ing from the different exposure times has been investigated. It 
seemed possible that an optimal density for the identification of 
faint nebular images might exist, which would affect the counts 
systematically. Fortunately, the extensive material in the 1-hour 
exposures with the Mount Wilson reflectors exhibits a range in plate 
densities sufficient to test this possibility. The results indicated that 
no appreciable systematic effects would be expected except in a few 
of the densest 2-hour exposures. On these latter plates the counts 
may be slightly smaller than normal, but the effect on log N for 
the survey would probably be negligible. 

The examination of possible sources of systematic errors strongly 
suggests that the data in Table IV exhibit a maximum range in m 
and a minimum range in log V. The conclusion has an important 
bearing on the interpretation of the data. 


RELATION BETWEEN LOG NV, m, AND Am 
In Figure 1, values of log NV, hereafter designated as “log N,” 

are plotted against the corresponding limiting magnitude m. The 
five points determine the linear relation 

log N = 0.501m — 7.371, 
probably within the uncertainties of the data. Since uniform dis- 
tribution is expressed by 

log N = 0.6(m — Am)+C, 
where Am is a correction to the observed magnitude required by 
red shifts, etc., the deviation from uniformity, expressed by Am, is 


Am = 0.165 m+ Le ; 
0.6 
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This relation represents the differential corrections over the ob- 
served range but is not suitable for extrapolation much beyond that 
range. The constant C is known to be of the order of —9; hence the 
correction is zero in the vicinity of m= 16.5 and negative for brighter 

















16 mi(limiting mag) 19 20 21 


Fic. 1.—Apparent distribution of nebulae in depth. The straight line marked 
“Uniform Distribution” is the relation log N (per square degree) = 0.6m — 9.052, 
which represents uniform distribution with the density equal to that observed in the 
vicinity of the galactic system. The curve marked “Observed Relation” is a least- 
squares solution representing the mean results of the five surveys, on the assumption 
that apparent departures from uniformity, when expressed as corrections to the ob- 
served limiting magnitudes, are linear functions of the distance. The solution is log N = 
0.6 (m — Am) — 9.052, where log Am = 0.2 (m — Am) — 4.239. The curve marked 
“Departures” exhibits the corrections, Am, as a function of the observed limiting mag- 
nitudes, m. 


limits. For the Harvard survey of bright nebulae (m=12.9) the 
correction would be about —o.6 mag., which seems to be definitely 
inconsistent. 

A more significant representation is derived as follows: the cor- 
rection to the observed magnitudes introduced by uniform space 
absorption would be a linear function of the distance. Corrections 
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due to red shifts will also approximately follow a linear function, 
because red shift itself is a linear function of distance; and the cor- 
rections, expressed in terms of magnitudes, approximate linear func- 
tions of red shifts. Finally, any systematic deviations from uniform 
distribution would be partially compensated by corrections to the 
magnitudes, which were linear functions of distance. Hence, as a 
first approximation, we may assume that Am is a linear function of 
distance.” Then 
log Am = 0.2(m — Am) +), 


and the corresponding expression for uniform distribution is 
log N = 0.6(m — 10°-2(™— 4m) +b) 4. C, 


The values of b and C are readily found by successive approxima- 
tions when the equation is written in a form suitable for least- 
squares solutions. As a first approximation, Am in the exponent 
may be ignored. As a second step, the values of Am derived in the 
first step are introduced in the exponent. The third approximation, 
with values of Am derived in the second, appears to be sufficient; 
further approximations do not change the significant figures in the 
results. The solution, expressed as the mean of the regression- 
curves, is*$ 

C = — 9.052 + 0.005, 


b = — 4.239 + 0.008. 


2 The term “assume” in the statement should be stressed because the empirical 
data are insufficient to distinguish definitely a linear relation between Am and distance, 
from a linear relation between Am and m. As usual, curves of many different forms 
would satisfactorily represent the observed points in Fig. 1, and the selection is guided 
by information from other sources. The apparent distribution must evidently be cor- 
rected for known effects on apparent luminosities. The chief, if not the only, effects 
known at the moment are those due to red shifts, whose calculated values closely ap- 
proximate the selected relation Am=Constant X Distance, over the range of the sur- 
veys. Later it will be shown that these effects alone account for apparent departures. 
If other effects are postulated, they must be arbitrarily adjusted to compensate one 
another. This evidence justifies the selection of the particular expression for departures, 
although the expression does not emerge directly from the observational data of the five 
surveys. 

13 Values of C previously derived were —9.12, by the writer, and —9.2, by Mayall. 
The present revision is influenced by the method of analysis rather than by the addi- 
tion of new data. 
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The solution adjusts C and 6 to fit a particular set of data; errors 
in one constant are compensated by adjustment of the other. Inde- 
pendent confirmation of either constant not only would be im- 
portant in itself but would inspire greater confidence in the other. 
A check on C would be afforded by counts of bright nebulae, for 
which the corrections Am are negligible. Surveys to m=14 to 15 are 
under way with this end in view, but the only significant data avail- 
able at the moment are those in the Harvard survey." 

These counts are believed to be complete over the entire sky to 
the tabulated magnitude, 12.9. Accurate photoelectric-cell measures 
by Stebbins and Whitford's indicate that the limiting magnitude 
should be corrected by —o.1 mag. and that the dispersion is not 
more than o.2 mag. The polar caps offer two samples, but one (the 
northern) is dominated by the great Virgo cluster and the group 
south-following, which Shapley has called the ‘Extension to the 
Coma-Virgo cluster.’’ The influence of the aggregation may be 
avoided by omitting the regions, R.A.=12" to 13", Dec.=o° to 
+20°, and R.A.=12"5 to 1355, Dec.=o° to —15°. The two samples 
combined then contain 592 nebulae (m= 12.8) scattered over 14,270 
square degrees. Therefore, log N = —1.38 and C= —9.06. 

The uncertainties are difficult to estimate because the data are 
rather meager. The mean latitude correction, which would reduce 
the numerica] value of log N by about 0.03, is probably compen- 
sated by clustering outside the omitted regions. The evidence for 
such clustering is fairly strong. The polar caps furnish two samples 
of the general order of those in the 2-hour survey, in which o= 
0.084. Therefore, the probable error for the mean of the two caps 
might be estimated as of the order of 0.04 or slightly less. 

Another way of presenting the data is to say that the solution 
from the five surveys indicates a limiting magnitude of 12.78 for 
the Harvard counts, which agrees with that actually observed, well 
within the uncertainties of the data. It follows that the inclusion 

14 Shapley and Ames, ‘‘A Survey of the External Galaxies Brighter than the Thir- 
teenth Magnitude,” Harvard Ann., 88, No. 2, 1932. 

S Annual Report of the Mount Wilson Observatory, 1935-1936. 
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of the Harvard survey in the general solution would not materially 
alter the values of C and 6 as derived from the five surveys alone.” 


RELATION BETWEEN Am AND RED SHIFTS 

The corrections Am represent observed departures from apparent 
uniformity arising from all sources. Nevertheless, they are so small 
that, for a preliminary discussion, space absorption and real de- 
partures from uniformity may be ignored and the corrections re- 
ferred to effect of red shifts alone. Since red shifts, dA/A, and the 
corrections, Am, are both approximately linear functions of distance 
(over the observed range), one can be expressed as a linear function 
of the other. The velocity-magnitude relation”? is 


log v = 0.2(m — Am) +0.77, 


or 
log o = 0.2(m — Am) — 4.707, 
and 
log Am = 0.2(m — Am) — 4.239. 
Hence 
dd 
Am = 2.94 x” 


The probable error of the coefficient in the last relation is rela- 
tively small; but, in view of the possibility of systematic errors, a 
rough estimate of the uncertainties is more significant. The uncer- 
tainty of the constant in the first two equations has been previously 
estimated as not more than 0.02; and, if the same reliability is as- 
signed to the constant 0, the uncertainty in the coefficient will be 
not more than 1o per cent (about three times the probable error). 

The most questionable datum is probably the zero point of the 
magnitude scale; but, fortunately, errors in this quantity do not 
seriously affect the coefficient. For instance, arbitrary revisions by 

© With log N= —1.38 and m=12.8 for the Harvard counts, a least-squares solution 
for the six groups leads to C= —9.053, b= — 4.241. 

17 Mt. W. Contr., No. 549; Ap. J., 84, 270, 1936. 
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as much as 0.3 mag. in either direction change the coefficient by less 
than 15 per cent. In the past, revisions of nebular magnitudes have 
rather consistently reduced the numerical values. A similar revision 
of the present data would increase the coefficient; but the change, 
it is believed, could not be very great. 

Other factors involved in the calculations have already been men- 
tioned. The observed range in limiting magnitudes is probably a 
maximum, and any revision would be expected to reduce the co- 
efficient. The completeness of the nebular counts appears to be com- 
parable in the various surveys. The linear velocity-distance rela- 
tion has been extrapolated well beyond the observed range; but 
three of the surveys fall within the observed range, and these alone 
furnish a value of the coefficient which is comparable with that 
derived from all five surveys as a group. 

The relations 

log N = 0.6(m — Am) — 9.052, 
log Am = 0.2(m — Am) — 4.239, 


Am = 2.94 cia 
summarize the observational results. They are formally consistent 
with the assumption that apparent departures from large-scale uni- 
formity, Am, are linear functions of distance and represent effects 
of red shifts on apparent luminosity. Therefore a comparison of the 
observed value of B in the relation Am =B dX/X with values calcu- 
lated on various possible interpretations of red shifts may be ex- 
pected to indicate which interpretation is the most probable. 
II. EFFECTS OF RED SHIFTS ON APPARENT LUMINOSITY 
The apparent luminosity of a nebula is measured by the rate at 
which its light-quanta reach the observer and by the energy in the 
quanta. Each of these quantities may be influenced by red shifts, 
and the results are known as the “number effect” and “energy 
effect,” respectively. 
Number effect—lIf red shifts are velocity shifts, the recession of 
the nebulae reduces the rates at which the quanta reach the observer, 
and hence the apparent luminosities, by the factor 1+v/c. This fac- 
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tor, the number effect, operates only in the case of actual motion. 
It is constant throughout a given spectrum because it merely re- 
duces the vertical, or intensity, scale of the energy distribution- 
curves. Since it is non-selective, the number effect increases ap- 
parent magnitudes on any system—bolometric, visual, photo- 
graphic, etc.—by the same amount, namely, 


Amy = 2.5 log (: oo 4 : 
Cc 


In expanding models of the universe obeying the relativistic law of 
gravitation, the ratio v/c is replaced by dX/X. 

Energy effect—The energy effect follows from the familiar rela- 
tion EX=constant. Red shifts, regardless of their interpretation, 
evidently reduce the energy of individual quanta and, consequently, 
the energy of the sum total of all the quanta that reach the earth’s 
atmosphere, by the factor 1+dd/X. The energy effect is selective, 
for it changes the horizontal or wave-length scale of the energy- 
distribution curves. Bolometric magnitudes are increased by the in- 
crement Amg=2.5 log (1+dX/X), but the effect must be traced 
through the atmosphere (selective absorption) and the telescope 
(selective reflection) to the photographic plate (selective sensitivity) 
in order to be evaluated in terms of photographic magnitudes. For 
convenience, the reduction factor, expressed as a magnitude incre- 
ment, may be termed K. Then, the total effect of red shift on ap- 
parent photographic magnitude is to increase the magnitude by 
the increment 


! dd : 
Am slog (r+) +4 
or by 


Am 


2.5 log (2 ft *) +K, 


depending upon whether or not red shifts are velocity shifts. 


EVALUATION OF K 
The evaluation of K depends upon the distribution of intensity 
in nebular spectra. The precise forms of the curves are not well 
known, but throughout the photographic region they resemble that 
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of the solar spectrum. The agreement suggests that nebular radia- 
tion may be discussed in terms of black-body radiation with about 
the same degree of confidence that is inspired by similar treatments 
of solar and stellar radiation. Red shifts distort an original black- 
body curve for temperature 7, into new black-body curves cor- 
responding to T=7,/(1+dd/Xd). The lower temperatures introduce 
increments to the effects of atmospheric absorption, to the heat 
index, and to the color index corresponding to the original radiation 
at 7,; and the sums of the increments are the desired values of K. 


Thus, 
K = A(Am,) + AHI + ACI. 


The increments are readily evaluated from curves representing 
Am, (reduction to no atmosphere, including reflection in the tele- 
scope), HJ, and CI, for stars, as functions of effective temperature."® 

The values of Am, vary less than 0.1 mag. over a range in black- 
body temperatures from 6000° to 4000°; hence uncertainties in T, 
will introduce only small uncertainties in the increments A(Am,). 
The temperature scale for color indices differs widely from the scales 
derived from heat indices and from water-cell absorptions; but the 
spectral type is a common denominator by means of which 7,, the 
zero point for the present purpose, may be transferred from one scale 
to another. The increments, AH/J and ACT, are then relatively in- 
sensitive to the particular scale employed. Actually the scale de- 
rived from water-cell absorptions was used for deriving AHJ; and 
another, the familiar black-body scale,'’® was used for deriving AC/ 
in order that the results for K might, with some justification, be 
regarded as minimum values. Table V lists values of K determined 
in this manner for various values of 7, and of the red shift. 

Another, and more purely computational, method of evaluating K 
may also be employed.” A black-body curve for temperature 7, 

18 The data for Am, and HI are given in Tables IV and V of “Stellar Radiation 
Measurements,” by Pettit and Nicholson, Mt. W. Contr., No. 369; Ap. J., 67, 279, 1928. 
Color indices are conveniently listed on p. 734 of Astronomy (1927) by Russell, Dugan, 
and Stewart. The values of Am,, reduction to no atmosphere, include the reflections 
in the telescope—two silver surfaces. 

19 Russell, Dugan, and Stewart, op. cit., Table XXX. 

20 The method is described by de Sitter, B.A.N., 7, 205 (No. 261), 1934, and by 
Hubble and Tolman, Mt. W. Contr., No. 527; Ap. J., 82, 302, 1935. 
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is distorted by factors representing atmospheric absorption, reflec- 
tions in the telescope, and plate sensitivity. The area under the dis- 
torted curve represents the intensity actually recorded on the plate. 
The original curve is then displaced by red shifts and later distorted 


TABLE V 
K AS A FUNCTION OF J, AND dX/X 











To 
dr/r ea a a a etic a 
| 
6500° | 6000 6000°* | 5500 
| 
: = 
0.05 0.08 0.09 0.09 0.12 
°75 Il | 13 14 18 
EO... 15 1d 20 25 
125 19 23 25 31 
15 .24 28 31 .38 
175 28 | 34 30 44 
20 33 | 40 2 $2 
225 .38 47 .48 60 
0.25. 0.44 0.53 0.55 0.68 











* De Sitter, B.A.N., 7, 205 (No. 261), par. 4, 1934. 


by the factors previously applied. The ratios of the various recorded 
intensities to the intensity for zero red shift indicate the energy 
effects, 


dx 


Amr = 2.5 log (x +  ) +e, 


from which the values of K are readily derived. 
This method was used by de Sitter?" in computing Am, in terms 
° . . ° ry ° sf Oo om ath 
of pg, pv, and vis magnitudes for T,=6000° and 5000’. The results 
were expressed in the form of power series, those for photographic 
magnitudes being 
dx 
7, 


dX 


: iiciees 


+ | 


(T. = 6000°) Am = 2.90 


(T. = 5000°) Am = 4.05 40.8 (Sy + ows 


21 B.A.N., 7, 205 (No. 261), 1934. De Sitter assumed one reflection from silver and 
one from aluminum, but the difference from the two silver reflections used in the 
surveys is negligible. 
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The two methods of determining K lead to results which are in 
substantial agreement. The greatest uncertainties are those intro- 
duced by the assumption of black-body radiation and the neglect 
of line absorption. Appropriate revisions will be suggested after 
the characteristics of integrated nebular radiation have been dis- 
cussed. 

The data in Table V could be expressed as power series, following 
de Sitter’s example; but for the present purpose, only the first terms 
are required. Out to the limits of the surveys, say to d\/A=0.22+, 
the three quantities, K, Amz, and Amy are all nearly linear func- 
tions of d\/. Approximate values of the coefficient B in the relation 


TABLE VI 
B* AS A FUNCTION OF 7, 

















B 
To 
K+Ampt K+Amp+Amyt 
GROO™ & os.s1¢, 00 sev 2.0 3.0 
7000 2.4 3:4 
ne 2.05 3.05 
> 3.0 4.0 
5500. 3.6 4.6 
5000........-.- | 4-45 5-45 





* B is the coefficient in the relation Am=B dd/2. 
tIf red shifts are not velocity shifts, Am=K+Amp,= 
K+2.5 log (1+dd/)). 
tIf red shifts are velocity shifts, Am=K-+Amp+Amy = 
K+5 log (1 +dX/a). 
Am=B dd/X are readily obtained from the values of Am for a 
given red shift; at d\/A=o.2, for instance, B=5 Am. These values 
are listed in Table VI for red shifts interpreted both as velocity 
shifts and not as velocity shifts and for various effective tempera- 
tures. 
EFFECTIVE TEMPERATURE OF NEBULAR R:DIATION 
The choice of an effective temperature for integrated nebular 
radiation is rather arbitrary. Spectra, in general, are recorded on 
small scales and are confined to the nuclear regions. Spectral types 
vary systematically with nebular types; but the range is small, the 
large majority lying between F5 and G5. Dwarf characteristics are 
conspicuous in the two spectra (M 31 and M 32) which have been 
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recorded on considerable scales, and are presumed to be normal 
features.” 

Colors of nebulae also vary systematically with the nebular types, 
and, except for the later spirals, indicate pronounced color excesses, 
only a fraction of which can be attributed to selective absorption 
within the galactic system.”? The measures, by Stebbins and Whit- 
ford, were, in general, confined to the central regions of the nebulae 
and did not include the outer portions. 

Mean spectral types (nuclear regions), mean color classes (cen- 
tral regions), and relative frequencies for different nebular types are 


TABLE VII 


SPECTRAL TYPES AND COLOR CLASSES OF NEBULAE 




















Nebular Relative Spectral | Tot Color | T 

Type Frequency* Typet si Class] | ” 
SG 2) Seas 17 G3.6 5710 26 | 4580° 
Rc RURIRE Sony 0)-2:'gnasteb tee 19 G3.4 5730 g5 4700 
Ao) SER | a A 26 G1.6 5870 24 4860 
BG) Ie. ci lenis greed 38 F8.8 6120 {7 | 5800 
RGA: ui ee ay oe G1.2 5910 g1.9g | 5140 











* The frequencies are derived from types estimated on large reflector plates of about 650 objects in the 
Harvard survey of bright nebulae. The rare Irr nebulae are omitted. 

t The spectral types are presumably dwarfs. Data by Humason, Mt. W. Contr., No. 531; Ap. J., 83, 
10, 1936. 

tT. in the fourth column indicates the effective temperature corresponding to the spectral type, 
while 7. in the last column indicates the temperature corresponding to the color class. 

§{ The color classes are on the scale of giant stars. Data by Stebbins and Whitford, Annual Rept. of the 
Mt. Wilson Obs., 1935-1930. 


listed in Table VII, together with the corresponding effective tem- 
peratures. The average nebula in the general field is about Sb, 
spectral type dG1.2, and color class g1.g (on the scale of giants). 
The mean effective temperatures, on the particular scale employed, 
are 5910 from the spectral types and 5140° from the color classes. 
These data permit a second approximation to the equivalent black- 

22 Exceptions are found in the rare nebulae whose nuclei appear stellar or highly 
condensed and give emission or early-type spectra. 

23 Investigations of the colors of globular clusters by Stebbins and Whitford (M?#. W. 
Contr., No. 547; Ap. J., 84, 132, 1936) suggest a cosecant relation between color excess 
and galactic latitude, but the selective effect would be a small fraction of the total 


obscuration. It is not impossible that the correlation may be with position within 
clouds rather than with latitude. 
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body temperature that should be used in calculating effects of red 
shifts. 

The data in Tables V and VI indicate the qualitative behavior of 
Am as a function of 7,, but purely empirical revisions may be made 
in order to adapt the results to the actual conditions of the surveys. 
The revisions involve variations in reflectivity of silver surfaces, 
departures from black-body radiation, and color excess. 

The tabulated values of Am were based upon reflectivity of freshly 
burnished silver surfaces, while the surveys were made with mirrors 
which were resilvered about every six months and burnished about 
every month. Reflectivity in the violet falls rapidly for a short time 
and then approaches a fairly steady state. Precise information on 
variations of the Mount Wilson mirrors (during the era of silver sur- 
faces) is not available, but measures of laboratory mirrors by Strong”4 
indicate the order of the deterioration. Curves for reflectivity as a 
function of X for silver surfaces 13 and 19 days after burnishing are 
closely parallel, and the curve for the longer period is believed to 
approximate the actual conditions of the surveys.” 

Departures from black-body radiation have been determined in 
detail only in the case of the sun. Measures of the spectral region in 
question—from about 6000 A to 3000 A—by Abbot, Pettit, and 
others furnish an empirical curve for the intensity distribution, 
which is the curve for the continuous radiation as depressed by line 
absorption.” In the nebulae, the line absorption is presumably com- 
parable with that in the sun, but the distribution of continuous 
radiation is altered by the observed color excess—of the order of 

24 Ap. J., 83, 40, 1936, supplemented by unpublished data very generously furnished 
to the writer. 

25 It makes little difference whether the 13- or the 19-day curve is used, because the 
significant feature, for the present purpose, is the form of the curve rather than the 
precise values of the ordinates. Some evidence in favor of the longer interval was de- 
rived from a silvered mirror on which Humason selected a region that, in his opinion, 
appeared to resemble the average state of the large mirrors on Mount Wilson. Strong 
very kindly measured this selected region and found a curiously close agreement with 
the 19-day curve. 

26 Abbot, Fowle, and Aldrich, Smithsonian Misc. Coll., 74, No. 7, 1923. Pettit, Mt. 
W. Contr., No. 445; Ap. J., 75, 185, 1932. See also a discussion of these and other 
data by Mulders, Zs. f. Ap., 11, 132, 1935. 
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0.15 mag. (see Table VII). The source of the excess”? is unknown, 
but the effect will clearly reduce the equivalent black-body radia- 
tion for the nebulae. The effect may be approximated, following a 
suggestion by Baade, by using observed intensity-curves near the 
limb of the sun instead of the curve for the integrated light. 

The graphical analysis, based on Pettit’s measures for integrated 
sunlight and Abbot’s reduction factors for a distance from the center 
of 0.95 times the radius (which approximates the desired color 
excess), together with Strong’s reflectivity-curve and a table of 
sensitivities of Eastman 4o plates furnished by the Eastman Kodak 
Laboratory, gives a relation between the energy effect, Amp, and 
dd/X, which is not precisely linear but which fairly approximates the 
relation for black-body radiation at 7,=5750°, over the range in 
d\/X represented in the surveys. The value of Bz is of the order 
of 3.25. 

Abbot’s measures of intensities lead to a slightly larger value of 
Br, while Mulders’ data give a considerably smaller value—less than 
3.0. For these reasons it is concluded that the correct value lies be- 
tween 3.0 and 3.25, and, in view of the prevalence of early-type 
giant stars in the outer region of open spirals, is probably nearer the 
lower limit than the upper. Therefore, the round number 3.0 is 
adopted, corresponding to an equivalent black-body temperature of 
6000°, with a recognized uncertainty of the order of 200°. 

Thus, for provisional discussions, the data in Tables V and VI may 
be used directly with some assurance that the errors are within the 
uncertainties of the investigation. The uncertainties refer mainly to 
the unknown origin of the color excess and to the possibility that the 

27 The natural assumption of scattering by finely divided material in the nebulae 
leads to difficulties in interpreting the appearance of the systems as seen from a dis- 
tance. Whipple’s ingenious attempt to account for the color excess as an effect of a 
mixture of stellar types in nebulae (Harvard Circ., No. 404, 1935) is not confirmed when 
the method of calculation is applied to globular clusters. Whipple concluded that the 
data for the one cluster then available were not inconsistent; but several clusters have 
been observed by Humason, using the standard nebular spectrograph, and these data, 
together with colors by Stebbins and his associates, are inconsistent. Moreover, Smith 
(Annual Report of the Mount Wilson Observatory, 1935-1936) reports that the dis- 
tribution of intensities over the wide range of 8000 A to 3100 A, in elliptical nebulae 
(where the color excess is most pronounced), show no appreciable indications of mix- 


tures of types. 
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effect of early giants in the outer regions of spirals has been under- 
estimated—that abnormal intensities exist in the ultraviolet which, 
when shifted into the photographic region, would materially reduce 
the tabulated values of K. The former source of uncertainty is prob- 
ably not important, because different assumptions lead to compa- 
rable results. The latter source has been investigated in two ways: 
First, an ultraviolet spectrum (by Smith) of a nearby, intermediate 
spiral, extending far out from the nucleus, showed no evidence of 
conspicuous excess intensities. Second, the relative colors of elliptical 
nebulae and spirals, in clusters as distant as d\/A=0.13, were found 
to be appreciably the same as those of similar nebulae in the neigh- 
borhood of the galactic system. Since the elliptical nebulae are 
known to be homogeneous and their spectra have been recorded in 
the distant clusters, it follows that the spirals also behave normally 
—that the excess intensities in the ultraviolet are not sufficient to 
affect seriously the integrated radiation. The observations do not 
eliminate the possibility of small effects at the limit of the deepest 
survey; they merely indicate that such effects, if present, are prob- 
ably not greater than a fraction of 0.1 mag. 


III. COMPARISON OF OBSERVATION AND THEORY 


It is evident that the observed result, B = 2.94, is accounted for if 
red shifts are not velocity shifts. The comparison is based on an 
effective temperature, T,, of 6000°, but the uncertainties cover the 
range down to about 7,=5750°. The interpretation is consistent 
with the data whether the nebulae are pictured as scattered through 
an old-fashioned, infinite universe or whether a homogeneous model 
obeying relativistic laws of gravitation is adopted, in which the 
expansion and spatial curvature are either negligible or zero. The 
conclusion assumes (a) that the large-scale distribution is uniform, 
(b) that nebulae maintain fairly constant luminosities over the time 
interval involved in the observations, and (c) that red shifts are an ap- 
proximately linear function of distance over the range of the surveys. 

The argument for real uniformity follows from the small values of 
the observed departures from apparent uniformity. They are no 
greater than the effects of the red shifts calculated under minimum 
conditions. These calculated corrections, when applied to the obser- 
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vations, lead to uniformity; further corrections would indicate that 
the real distribution density increased outward, steadily and sym- 
metrically, in all directions. The latter picture is unwelcome be- 
cause it implies a curiously unique position for the observer.”* For 
this reason uniform distribution is postulated and space absorption 
is assumed to be negligible. The whole of the observed departures 
is then referred to effects of red shifts alone. 

As for the constancy of nebular luminosities, the question is 
whether or not Juminosities of spirals change materially (say 10 per 
cent, or o.1 mag.) during the time required for light to travel from 
the limit of the deepest survey to the limit of the shallowest survey 
an interval of the order of 250 million years. No definite informa- 
tion is available, of course; and probabilities depend upon specula- 
tive arguments. Nevertheless, very few students will hesitate to 
adopt the assumption that systematic variation in so short an inter- 
val will be inappreciable. Nebulae are enormous systems, and it is 
reasonable to suppose that their evolution is correspondingly slow. 
Theories of stellar radiation, especially of stars with spectral types 
comparable with those of nebulae, suggest relatively little loss of 
mass or luminosity over such intervals. Moreover, the sequence of 
nebular classification may be an evolutional sequence; and it can be 
stated with some confidence, from the study of nearby nebulae, that 
the range in mean photographic luminosities, from one end of the 
sequence to the other, is not more than o.2 mag. Finally, the 
velocity-distance relation is sensibly linear over an observed range 





comparable with the interval in question. 


28 If the nebulae are assumed to be receding with constant velocities in Euclidean 
space (as a cloud expands in a vacuum), those counted in a given survey would be 
scattered through a volume V> when the light left the limit of the survey, but would 
be scattered through a larger volume, Vn=Vo(1+2/c)3, when the light reached the 
observer. Appropriate corrections to the various surveys would then be necessary in 
order to derive the momentary distribution at a given epoch. The momentary dis- 
tribution could be uniform at one, and only one, instant. The fact that the data indi- 
cate moderate uniformity near the epoch of the surveys—thus formally accounting for 
an outwardly increasing apparent density—is a remarkable coincidence that immediately 
suggests the operation of some principle of relativity. This circumstance encourages the 
recourse to models obeying the relativistic laws of gravitation, the justification for 
which is discussed in the preliminary contributions on principles and methods by Tolman 
and the writer, Mt. W. Conir., No. 527; Ap. J., 82, 302, 1935. 
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The most impressive evidence for the linearity of the velocity- 
distance relation is that offered by the great clusters at distances 
ranging out to about 250 million light-years, or approximately two- 
thirds of the distance to the limit of the deepest survey. The ex- 
pression 

dyn 
log aes o.2(m — Am) — 4.457, 
where m is the observed magnitude of the fifth brightest nebula and 
Am=3 d/X, very closely represents the data for the ten clusters 


TABLE VIII 


DEPARTURES FROM LINEARITY IN THE VELOCITY-DISTANCE RELA- 
TION WHEN RED SHIFTS ARE INTERPRETED 
AS DOPPLER EFFECTS* 





m,—Amg | log dA/A—o.2m, 
WIEN oss senkte at 0.0041 10.50 —4.483 
POU C., ciy asx sles 2% | .O127 12.92 .470 
PCCMOUR oe sien laps tier oll .O174 F063 452 
ear 0245 14.30 451 
U2" OS rere 0517 16.27 | 408 
MOO sn so in ese | 0653 16.53 439 
COESBOR.) o. 52 6555 | .0707 16.75 .445 
LS ee ery 1307 18.28 436 
0 | ee a 0.1403 18.15 —4.371 











* The data are extracted from Table I, Mi. W. Contr., No. 549; Ap. J., 84, 270, 1936, 
where full details are given. The Gemini cluster (A =150°, 8 = +20°) has been omitted be- 
cause the latitude correction, Amg, is known to be very uncertain. Magnitudes refer to the 
fifth brightest nebulae in the clusters. The corrected magnitudes, m., include the effects of 
red shifts interpreted as velocity shifts: m,=m,—Amg—4dd/d. The two Ursa Major clus 
ters have minimum weight because they are each represented by a single spectrogram of a 
single nebula. 
that have been investigated. The only residuals as large as 0.1 mag. 
involve recognized uncertainties, and there are no indications of any 
systematic trend away from the linear relation. Therefore, a linear 
extrapolation to distances 50 per cent greater than the observed 
range will probably serve as a fair approximation. Otherwise, a 
rather sharp departure from linearity must be postulated, im- 
mediately beyond the limits of the observations.” 
If the corrections Am=3dX/) led to a precisely linear relation, 
29 The most recent discussion of the relation, in the form log v=0.2 (m—Am)-+t<, is 
found in Mt. W. Contr., No. 549; Ap. J., 84, 270, 1936. 
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the larger corrections, 4d\/\, which should be used if red shifts are 
assumed a priori to be velocity shifts, would evidently introduce de- 
partures from linearity. The data are now sufficient to demonstrate 
the existence of such departures, although the numerical values are 





log dA/A-0.2mM¢ 
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Fic. 2.—Departures from linearity in the velocity-distance relation when red shifts 
are interpreted as velocity shifts. The points represent mean red shifts, d\/A, and cor- 
rected magnitudes of the fifth brightest nebula in nine clusters (Table VIII). The ob- 
served magnitudes have been reduced to the galactic poles and corrected for red shifts 
interpreted as velocity shifts, me = mo — Amg — 4dX/X. If the velocity-distance rela- 
tion were linear (as it is when red shifts are referred to some effect other than recession, 
and the appropriate corrections are 3d\/A), the differences, log d\/A — 0.2m,, would 
be constant, and the points in the diagram would determine a horizontal line. The 
systematic variation of the differences with m, indicates the departure from linearity. 
The straight line corresponds to a linear velocity-magnitude relation, log d\/A = 
0.2064m, — 4.545, adjusted to the data, and the curved line to the power series, dA/A = 
1.75 X 10-9 + 2.7 X 107'r?, where the distance, 7, is expressed in parsecs. The 
two outstanding residuals refer to the two Ursa Major clusters, which have minimum 
weight because they are each represented by a single spectrogram of a single nebula. 


still uncertain. Mean red shifts and observed magnitudes (reduced 
to the galactic poles) of the fifth brightest nebulae in nine clusters are 
listed in Table VIII. The significant correlation is between dd/X 
and the corrected magnitudes, m.=(m,—Amg) —4dX/X. 
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Departures from the linear velocity-distance relation, log d\/X 
—o.2 m,=C, are shown in Figure 2, where C is plotted against m.. 
The straight line corresponds to a linear velocity-magnitude relation, 


dd 
log oo ©. 2064m, — 4.545 . 


The curved line represents the expression of red shifts as a power 
series of the distance, 
dy 


x = 1-75 xX 10° + 2.7 X 10°", 


where r is expressed in parsecs.3° The two outstanding residuals 
refer to the two Ursa Major clusters which have minimum weight. 

The power series will be adopted in the following discussions be- 
cause it conforms with the notation in theoretical investigations of 
expanding models of the universe. Since the second-order term is 
definitely positive, the possible models are restricted to those in 
which the rate of expansion has been diminishing during the past 
several hundred million years. The numerical value of the coefficient 
is rather narrowly limited (probably between 2.4 and 3.4X 107") 
over the range of the observed clusters, and extrapolations to the 
limit of the deepest survey (about m,=18.8 in Fig. 2) should not 
introduce serious errors. Further extrapolations rapidly become 
very uncertain. In view of the importance of the second-order term 
for theoretical investigations, it may be emphasized that the em- 
pirical value is rather precisely the effect of recession which has been 
introduced as an additional assumption not necessarily demanded 
by the observations themselves. 


HOMOGENEOUS UNIVERSE OF GENERAL RELATIVITY 


A detailed comparison may now be made between observational 
results and world-pictures for homogeneous models of the universe 
obeying the relativistic laws of gravitation. The principles and 

30 Since the mean absolute magnitude of the fifth nebula in clusters is about M pg = 
— 16.45, the distances of the clusters in Table VIII are log r=o0.2 mce+4.29. For isolated 
nebulae, where M = —15.15, log r=0.2 mc+4.03. The first term in the power series is 
essentially that derived in the earlier investigation from resolved nebulae in the general 
field. 
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methods have been discussed in a preliminary contribution** by 
Tolman and the writer, where the theoretical formulae derived by 
Tolman* were adapted for application to observational data. It 
was found that the number of nebulae, NV, to the observed limiting 
magnitude, m,, would be given by the relation 


log N = 0.6(m, — Am) + F+C, 
where 
dx 


Am = 5 log (, xa \ 


II 


yee 
or 


Am = 2.5 log (: + 5) +K, 


depending upon whether or not red shifts are velocity shifts measur- 
ing the rate of expansion of the model. Further, 


F=\log | = (sin? r-xVi-#)| ; 


where x =r/R, represents the distance, 7, expressed as a fraction of 
a constant, R,, commonly identified as the radius of curvature of 
space. In expanding models, R, refers to a particular epoch. 


TABLE IX 


EFFECTS OF SPATIAL CURVATURE EXPRESSED AS 




















7] 
r/Ro F/o.6 | r/Ro F/o.6 

- ipieniatntla || = — Sn 
GiGE ue Aik 0.001 | Oe -( eeee ae 0.075 
igs tears } .002 6 .OQ2 
rere citer .005 .65 II2 
Re ee eee .O10 By 2 135 
25 .O14 Ry 162 
3 .020 8. 196 
: ncaa .028 85. 238 
4 .037 9 292 
45 | 048 || 0.95 -373 
0.5 | 0.001 0.620 





3 Mt. W. Contr., No. 527; Ap. J., 82, 302, 1935. The discussion includes a pre- 
iiminary statement of the observational findings, which are now presented in detail. 


32 Relativity, Thermodynamics and Cosmology, Clarendon Press, 1934. 
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The term F is equivalent to a magnitude increment, F’/o0.6, the 
values of which are listed in Table [X for successive values of x= 
r/R,. An inspection of the table indicates that the effects of curva- 
ture may be neglected for values of « less than about one-third and 
that they do not become important until x is greater than one half. 
Moreover, for x=c.9 and greater, the effects dominate the apparent 
distribution. These features suggest caution in resorting to curva- 
ture to account for observed residuals. No significant contributions 
can be expected unless R, is of the same general order as the pene- 
trating power of existing telescopes. 

TABLE X 


RED SHIFTS NOT VELOCITY SHIFTS 


log N —o.6 (mo— 3d/X) = Constant 








| CONSTANT 
0.6% . | 
Me dd/d | - | log N : E as 

| (Mo— 3dr d)| | | 

| | Cal. Obs. | CcC-O 
21.03.. 0.230 | 12.204 3.162 | —g.042 | —9.050} +0.008 
We crass ca 158 | 11.728 2.686 | 042 | 033 | —0.009 
oe eee 125 | IE.Ar4. | 2.342 | 072 | 077 +0.005 
Oc eer 107 | Ir.207 | 2:101 | 046 | O51 | -+0.005 
“St ty ara ge 0.086 | 10.927 1.886 | -—9.041 | —9.044 | -+0.003 








* Calculated from the velocity-magnitude relation, log v=0.2 (m—Am) +0.77, from which 5 log dA/A+ 
3dA/X = Mo — 23.535. 


Case I. Red shifts not velocity shifts —The model employed is 
Einstein’s static universe which, for the purpose, is equivalent to an 
expanding model whose rate of expansion is inappreciable. The 
data in Table X show that the model is thoroughly consistent with 
the observations (B = 3.0, as compared with 2.94), provided curva- 
ture is negligible. The comparison is based on 7,=6000°; but, as 
previously mentioned, the uncertainties permit a considerable range 
ie 

Case IT. Red shifts are velocity shifts and measure the rate of ex- 
pansion of the universe.—lf curvature is negligible, these models are 
definitely inconsistent with the observations. The values of B are 
4.0 (calculated) and 2.94 (observed). The residuals, C—O=dyd/X 
(approximately), range up to 0.25 mag. and are equivalent to the 
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number effects introduced by the assumption that red shifts measure 
actual motion. 

If the discrepancies are due to curvature,’ the effects must be 
just sufficient to balance the number effects. The curvature must be 
real and positive, and very great—-the model must be closed and 
small. The calculation of R, starts from the value, C=—g9.05, 
derived from the direct, least-squares solution of the empirical] data, 
supported by the Harvard survey of bright nebulae and consistent 
with the static model with negligible curvature. The empirical cor- 
rections, Am,, are then compared with the corrections for red shifts, 


calculated from the relations 


dd 
Am, = * ; 
dd 
as me SS Oy my SZ —~—18,02 
a ‘78 Mid ¥-- 2.7% ta". 


log r = 0.2(m — Am.) + 4.03. 


The differences, Am.—Am,, which are included in Table XI, are 
the residuals which must be compensated by curvature. 

The next step is to plot the curvature increments, F/o.6 (Table 
IX), as a function of log r/R,. The curve is shown in Figure 3. 
The range in log r for the limits of the various surveys (Table XI) 
is about 0.40. Therefore, the curve in Figure 3 may be searched for 
some range of the order of 0.40 in log r/R, which will furnish the 
desired series of corrections, F/o.6. 

No segment of the curve will fit the tabulated residuals exactly. 
Nevertheless, a segment does exist, log r/R, = —o.11 for the deepest 
survey, which furnishes increments about equal to the residuals 
minus a constant, 0.075 mag. This fit (Table XI), which is rather 
narrowly limited, may be accepted as a possible solution on the as- 
sumption that the small constant, 0.075 mag., reflects a considerable 

33 Since the formulae for the world-pictures were derived by means of co-moving 
space-like co-ordinates, with respect to which the nebulae have no relative systematic 
motion, the differences in time required for light to reach the observer from the limits 
of the various surveys is already taken into account (Mt. W. Contr., No. 527; Ap. J., 


82, 302, 1935). 
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error in the zero point of the magnitude scale (true magnitudes 
brighter than observed magnitudes). 


TABLE XI 
OBSERVATIONS CORRECTED FOR SPATIAL CURVATURE 
IN AN EXPANDING UNIVERSE 





























D= 
Me dd/d Am,* ( z V/o. | D—F/o.6 
ri m Amo ie ailiii log rt F/o.6t | D—F/o0.6 
21.03 0.231 0.924 0.676 0.248 8.051 0.180 | 0.068 
i .158 630 .468 162 7.904 77 085 
[0.8 ;.. 124 500 . 368 £32 7.810 047 085 
10... 106 .423 .314 109 7.746 034 075 
18.47 0.085 0.338 0.253 0.085 7.656 0.022 | 0.063 
Mean | eT | eee OPI SO | (er Rena eae | teeeeeee| 0.075 
| 





* Am.=4 dd/X. 

t log r (in parsecs) refers to the observed limiting magnitude of the surveys. 
_ _} Corrections for curvature are taken from Table VIII on the assumption that log r/R. = —o.11 for the 
limit of the deepest survey (m.=21.03). 
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0.2 


10/6 F 














-0.8 0.6 0.4 0.2 Log F/Ro_ 
Fic. 3.—Corrections for spatial curvature. The curve shows the corrections to nebu- 
lar counts, expressed as increments to the limiting magnitudes, required by spatial cur- 
vature whose radius is Ro, in a homogeneous universe obeying the relativistic laws 
of gravitation. The magnitude increments (Table IX) are exhibited as functions of 
log r/Ro. The section from log r/Ro = —o.11 to —o.51 has been selected as furnishing 
the best fit with the observational results when red shifts are interpreted as velocity 
shifts measuring the rate of expansion of the universe. The corresponding value of Ro | 
(at the present epoch) is 4.7 X 108 light-years. 
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The positive curvature, as previously mentioned, indicates a 
closed, finite universe. The radius, R,, would be of the order of 
1.45 10° parsecs or 4.7X10° light-years—just less than the esti- 
mated mean distance of the faintest nebulae recorded with the 
100-inch reflector under optimal conditions. If the radius of curva- 
ture were of this order, it would follow that we now observe a large 
fraction of the universe.*4 

The mean density of the model would be of the general order of 
10°” gm/cm3, Any significant revision of the empirical data would 
probably reduce R, and, consequently, increase the density. Thus 
the validity of the solution, as a representation of the actual uni- 

34 The spatial curvature and the two terms in the power series for red shifts deter- 
mine a particular model of the universe among the many homogeneous expanding 
models which obey the relativistic laws of gravitation. Tolman’s calculations of the 
possible ranges within which the various constants must lie (Relativity, Thermodynamics 


and Cosmology, pp. 472-474) can now be pushed to a much closer approximation. In 
his notation, the unit of distance is the light-year and 


dx k 


~ = BF IP + (5 


ee vaestd 


The co-moving co-ordinate, 7, is equivalent to the co-ordinate r used in the present 
discussion, namely, the distance derived from apparent magnitudes when the latter 
are corrected for red shifts. The new data give the following values of the constants: 


k = 5.37 X 107% (It. yr.)", 


j= — 2.54 X 10 ™ (it. yr.) 
I 
R 4:52 X 107°" (It. yr.)*; 
A = (4.4 to 6.9) X 10°? (It. yr.)7?. 


The cosmological constant, A, is determined between upper and lower limits by the 
assumptions: (a) that pressure cannot be less than zero, and (6) that the density of 
matter is greater than zero. The constant m is undetermined; but if it counteracts, 
rather than adds to, the second-order term in the power series, it will be negative, and 
the numerical value will be not much greater than 


2.7 X 10778 (It. yr.)~3. 


Since the model is closed and the value of A is evidently greater than that for a 
possible Einstein universe, the model is probably of the type known as a “monotonic 
universe of the first kind, type M,.” It “would be an ever-expanding type which pro- 
ceeds from some singular state at Rs=o to the final state of an empty de Sitter uni- 
verse as R>”’ (loc. cit., p. 399), in other words, the type that will always be associ- 
ated with the name of Lemaitre. 
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verse, may depend rather nicely on the question of density. In the 
neighborhood of the sun the mean density is believed to be of the 
order of 1o~*55, and is probably abnormally high because of the 
local cluster. Since the galactic system, by analogy with other 
spirals, fades outward to undefined boundaries, it is reasonable to 
suppose that the mean density of matter in internebular space is 
thousands of times less than the local density. The smoothed-out 
density of material concentrated in nebulae is believed to lie between 
10-8 and 10~%°; evidence from various sources* is fairly consistent 
within this range. There is evidently no reason for assuming a mean 
density of 10~ unless the nebulae are imbedded in a relatively 
dense, pervading medium which escapes detection. 

Photometric evidence is summarized by the statement that space 
absorption over a light-path of 3.510” cm (limit of the deepest 
survey) is less than 0.1 mag., and probably less than 0.05 mag.*° A 
density of 10~*° corresponds to an absorption of less than 1 mag. by 
about 35 (or 70) grams of matter per square centimeter cross- 
section. This condition immediately rules out many forms of matter 
(for instance, dust, meteoritic material, and highly ionized gas), but 
certain forms (for instance, large chunks and non-ionized gas similar 
to the earth’s atmosphere) would still be possible. The permissible 
forms, however, are probably not those which would be a priori pos- 
tulated for the contents of internebular space. 

CONCLUSION 

The observations may be fitted into either of two quite different 
types of universes. If red shifts are velocity shifts, the model is 
closed, small, and dense. It is rapidly expanding, but over a long 

35 The limits are set by Smith’s value, Mass of the Virgo cluster + Number of nebulae 
in cluster= 2X 10" suns per nebula (Mt. W. Contr., No. 532; Ap. J., 83, 23, 1936), and 
the writer’s value, 2109 suns, derived from spectrographic rotations (The Realm of 
the Nebulae, Yale University Press, 1936). Since the present investigation leads to a 
density of about one nebula per 1.4 X 1017 cubic parsecs, the smoothed-out densities are 


10°78 and 1073, respectively. 


which increases systematically with distance. The smaller limit is probably permissible 
because any appreciable absorption would increase the curvature, and hence the 
density, required to compensate the residuals in Table XI. 
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period the rate of expansion has been steadily diminishing. Existing 
instruments range through a large fraction of the entire volume and, 
perhaps, through a considerable fraction of past time since the 
expansion began. 

On the other hand, if red shifts are not primarily due to velocity 
shifts, the observable region loses much of its significance. The 
velocity-distance relation is linear; the distribution of nebulae is 
uniform; there is no evidence of expansion, no trace of curvature, no 
restriction of the time scale. The sample, it seems, is too small to 
indicate the particular type of universe we inhabit. 

Thus the surveys to about the practical limits of existing instru- 
ments present as alternatives a curiously small-scale universe or a 
hitherto unrecognized principle of nature. A definitive choice, based 
upon observational criteria that are well above the threshold of un- 
certainty, may not be possible until results with the 200-inch re- 
flector become available. 

Meanwhile, certain features of the two solutions may be empha- 
sized. The one outstanding objection to the assumption that red 
shifts are not velocity shifts is, of course, the fact that no other 
satisfactory interpretation has yet been formulated. Nevertheless, 
the large-scale characteristics of nature are known only from the 
kind of data whose interpretation is now in question. The assump- 
tion that red shifts measure the rate of expansion of the universe is 
a long extrapolation from the familiar small-scale Doppler effects. 
Clearly, the data should be analyzed with the aid of as few assump- 
tions as possible. The observations must be corrected for energy 
effects, regardless of the origin of red shifts. The interpretation of 
the corrected data is then direct, economical], and very simple. The 
observable region is homogeneous and is evidently an insignificant 
fraction of the universe. 

The unexpected and truly remarkable features are introduced by 
the additional assumption that red shifts measure recession. The 
velocity-distance relation deviates from linearity by the exact 
amount of the postulated recession. The distribution departs from 
uniformity by the exact amount of the recession. The departures are 
compensated by curvature which is the exact equivalent of the 
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recession. Unless the coincidences are evidence of an underlying, 
necessary relation between the various factors, they detract ma- 
terially from the plausibility of the interpretation. 

The high density required by the curvature is another disturbing 
feature. Although certain specific forms of matter are not ruled out 
by the photometric data now available, it may be emphasized that 
the observations merely set upper limits to possible densities and 
otherwise have no bearing on the actual amount of matter. Dy- 
namical arguments may eventually be developed which will be valid 
regardless of the form of material. Finally, the small scale of the 
expanding model, both in space and in time, is a novelty, and as 
such will require rather decisive evidence for its acceptance. 
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COMPUTATION OF MEAN PARALLAXES, MEAN ABSO- 
LUTE MAGNITUDES, AND MEAN DISTANCES FOR 
GROUPS OF STARS FROM COMPONENTS OF THE 
PROPER MOTIONS* 


GUSTAF STROMBERG 


ABSTRACT 

Formulae are derived for the computation of mean parallaxes, mean absolute magni- 
tudes, and mean distances from the proper motions of stars, projected on the right 
ascension and declination axes, on the antapex direction and the perpendicular axes, 
and on the axis of galactic longitude. Formulae for both the parallactic and the peculiar 
motions are derived for all the co-ordinate systems mentioned above. The general 
formulae are given for the derivation of the mean errors in mean parallaxes. 

Formulae for computing the dispersion in parallax, distance, and absolute magnitude 
are also derived. 

Mean parallaxes are in general computed from v- and 7-com- 
ponents of the proper motions, combined with values of the velocity 
of the group of stars relative to the sun and with the velocity disper- 
sion. The usual derivation of v- and r-components involves, in suc- 
cession, the determination of total proper motion and its position 
angle, the angular distance from the star to the apex, the position 
angle at the star of the great circle through the antapex, and the 
projection of the proper motion on this great circle and on the per- 
pendicular direction. Although usually facilitated by graphs or 
tables, the computation is, nevertheless, rather laborious. The writer 
suggested to the Commission for Stellar Statistics reporting at the 
Paris meeting of the International Astronomical Union in 1935 that 
tables be prepared for the derivation of mean parallaxes directly 
from the proper-motion components in right ascension and declina- 
tion.' Apart from the simplicity of the computation, this method 
has the advantage that we can give proper weights to the equations 
of condition and easily correct the mean parallaxes for the effect of 
systematic errors in proper motion and of galactic rotation. The 
proposal was adopted by the Commission, and the tables will be 
constructed at the Lund Observatory. 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 558. 


1 Trans. I.A.U., §, 230, 1935. 
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The following pages give a general derivation of the formulae 
needed in the computation of mean parallaxes, mean absolute mag- 
nitudes, mean distances, and the dispersions in these quantities. 


GENERAL FORMULAE 


The equatorial rectilinear co-ordinates relative to the sun for a 
star of right ascension a, declination 6, and distance KR are 


& = Rcosacosé, 


(1) 


n = Rsinacos6é, 


| 
¢=Rsin6. 


Differentiating with regard to the time and introducing the radial 
velocity, V, expressed in kilometers per second, the parallax 7, and 
the components yu, and pw, of the proper motion, expressed in seconds 
of arc per year, we obtain 


x = V cosacosé6 — — (yu, sina +p, cos a sin 6) , 


a i 


y= Vsinacos6+~— (yu, cos a — pw, Sin a sin 4) , 


(2) 
k = 4.738 km year/sec , 
da 
= cos 6 
- dt : 
_ d 
‘a dt © 
From (2) we find 
Vv =xcosacosé6+ysinacosé+zsin6, | 
k ; | 
Mt _ _xysina+ycosa, 
T ‘ (3) 
ku: | 


= —xcosasinéd — ysinasinéd+zcos6. 
T 


1 
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The observed quantities are V, u;, and w.. From the first of equa- 
tions (3) we can determine mean values of x, y, and z, that is, the 
components of the group motion; and also the residuals in the radial 
velocities, that is, the peculiar radial velocities V’. 

The sun’s motion relative to the stars was formerly thought to be 
a fundamental constant and was early introduced in statistical as- 
tronomy; and we shall use this velocity instead of the opposite vec- 
tor, the group motion, realizing, however, that the sun’s motion de- 
pends upon the group of objects to which it is referred. Using a 
standard solar velocity V., with the components 2, Yo, and Zo, and 


putting 


x= —pPX + 1; y= —py + yy; = —p%+n, (4) 


we obtain from equations (3) 


V + pAV =7,=V’', 





Mr Na 
_ pa a » | 
T k i (5) 
M2 UT) 
= OD + —, 
T k 
where 
AV = x, cos a cos 6 + y sin a cos 6 + 2 sind = V, cos dX, | 
kA =x sina — y cosa, > (6) 


| 


kB = x, cos asin 6+ y sin a sind — 2% cosé. 


The quantity p is determined by a least-squares solution of the 
equations of condition pAV + V = o. It is equal to unity if the 
group motion projected on an axis toward the standard antapex is 
numerically equal to the standard solar velocity. The quantities 7 in 
equations (4) and (5) are the velocity components along the axes 
indicated by the subscripts, which remain after the effect of the 
standard solar velocity, multiplied by the factor p, has been sub- 
tracted. \ is the angular distance of the star from the standard apex. 

In the tables to be published, AV, A, and B will be given as func- 
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55 
tions of a and 6 and referred to the standard solar apex, a, = 270°, 
56, = +30°, with V, = 20 km/sec. In other words, 
% = 0.00, Yo = —17.32, 2 = +10.00 km/sec. (7) 
Equations (5) can be written 
Ta ) 
pAr = wr — k aes | 
(8) 
15 
pBr = wh — : = U2 — 6. 


If the stars are well distributed over the sky, the condition 2V” = 


minimum, used in deriving p, implies that nz and Ynj also have 
minimum values. This justifies the use of the conditions Ye? and 
Ye} equal to minima, even if the true apex of the solar motion for 
the group should differ appreciably from the standard apex. 

If the number of stars in the group is small, there is danger that 
an isolated large peculiar radial velocity may greatly influence the 
value of p. In some cases we know the value of p for the class of 
stars studied. In other cases we may determine / from the equation 


pr=|AV| = —=V|, (9) 


where V | is equal to the radial velocity if AV is positive, but has its 

sign changed when AV is negative. The condition used is then equiv- 

alent to AV’ = o. Double vertical bars here and in what follows 

mean that the inclosed quantities are regarded as positive. 
Projecting the proper motions on the great circle through the 

antapex and on another great circle perpendicular to the former, we 

obtain expressions for the v- and the 7-components and for X, the 

angular distance to the apex: 

oa Ap, + Bu, T aes 

+ a +B ae (pV. sind+ m), 


Ap, aa Bu, Tr 





"Vat BOR’ (x0) 
. .  RVA?7+ 3B AV\2 
aa 7m ” tal Goa ° 
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These formulae give uncertain values of v and + when A? + B? is 
small, that is, for stars near the apex or the antapex. But in this 
case uw, and w, can be regarded as peculiar motions equivalent to 7 
and can be used directly for computing mean parallaxes, as will be 
explained later. 
MEAN PARALLAXES 

The equations of condition for determining mean parallax by 

the use of the tables here referred to are 


T; 
pAr — m= m—*) Cri 
| 
f (11) 
rl, | 
pBr - ws = m — es 


These equations, apart from a change in notation, are identical with 
equations (8). We have here introduced 7, and 7., the true errors 
in the measurements of uw, and w., respectively. 7, and T, are the 
“peculiar,” linear tangential-velocity components remaining after 
the effects of the standard solar motion have been subtracted from 
M, and yw., respectively. Equations (11) can be used when only one 
of the proper-motion components is known. If uw, and pw, are de- 
termined with different accuracy, proper weights can be introduced. 
To simplify the following equations, we will write them as if we used 
only the proper motion in right ascension, although, of course, simi- 
lar equations hold for the proper motion in declination, and, in gen- 
eral, the results should be combined. 

In the general case we regard the parallax as differing from one 
star to another in the group studied, and wish to determine the mean 
parallax, 7. 

1. Parallactic motions. —We assume Ye? = min. We have then 


Le? = p?TA?n? + Lu? — 2pLAmrp: = min. 


If + is nearly constant or if there is no correlation between A and 


7, we can write 


p?r?dA? + Dwi — 2pLAmrp: = min. 
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There is always a correlation between A and u,, and also between 
m7 and yu,, and hence between Ay, and 7. Denoting the coefficient 
of correlation between Ay, and z by r, and the number of stars by n, 


we have 
LAmp, = 2(Ap, — Aps)(e — F) + FLAp, = FLAp, + nrowo,, 


= L(4 = 1)? ee D(Apy = Au;)? 
— -—_<— * —_— ef * 


Further, 


Hence 
— . > . > —~A ° 
Lei = P(e + of)DTA? + Lyi — 2pT TA. — 2pnro,o, = min. 


By varying 7, we find 


mwprA? = TAp.. (12) 


If both uw, and yw, are given for all the stars and if they have the same 
mean errors, we obtain from equations (10) and (12), 


_ 2(An + Buz) 


ey kv sin X 
= alls 


ps(A? +B’) pV.2 sin? d* (13) 


This is a convenient form of the formula generally used, which is 
ordinarily derived on the assumptions that 7 is practically constant, 
that V, is the actual group motion, and that the standard apex is 
identical with the true apex for the group motion. 

2. Parallactic motions—-We assume Ye, = 0. Since we can re- 
arrange the signs of e, in any unbiased way we please, we get from 
equation (11), provided that there is no correlation between z and A, 

tpr| A = =u,|, (14) 
where the bar after yu, indicates that the sign of yu, is reversed if A 
is negative. If the v-components are derived for each star separately, 
the same condition gives 

kdv Zu 


7 = —— (15) 


pV. sin » id prov A? + B? 
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| 


3. Peculiar motions—We assume |7,| = | V’| = 0, the average 
peculiar radial velocity. From (11) we obtain 


lex} = m+ = fx — PA]. 
To solve this equation for 7 we must know an approximate value 


of this quantity, which we will denote by 7’. We can now write 


the approximate equality 


_3 zA\ 2 2 
M1 | + (=) = |u: — pA 

Hence, 
a=" lina par (lal (16) 
T= gules pAt | - I 


As a second approximation, we replace 7’ by 7 and recompute 7. 
For stars near the apex or the antapex, A and B are small and errors in 
r’ have little effect, and we might even neglect the term pA7’ in the 
first approximation. For v- and r-components, respectively, equa- 
tion (16) becomes: 





p. | F Ks . , "i 
ia Ulm Sk) ~ ee, 
6 k | 
(17) 
ae 
oa — | 


In the last equation we recognize the usual form for deriving 7 
from r-components, obtained by neglecting the effect of errors of 
observation. 





° : r Tr > 7, ° 
4. Peculiar motions—We assume 7,? = V” and obtain 
a 7 
a ee 
eat 2 eG pAn). 


This equation can be solved by approximations. We find 





a a. 
wae [Gar — paw’ — ni] =F +7. (18) 
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The last equation gives us, after combination with 7 from other 
determinations, a value of o,, the dispersion in 7. 
The equation for the 7-components equivalent to (18) is 


Pan one tat. (19) 


V7 





If we express all linear motions in astronomical units per year, 
instead of in kilometers per econd, we can omit the factor & in all 
the preceding formulae. 

The advantages of the foregoing formulae lie in the many different 
values for 7 which we can obtain by different assumptions about the 
errors. Further, the peculiar motions are here utilized in a much 
more general way than is possible when they are measured by r-com- 
ponents alone. 

MEAN ABSOLUTE MAGNITUDES 

In many cases the mean parallax of a group of stars is in itself of 
little physical significance. Often we want to determine the mean 
absolute magnitude of a class of stars and obtain an estimate of the 
spread of intrinsic brightness. To reduce the spread in parallax due 
to inequality of the apparent magnitudes of the stars in the group, 
we must first reduce the proper motion and parallaxes to a common 
apparent magnitude. 

The absolute magnitude M, as defined by the International As- 
tronomical Union, the apparent magnitude m, and the parallax are 
related by the equation 

M=m+5+ 5logz. (20) 


We introduce a parallax 7, reduced to zero apparent magnitude and 
defined by the equation 

T™ = mw 10°?", (21) 
We have then 


M=5+5log™, ) 
or ? (22) 
Tm =0.1°10°?%, | 


From equations (11) and (21) we find 


Tol > 
pAm — m= 1 — = €, 10°" = ¢i, ( 


k 


nN 
Ww 
— 
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where yu; and ¢; are the proper-motion component and its residual, 
respectively, both multiplied by the factor 10°”, that is, reduced 
to zero apparent magnitude. 

We can now determine 7, in the same way as we determined 7 by 
use of the equations (12) to (17). From equation (18) or (19) we 
can also obtain a value of z?, which is needed in order to find the 
dispersion in absolute magnitude, which we will denote by g. We 
will assume that the absolute magnitudes are distributed according 
to a normal error-curve about a mean value M. We have then 


M=M+4A, 
10° 7 a 10° 2M * 10° 2A ' 


0.4M oa» 0.4m a 0.4 
to’ 4M = yor 4M . 70-44 | 


l "100 Po: 0.02q? (23) 
10°:24 _ ; e 2q? 10° 24d A —_ eMod? — C ‘ 
V 21g J —« 
I “409 _ 4 0.08¢q? | 
10°) = — ( e 10° 4dA = eMed? = (4, 
Vv 277. — oO ) 
Thus 
7 =0.1- 10° =o0.1- 10°C, 
m= 0.01- 10° = 0.01 - 10° C3, 
The final formulae now become 
= We 
Oak, 
To 
(24) 


g = (M —- M)? = 50 Mod log C = 21.71 logC, 





M=5+5logm — slogC. 


All the logarithms used here are to the base ro. 

Some of the formulae given here were published in Mount Wilson 
Communications, No. 103; Proceedings of the National Academy of 
Science, 14, 834, 1928, where references to and comparisons with 
‘arlier formulae by Gyllenberg are given. 
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MEAN DISTANCES 
The distance R to a star expressed in parsecs is the inverse value 
of the parallax. Since the parallaxes of the stars in a group in gen- 
eral are different, the mean value of R is always greater than the in- 
verse of the mean parallax. If the dispersion in 7 is o,, we obtain 


by expansion in terms of ¢,/7, 


i) 
On 
~~ 


Rr=1+%, ( 
7 
where o, can be found from equations (18) or (19). 

If the stars in the group have the same apparent magnitude, the 
spread in distance is due solely to a spread in absolute magnitude. 
We have then 

C 


R _—w 
a(m, M) 


R(m, M) - 107974 = 


7 = m(m, M) - 10°74 = x(m, M)C 

Here R(m, M) and x(m, M) denote distances and parallaxes, re- 
spectively, computed from m and M. We thus obtain, using one of 
the equations (23), 


0.04q? 


Rx = C? = eMot | (26) 


This formula was first given by Oort.? 
If there is a small variation in m, we can write 


=z —- _ | =r ] Der 
R = 10+ 107°? 19°2™C = — 10°C = R(m, M)C, (27) 
To 


where C is given by equation (23). 
In the general case we can only use equ. (25). 
MEAN ERRORS OF PARALLAXES 


In the preceding study we have regarded z and M as quantities 
differing from star to star. The means of m7 and of M are derived 
from a limited set of data, however, and hence are affected by errors 


2 B.A.N., 4, 81, 1927. 
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due to the use of relatively small samples. Such errors would be 
present even if 7 or M were constant. These statistical errors can 
be computed in the same way as when z and M were constant for 
the group. 

Equations (12) and (13) are based on a least-squares condition. 
Denoting by the symbol E(7) the mean error in 7 and by u the num- 
ber of stars, we therefore have 


2( pr) = ——>et_ = 2(A® = my)? 
seth ta ae (n—1)DA2 (n — 1)DA?’ (28) 
E(x) = -_ es (7 +4 )E*(p) | 


Some estimates of the mean error in ~, which enters in equations 
(28), can in general be made. 

In the case of equation (14) we might use the same mean error 
as given by the least-squares condition. We may also use the equa- 
tion 

>> = 
E(pr) = 1:752le oN (29) 
Vn(n — 1) Z|A| 
which is analogous to Peters’ formula and has the advantage that 
no squared values occur. 
Equation (15) gives the similar equation 





V 
/ .2sk>|v — —°F sin AX 
Yio = ee nance neice (30) 
: Vn(n — 1) 3V A? + B VoVn(n — 1) sind ”% 


When we are dealing with means taken regardless of signs, as in 
equations (16) and (17), we can use the equation? 


42617 
B(|#|) = 2-751) G2) 


for determining the mean error in an arithmetic mean. 
Similar equations hold for | V’| = 6, and for | u, — pA7|. 


3 Cf. Russell, on the accuracy of mean parallaxes, Mt. W. Conir., No. 215; Ap. J., 


54, 140, 1921. 
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The mean error in 7 determined from equation (16) can then be 


found from the equation 


Ex(n)_ E60), By, 


a p= (2 a _2 E( Mr — pAr ). (32) 
7 O3r 
From equation (17) we get the analogous equation 
F( E(0 ae 
7) _ ai E+(|71) bass 


In those formulae which depend on a comparison between peculiar 
angular and linear motions, some improvements can be made by 
correcting the spread in these quantities for the effect of an ellip- 
soidal velocity distribution. If the stars are well distributed over the 
whole sky, no correction is needed. In other cases a knowledge of 
the ratio of the three axes is needed to apply such a correction. 


GALACTIC STARS 
If the stars are distributed along the galactic circle,only the proper- 
motion components along this circle should be used, since in general 
we have insufficient knowledge of the linear motions perpendicular 
to the galactic plane. A and B, as well as uw, and w., should then be 
projected on the galactic circle with the aid of the galactic parallactic 
angle y, which is given in Ohlson’s tables.4 We have 


Mi =u.cosg+yu.sing =ywsin(g+y), 
A,=Acosg+Bsing, (34) 


Aipr — wi =e. 


In the first of these equations uw and y are the total proper motion 
and its position angle, respectively. 

A study of the errors in the mean parallaxes shows immediately 
that the relative accuracy of the results from parallactic and peculiar 
motion depends upon the size of the linear group motion relative to 
the linear peculiar velocities. The quantitative formulation of this 
fact was first given by Russell. Except for extreme cases, it is ad- 


4 Annals Obs. Lund, No. 3, 1932. 
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visable to use all the methods here given for determining mean 
parallaxes and then weight them properly according to their mean 
errors when forming the final means. 


CORRECTIONS FOR GALACTIC ROTATION 

The formulae for computing mean parallaxes are based on the 
assumption that the velocity distribution is the same everywhere 
within the region in which the stars studied are situated. If the stars 
are very distant, this assumption is known to be invalid. The effects 
of differential galactic rotation can be allowed for by applying correc- 
tions to both the proper motions and the radial velocities. In this 
case we need consider only stars near the galactic plane and use 
only the proper-motion component parallel to the plane. The correc- 
tions to be applied to the proper motion in galactic longitude and to 
the radial velocities are then 


/ 


Au = — cos 2(1 — J,) cos 6 + a cos | , 


AV = — RA’ sin 2(/ — 1.) cos? b. 


In these formulae / is the galactic longitude of the star, /, that of 
the center of the galaxy; 6 is the galactic latitude, R the distance in 
parsecs; and A’ and B’ are Oort’s constants for galactic rotation. 
There are many determinations of these constants, and I give here 
only the latest determinations of Plaskett and Pearce from the mo- 
tions of O and B stars and the calcium clouds: 


A’ = +0.0155 + 0.0009 km/sec parsec , 
B’' = —o0.0120 + 0.0027 km/sec parsec , 
lL, = 324-4 4 2°54. 
Mean values of R can be used, the means being computed with the 
aid of equations (25), (26), or (27). 
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A SPECTROGRAPHIC STUDY OF g PERSEI 
HERBERT F. SCHIEFER 


ABSTRACT 


The spectral variations of g Persei were measured on 32 two-prism spectrograms. 
The radial velocity of the central absorption of hydrogen shows cyclic variations in a 
period of 126.6, days with a secondary variation on the descending portion of the veloc- 
ity-curve which increases progressively with decreasing wave-length. Microphotomet- 
ric studies of the intensities of the hydrogen absorption and of the emission components 
of hydrogen show similar variations. 

At maximum positive velocity of the hydrogen absorption the red edges of the hydro- 
gen emission show a similar displacement to the red, while the violet edges show a violet 
displacement, resulting in a maximum width for the total emission and in a maximum 
intensity of the violet emission component relative to the red emission component. At 
maximum negative velocity of the hydrogen absorption the violet edges of the hydrogen 
emission show a similar violet displacement, while the red edges show a red displace- 
ment, resulting in a second maximum width for the total emission and in a maximum 
intensity of the red emission component relative to the violet emission component. The 
widths of the hydrogen emission and the components vary therefore in one-half of the 
star’s velocity period, while the relative intensity of the emission components and of the 
central absorption of hydrogen relative to the continuous spectrum show periodic and 
secondary variations similar to the star’s velocity of hydrogen absorption. 

The truly periodic spectral variation of ¢ Persei is interpreted satisfactorily by the 
hypothesis of orbital motion and rapid axial rotation of two stars, each having a very 
extensive atmosphere. The hydrogen emission is assumed to originate in the projecting 
atmospheres of the stars at right angles to the line of sight. The hydrogen absorption 
is assumed to originate mainly in the outer atmosphere between the observer and the 
fairly well-defined bright photosphere of the primary star. The pronounced secondary 
spectral variations, which vary progressively with wave-length, appear to be explained 
most readily by an apparent atmospheric expansion which is produced in the interven- 
ing extensive, rarefied, and transparent atmosphere of the companion star due to excita- 
tion by radiation from the hot photosphere of the primary star at the time it is partially 
eclipsed by the extensive atmosphere of the companion star. 


INTRODUCTION 

The spectrum of ¢ Persei, class Bope, is of unusual interest be- 
cause of the great variations of the hydrogen lines. Each line con- 
sists of a very broad absorption upon which is superimposed a nar- 
rower emission. The emission is divided into two components by a 
sharp central absorption. The emission is very strong for HB but 
vanishes rapidly for lines below H6. 

The spectral variations of this star have been extensively investi- 

t A dissertation submitted in partial fulfilment of the requirements for the degree 
of Doctor of Philosophy in the University of Michigan (1928). The general results were 
presented at the annual meeting of the Iowa Academy of Science, April, 1929. In pre- 
paring this paper for publication several additions have been made, as well as many ex- 


cisions. 
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gated by Campbell and Reese,? Cannon,’ Ludendorff,* Jordan,5 Fré- 
dette,° Lockyer,’ and Dustheimer.* The radial velocity of the cen- 
tral absorption of hydrogen has been found to vary in a period of 
126.6 days. The relative intensities of the hydrogen emission com- 
ponents vary in a similar period. The results of Lockyer in 1925 in- 
dicate a cyclic variation of intensity in a period of 21 days in addition 
to the 126-day period. Subsequent data published by him in 1926, 
however, do not confirm this short period. The widths of the hydro- 
gen emission lines show a cyclic variation in one-half of the star’s 
velocity period. Photoelectric measurements by Guthnick? during 
1916-17 indicate a cyclic variation of about o.11 mag. in a period of 
18.1 days. This short period corresponds to that determined from his 
visual observations in 1911-12. The photoelectric measurements by 
Giissow’’ during 1925~—28 indicate variations of about o.1 mag. which 
do not satisfy the previously announced periods. 

With the attachment of the two-prism spectrograph to the Ann 
Arbor 373-inch reflector in the summer of 1927, it seemed worth 
while to obtain a series of two-prism spectrograms of ¢ Persei for a 
more accurate study of its spectral variations. Thirty-two spectro- 
grams were obtained over an interval slightly more than the star’s 
velocity period. On each two-prism plate ten standard intensity 
spots were photographed by means of a tube-sensitometer. 

The results which will be presented in this paper include a study 
of the intensity variations of the hydrogen lines with the Moll micro- 
photometer, of the widths of the hydrogen emission lines and com- 
ponents, and of the measures for radial velocity of the following fea- 
tures: the central absorption of HS, Hy, Hé, He, and H¢; the cen- 
tral absorption edges of HB, Hy, and 6; the emission edges of HB, 
Hy, and H5; K of calcium; helium at AX 3964, 4026, and 4471; and 
the central absorption of ionized iron at AA 4233, 4352, 4584, 4630, 


2 Lick Obs. Bull., 1, No. 20, 158, 1902. 

3 J.R.A.S., Can., 4, 195, 1910. 

4A.N., 186, 17, 191T. 6 J.R.A.S., Can., 19, 185, 1925. 

5 Pub. Allegheny Obs., 3, 31, 1913. 7M.N., 85, 580, 1925; 86, 474, 1926. 
§ Dissertation, University of Michigan, 1927, unpublished. 

9 Veréff. Kgl. St. Babelsberg, 2, 89, 1918. 

10 4 .N., 237, 321, 1920. 
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and 4924. Some results of the intensity variations of the hydrogen 
emission, obtained from 1og unstandardized single-prism spectro- 
grams, are presented for comparison. 


A SPECTROGRAPHIC STUDY OF ¢ PERSEI 


Relative intensities—The variations with phase in the intensity 
of the hydrogen lines of g Persei are shown in Figure 1. The intensity 
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Fic. 1.—Microphotometric tracings showing the periodic variation in the intensity 
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of the hydrogen lines of ¢ Persei. 


variations were determined quantitatively with the aid of the stand- 
ard spots. The results are expressed in an arbitrary unit which is 
equal to the logarithm to the base 1o of 1.5, that is, the logarithm of 
the ratio of the aperture areas of the tube-sensitometer. 

The relative intensities based upon the individual two-prism 
plates are given in Figure 2. The period used in computing the 
phases and the epoch to which they are referred are those given by 
Dustheimer, namely, a period of 126.626 days and zero phase at 
November 18.000, 1925, G.M.T. 

The curves for the relative intensities of the emission components 
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of HB, Hy, and Hé, at the bottom of Figure 2, show that at zero 
phase the violet emission component is slightly more intense than 
the red emission component. This difference in intensity increases 
with phase and becomes a maximum at approximately 28 days. On 
the descending portion of the curves, in the region where the in- 
tensity of the red emission component is approximately equal to that 
of the violet emission component, the curves show a secondary varia- 
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Fic. 2.—Relative intensity of the hydrogen emission components and of the central 
absorption of hydrogen relative to the continuous spectrum of ¢ Persei. 


tion. Beyond phase 60 days the intensity of the red emission com- 
ponent is greater than that of the violet emission component. The 
difference is a maximum between phases 70 and 80 days and then 
decreases again to complete the cycle in 126.6 days. 

The range of variation appears to increase with the wave-length. 
The secondary variation which is pronounced in Hf is only indi- 
cated for Hy and is practically absent in Hé. The absence of a dis- 
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tinct maximum at phase 28 days for Hé and the shape of the second- 
ary variation for Hy suggest that the secondary variation occurs 
earlier and may be of greater magnitude in lines of shorter wave- 










































































length. 
TABLE I 
RADIAL VELOCITIES OF ¢ PERSEI IN KM, ‘SEC 
—— RED VIOLET 
——. . EMISSION EMISSION 
DATE PHASE — . EDGE EDGE K oF} MEAN| MEAN 
(G.M.T.) | (Days) Ca HE FE 
| 
AB | Hy us| He| Ht | HB| Hy| Hi) HB | A Hs 
| | | Z 
93 Il oclsael 
July 29.780 |112.276 |—17 |—16|—21 - hci 174| 152] 153|/—219 — 209] —182 —32) —21 —2I 
Aug. 1.803 |115.299 —18| 57) —18]— 9| —53 174| 173) 144] 245| 189) 218) + 5 5 | —12 
10.798 |124.294 |—12| —11|— Q|— 5/—31| 182] 143] 155] 218 176| 172] — 4| —28 | —16 
11.802 |125.298 |— 8 —10| — JI— 0) —24 188) 164] 154] 208 101) 196] + 4] + 3 + 7 
14.847 1.717 — 7\— 2/+1 Li Ol...) 202] ZISi- >»: EBT) «FOSl6v.cas +3, —-3]-— 
21.707 8.667 |+ 4/+ 6} +11} +16} _ 197} 180} 167; 206; 188} 149} + 5} +19 | + 8 
: 25.816 | 12.686 |+ 7 +13|+15/+22) aa 6 181} 179] 172} 199] 188] 160) +14] +22] + 2 
Sept. 1.728] 19.598 |+33|+37}+35|+38) +49 203| 1. 2I§| 20f|...-.| +34] +60 | —12 
2.753 | 20.623 |+31|/+20]+29/+44|+24] 212] 165] 149] 200] 192} 132! +43] +54 | —13 
5 770 23.640 |+26/-+26)+30/+35|+ 6 178| 196] 191] 203| 1381} 173) — 6) +30 + 8 
15.782 33.652 sal +26) +21/+29/— 5 206 IQ5| IQI| 212} IoQr BOS. <ai0dh 4S +13 
18.775 | 36.645 +27) +30] +26) J hes 196) 182| 182] 206] 202] 161 ore’ —I19 
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* The values for HB, Hy, and Hé are weighted means of measures on the central absorption edges and 
ae CC Cee 

The variations of the relative intensities of the hydrogen emission 
components shown in Figure 2 for the single-prism spectrograms ap- 
pear to be in substantial agreement with those of the two-prism 
plates. 

The variations in intensity of the central absorption of Hy and H6 
relative to the continuous spectrum at these lines are shown by the 
curves at the top of Figure 2. The variation is similar for the two 
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lines. The absorption is a maximum relative to the continuous spec- 
trum at a phase of approximately 18 days. This precedes the phase 
of maximum positive velocity, shown later in this paper, by about 
10 days and confirms the results of Lockyer. The absorption is a 
minimum at a phase of approximately go days. The secondary vari- 
ation is marked, with a minimum at phase 45 days and a maximum 
at phase 56 days. 

Radial velocities and widths of hydrogen emission—The radial 
velocities are given in Table I. The velocity-curves for the hydrogen 
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Fic. 3.—Radial velocity of spectral lines of ¢ Persei 


absorption lines, based upon weighted mean places, are given in 
Figure 3. The curves have, in general, the same form but differ in 
some respects. A striking difference is the secondary variation, 
which occurs on the descending branch after maximum positive 
velocity. It becomes progressively more pronounced with decreasing 
wave-length. It appears earlier, increases in magnitude, and per- 
sists longer for lines of shorter wave-length, which is a strong con- 
firmation of the results noted for the secondary variation of the rela- 
tive intensities of the hydrogen emission components. These pro- 
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gressive variations are shown in Figure 4 by curves representing the 
isolated secondary velocity. They suggest that this variation is 
caused by a disturbance which affects the lower atmospheric levels 
of the star more than the higher levels, assuming that the lines of 
shorter wave-length originate at lower effective levels. 

In isolating the secondary variation from the main velocity-curve, 
smooth curves were drawn in such a manner that the secondary vari- 
ations were entirely below these curves. This procedure is suggested 
by the natural trend of the curves and is justified since the smooth 
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Fic. 4.—Isolated secondary velocity of the hydrogen absorption lines of ¢ Persei 





curves thus formed are essentially the same for H6, Hy, H6, and 
also for He, provided a constant correction of — 8 kilometers is 
first applied to this curve. It should be added that the measures of 
He are rendered uncertain by the proximity of the H line of Ca 1. 
The curve for H¢ is doubtless influenced by blending with the line 
d 3888.65 of Hei. Since both the hydrogen and the helium ab- 
sorptions vary in intensity, the position of the resulting blend is 
rendered uncertain. The main variation is undoubtedly produced by 
the same cause. Similar velocity-curves should therefore be expected 
for the different hydrogen lines. Since such curves are obtained by 
isolating the secondary velocity variation in this manner, the method 
appears to be a rational one. 
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The secondary velocity variation, interpreted as radial motion in 
the line of sight, indicates a mass of hydrogen gas which begins to 
have an apparent motion of approach toward the observer immedi- 
ately after maximum positive velocity (recession) of the main varia- 
tion. This apparent motion of approach reaches a maximum value 
and then decreases again to zero at the phase where the main positive 


velocity has decreased approximately to zero. It should be noted 
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Fic. 5.—Radial velocity of hydrogen emission and absorption edges of ¢ Persei 
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that this interpretation of the secondary variation does not require 
a subsequent recession of this mass of gas relative to the observer. 
The velocity-curves for K of calcium, for the mean of the helium 
lines, and for the mean of the ionized iron lines are given at the bottom 
of Figure 3. The K line exhibits the maximum positive velocity and 
the secondary variation. For the remainder of the period the veloc- 
ity is practically constant. The line appears strong between phases o 
and 16 days. At the maximum and during the secondary variation 
it appears weak. This and the fact that it shows little variation in 
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velocity during the greater portion of the star’s period, suggest that 
its origin is chiefly interstellar. 

The mean velocity-curve for the helium lines is similar to the 
velocity-curves for the hydrogen absorption lines. The helium lines 
differ from the hydrogen lines with respect to intensity, being strong 
between phases 50 and 86 days and weak during the remainder of 
the period. A secondary velocity variation is indicated near phase 
50 days. 

TABLE II 
MEAN WIDTHS OF HYDROGEN EMISSION LINES AND 


COMPONENTS OF ¢ PERSEI IN ANGSTROMS 
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The mean curve for the ionized iron lines shows only a negative 
velocity with respect to the center of mass velocity. The velocity is 
zero at phase 28 days when the velocity of the hydrogen lines is a 
maximum. The minimum is similar to that for the hydrogen lines, 
occurring at phase 96 days with a minimum velocity of — 24 km/sec. 

The displacements of the emission and absorption edges of the 
hydrogen lines expressed as radial velocity are shown by the curves 
in Figure 5. The velocity-curves for the central absorption of hydro- 
gen are also plotted for comparison. The widths of the total emis- 
sions, of the emission components, and of the central absorptions 
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are related to the distances between the respective velocity-curves. 
The widths in angstroms for the emission lines and components are 
given in Table II, and the changes are shown by the curves in Fig- 
ure 0. 

The velocity-curves for the red emission edges of H8, Hy, and Hé 
are similar in form during the last half of the period but differ in 
the first half. The secondary variation is pronounced in H9, is only 
indicated in Hy, and is absent in H6, thus showing the same pro- 
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Fic. 6.—Widths offthe hydrogen emission lines and components of ¢ Persei 


gressive change found in the relative intensity of the emission com- 
ponents and in the total emission width. The maximum displace- 
ment to the red occurs at a later phase than for the absorption. This 
lag in phase increases for lines of shorter wave-length, a fact also to 
be noted in the widths of the total emission. The maximum dis- 
placement toward the violet, between 70 and 80 days, corresponds 
to the phase at which the total width of the emission lines is a mini- 
mum, and at which the intensity of the red emission component is a 
maximum relative to the intensity of the violet emission component. 
Between go and 100 days the red emission edges show a displace- 
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ment toward the red, while the hydrogen absorption lines show a 
displacement toward the violet. At this phase the widths of the red 
emission components are a maximum. 

The velocity-curves for the violet emission edges of HB and Hy 
are similar in form, and during the last half of the period show varia- 
tions similar to the hydrogen absorption, namely, a displacement 
toward the violet. During the first half of the period, however, when 
the hydrogen absorption lines show a displacement toward the red, 
the violet emission edges of 78 and Hy show a displacement toward 
the violet. They show a similar secondary variation, but at a later 
phase than the hydrogen absorption. 

The velocity-curve for the violet emission edge of //6 is different. 
It resembles the curves for hydrogen absorption at phases of maxi- 
mum and minimum velocity. The secondary variation between 
phases 40 and 70 days is exceedingly great. The secondary minimum 
at 44 days corresponds to the phase of maximum emission width of 
Hé, while the secondary maximum at 64 days corresponds to the 
phase of minimum emission width. 

The widths of the hydrogen emission (Fig. 6) are greater for lines of 
longer wave-length. This fact was noted by Curtiss for Be stars" and 
was later verified by Dustheimer for g Persei. The range of variation 
in width is greater for lines of shorter wave-length. The curves show 
two maxima and two minima in the period of 126.6 days, which was 
also found by Dustheimer from his measures of the one-prism spec- 
trograms. The phase of the first maximum is greater, and of the sec- 
ond maximum is smaller, for lines of shorter wave-length. The 
phases of the two minima are the same for 18 and //6, while for Hy 
the first minimum occurs at an earlier phase and the second mini- 
mum occurs at a later phase. The secondary variation is pronounced 
in H8, is doubtfully indicated in Hy, and appears to be absent in H6. 
This progressive change of the secondary variation with wave- 
length confirms the two-prism microphotometric results discussed 
sarlier in this paper (see Fig. 2). 

The curves for the widths of the violet emission component of /78, 
Hy, and Hé6 are similar in many respects to the curves for the total 
emission widths. The progressive changes in the phase at the first 


1 Pub. Obs. U. of Michigan, 3, 1, 1923. 
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maximum and in the secondary variation are practically identical 
with those noted for total widths. 

The curves for the widths of the red emission component of H8, 
H+, and H6 show in many respects variations just opposite to those 
observed for the widths of the violet emission components. This is 
particularly true for 7. It is apparent, therefore, why the micro- 
photometric results for HB show such great changes in the relative 
intensities of the two emission components. On the other hand, the 
curves for the widths of the violet and red emission components of 
Hé show practically the same variations. These similar variations 
in the width of the emission components of H6 explain why the 
microphotometric results for this line show comparatively small 
variations in the relative intensity. 

At phase go days, when the violet emission edges show a violet 
displacement similar to that of the hydrogen absorption lines, the 
red emission edges show a red displacement; and at phase 30 days, 
when the red emission edges show a red displacement similar to that 
of the hydrogen absorption lines, the violet emission edges show a 
violet displacement. At these phases the widths of the total emis- 
sion lines are maxima. It may be said, therefore, that an emission 
velocity-curve for this star, formed by taking the average of the 
settings on the outer emission edges, has little meaning. Such a 
curve can be interpreted readily only when the two edges show iden- 
tical variations. The values obtained by taking the average of the 
settings on the outer emission edges are plotted in Figure 3 as broken 
curves. As expected, the range of the ‘‘apparent”’ emission velocity 
is small. The emission velocity is roughly zero at the phase where 
the absorption velocity is a maximum. The minimum occurs at the 
same phase as the minimum of the absorption and is approximately 
of the same magnitude. This behavior may be interpreted on the as- 
sumption that when the absorption velocity is a maximum, at phase 
30 days, the central absorption is greatly displaced toward the red 
relative to the emission, thus exhibiting a maximum width for the 
violet emission component, a minimum width for the red emission 
component, and a maximum intensity of the violet emission com- 
ponent relative to the red emission component. When the absorp- 
tion velocity is a minimum, at a phase of about go days, the emission 
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velocity is also a minimum and nearly of the same magnitude, thus 
indicating approximately equal widths for the two emission com- 
ponents, even though the width of the red emission component is 
then at its maximum. However, at this phase the intensity of the 
red emission component is greater than the intensity of the violet 
emission component, instead of being equal to it. Furthermore, the 
intensity of the red emission component is a maximum relative to 
the intensity of the violet emission component at an earlier phase, 
75 days, corresponding to the phase when the width of the emission 
is a minimum and when the width of the violet emission component 
is greater than the width of the red component. It is apparent, there- 
fore, that the intrinsic intensity distribution of the emission lines 
must vary greatly during the star’s period, in addition to variations 
caused by the relative displacement of the absorption with respect 
to the emission. 

It should be added that the foregoing spectral variations of ¢ Per- 
sei occur without conspicuous variations in the visual brightness of 
the star. 

DISCUSSION 

Any system which is proposed for ¢ Persei must account for the 
foregoing spectral variations. Simple orbital motion, as suggested 
by previous investigators, is not capable of explaining the observed 
spectral variations. If it is assumed, however, that the two stars of a 
binary system contribute to the emission of hydrogen, then changes 
in the width and in the intensity of the hydrogen lines would be 
produced by the relative motions of the superimposed lines from the 
two stars. These changes may be enhanced somewhat by the photo- 
graphic effects or shifts discussed by Frédette. Even such a modi- 
fied binary system cannot account adequately for all of the observed 
spectral variations of ¢ Persei. 

Curtiss” suggested radial pulsation for the spectrum variables of 
the y Persei group. Dustheimer’ adopted the pulsation hypothesis 
for g Persei and suggested a secondary pulsation of shorter period 
and smaller range to explain the secondary velocity variation. This 
hypothesis must assume that the inconspicuous light variation ob- 


12 Pub, A.A.S., 34th Session, 1925; Pub. Obs. U. of Michigan, 4, 174, 1932. 
13 Pop. Astr., 35, 484, 1927; dissertation University of Michigan, unpublished. 
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served by Guthnick and Giissow is explained by periodic changes in 
the effective temperature of the photosphere. It does not account 
satisfactorily for the cyclic variations in the widths of the emission 
lines and components in one-half of the star’s velocity period, since 
the range of velocity variation of the emission edges of hydrogen is 
of the same magnitude as the range of the central absorption (see 
Table III). The assumptions of a transparent atmosphere and no ap- 
preciable occultation by the photosphere are therefore not sufficient. 
The chief objection to the pulsation hypothesis is that it requires a 
change in the radius of the star which is greater than the radii of 
normal B-type stars. 
TABLE III 
RANGE OF VELOCITY VARIATION OF THE 
HYDROGEN LINES OF ¢ PERSE! 


| 
He | Hy Hé 





| km/sec | km/sec | km/sec 
Red emission edge | 34 | Gx 61 
Centralabsorption....| 51 | 52 50 
Violet emission edge. | 36 | 51 
| 








The writer suggested, at the time this work was completed in 
1928,"4 a star of very low density, pear-shaped or apioidal in form, 
with a very extensive atmosphere and a small and ill-defined photo- 
sphere, and which is rotating about a line greatly displaced from the 
effective center of the apparent photospheric area with a period of 
126.6 days. It was assumed, further, that the star’s atmosphere is 
pulsating in the same period, that it is contracting while the ap- 
parent photosphere is receding, and that it is expanding while the 
apparent photosphere is approaching; while a secondary pulsation 
was assumed for the secondary variation. Finally, it was assumed 
that the hydrogen absorption originates in the atmosphere between 
the photosphere and the observer, while the hydrogen emission 
originates in the atmosphere projecting beyond the photosphere at 
right angles to the line of sight, as suggested by Rosseland.'’ The 
amplitude of pulsation required in this hypothesis would be much 
smaller than is required by pulsation alone. The arguments ad- 


14 Towa Acad. of Sci. Proc., 36, 281-93, 1929. 5 Ap. J., 63, 218, 1926. 
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vanced in favor of the pulsation hypothesis apply equally well; and 
the effects of rotation, pulsation, and the shape of the star are such 
as to explain the observed variations of g Persei more satisfac- 
torily. 

The most promising hypothesis for explaining the spectral varia- 
tions of ¢ Persei consists in assuming a system of orbital motion in 
which each star is rotating rapidly about its axis, combined with an 
apparent atmospheric pulsation of expansion in the companion star 
at the phase of the secondary variation. Certainly the truly periodic 
variations of g Persei with respect to width and intensity of the hy- 
drogen emission and with respect to radial velocity of the hydrogen 
absorption can be explained satisfactorily by this system. It is only 
necessary to assume that both of the rotating stars of the binary 
system have very extensive atmospheres which contribute to the 
hydrogen emission, the emission originating in the atmospheres of 
the stars at right angle to the line of sight, and that the hydrogen 
absorption originates mainly in the outer atmosphere between the 
observer and the fairly well-defined bright photosphere of the pri- 
mary star. The great widths of the hydrogen emission lines and the 
variation in one-half of the star’s velocity period are explained readi- 
ly by the resultant motion produced by rotation and orbital motion 
in the line of sight. For example, at maximum positive velocity 
(phase 30 days) the primary star is receding in orbital motion. The 
red emission component of hydrogen is greatly displaced to the red, 
owing to rotation, since it originates in the receding atmosphere of 
the primary star, which itself is receding because of orbital motion. 
The violet emission component of hydrogen is displaced toward the 
violet, owing to rotation, since it originates in the approaching atmos- 
phere of the companion star, which itself is approaching because of 
orbital motion. Since the central absorption of hydrogen is produced 
mainly in the atmosphere of the receding primary star, it appears to 
be displaced toward the red with respect to the violet emission compo- 
nent. The violet emission component should therefore show maximum 
width and greatest intensity relative to the red emission component 
at this phase. Furthermore, the width of the total emission should 
be a maximum because of opposite displacements of the two emis- 
sion edges. It should be noted that these observations are in com- 
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plete agreement with the results presented in this paper. Obviously, 
the reverse occurs at the phase of maximum negative velocity; and 
it is to be noted that at this phase the red emission component shows 
maximum width and approximately maximum intensity relative to 
the violet emission component. The width of the total emission 
shows a maximum for the same reason as at maximum positive veloc- 
ity. At the intermediate phases of zero velocity for the absorption 
of hydrogen the widths of the emission components should be mini- 
ma and the relative intensities of the two emission components 
should be equal. Also, the widths of the total emission should be 
minima. Thus, the main and truly periodic variations of ¢ Persei 
can be explained satisfactorily by this binary system. 

McLaughlin” recently suggested a rotating and pulsating mecha- 
nism for Be stars. Dodson‘ has successfully applied this hypothesis 
to explain the irregular cyclic spectral variations of 25 Orionis. There 
seems to be no reason why a pulsation of this nature could not be 
operative in a rotating star of a binary system, such as 6 Lyrae, for 
example. The spectral variations produced by this type of pulsation 
could be quite independent from the truly periodic spectral varia- 
tions caused by orbital motion and axial rotation, and could be merely 
superimposed upon them. The pronounced secondary variations of 
yg Persei, which vary progressively with wave-length and probably 
vary also with time, would appear to be explained by atmospheric 
pulsation of this type. 

It should be noted, however, that immediately after maximum 
positive velocity resulting from the orbital system suggested above, 
the primary star, with the well-defined and bright photosphere, may 
be partially eclipsed by the extensive, rarefied, and transparent at- 
mosphere of the companion star. It is noteworthy that this is the 
phase at which the pronounced secondary variations occur in ¢g Per- 
sei. This condition of partial eclipse could produce spectral varia- 
tions which may be termed “‘secondary variations” and which would 
probably vary progressively with wave-length. At this phase an ap- 
parent atmospheric pulsation of expansion or approach along the 

16 Proc. Nat. Acad., 19, 44, 1933. 


17 Dissertation, University of Michigan, 1934, Ap. J., 84, 180, 1936; abstract in 
Pub. A.A.S., 8, 109, 1935. 
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line of sight may be produced in the intervening extensive, rarefied, 
and transparent atmosphere of the companion star, owing to excita- 
tion by radiation from the hot photosphere of the partially eclipsed 
primary star. Obviously, secondary variations would result also in 
the relative intensity of hydrogen absorption and emission. This 
type of excitation would necessarily be continuous, but its effect 
would be observed only during the partial eclipse. It is probable, 
however, that, owing to rotation, the excitation would be very much 
more effective in the receding atmosphere of the companion star 
at the beginning of the partial eclipse than in the approaching atmos- 
phere of this star toward the end of the eclipse, which is in agreement 
with the fact that the observed secondary velocity variation of hy- 
drogen absorption appears and practically subsides before the ab- 
sorption velocity has decreased to zero. This interpretation of the 
secondary variations of g Persei appears to be the most promising of 


any so far suggested. 


It is with great pleasure that the writer expresses his thanks to the 
late Professor Ralph H. Curtiss for making the investigation possible 
and for his kindly criticism, suggestions, and direction; and to those 
who have helped in procuring the spectrograms. 
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THE SPECTRUM OF H, FROM \X 3612 TO X 4224 
NORTON A. KENT 


ABSTRACT 

The wave-lengths of 512 lines of 72 have been measured on plates taken with a 21- 
foot concave grating at the California Institute of Technology. The results have been 
compared with the measurements by Gale, Monk, and Lee, and by Finkelnburg. 

In 1927 Dr. Millikan suggested to the writer that, although sev- 
eral determinations of the wave-lengths of the lines of the secondary 
spectrum of hydrogen had already been made, further work in this 
direction would be of value. While plates taken with this end in 
view were being measured, the extensive papers of Gale, Monk, and 
Lee,‘ and Finkelnburg? appeared. However, since there were present 
in numerous instances marked wave-length differences between their 
determinations, it seemed advisable to push to its completion this 
present set of values. 

THE APPARATUS EMPLOYED 

The grating used was the 21-foot concave (vertical Rowland 

mount) of the California Institute of Technology. The tube was of 


<2 
ee Se 
x — 
i Lt 
> LO a —— a 
ME — SS \ >a 
- oon eee JO} {+} 
CN l af einen y— 
= o! F / — : 
: Set ia ——_—_— Fiala ae 
oar / 
FIG. 1 


pyrex glass, the capillary, C (see Fig. 1), being 1 cm in internal diam- 
eter and 4o cm long, air-cooled by a fan. A removable quartz win- 
dow, W, permitted the introduction of molybdenum or tungsten (in 
the form of numerous small cylinders) to the capillary, resulting in 
enhancement of the secondary spectrum and partial submergence of 
tAp. J., 67, 89, 1928. 
2 Zs. f. Phys., 52, 27,1929. Attention may be called to the excellent bibliography in 


this article. 
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the primary. Quartz lenses were used throughout. A constant-flow 
method was employed, the hydrogen being generated by a Hoffman 
apparatus and dried by phosphorus pentoxide. The current, given 
by a 10-kw transformer (ratio of transformation 1:74), was always 
under o.6 amp (current density under 0.7 amp per square centi- 
meter), and the pressure was about o.1 mm. 

When the iron spectrum from the Pfund arc, A (length 15 mm and 
current 4~5 amp), alone was desired, a platinum mirror, M, was 
placed as shown. The magnification given by the lens, L,, being 10, 
the image of the arc upon the diaphragm, D, was 15 cm. Of this the 
central part was allowed to pass through a rectangular opening 1 cm 
high. When the tube spectrum alone was desired, M was removed. 
Both spectra could be photographed simultaneously by replacing M 
by a plane quartz plate, P,, and putting another quartz plate, P., 
between the arc and the lens, L,. By rotating P, about a vertical 
axis, the amount of transmitted light could be adjusted. The loss of 
light at P, by reflection, together with that at P, by transmission, so 
reduced the intensity of the arc that correct relative intensities of 
the hydrogen and iron spectra could be obtained, and the arc ex- 
posure could be completely simultaneous with that of the tube. 

Upon each photographic plate three strips were exposed, ob- 
tained by shifting the position of the plate. These contained (A) the 
arc alone, merely for identification, (B) the arc and tube taken simul- 
taneously, and (C) the tube alone. On (A) the exposure time was 
but a few minutes, on (B) some hours, and on (C) many hours—up 
to eleven in some cases. First- and second-order spectra were used. 
The definition on all the plates measured is excellent. The wave- 
lengths of the iron lines employed in calculating the hydrogen wave- 
lengths of strip (B) were the recommended wave-lengths of the 
Transactions of the International Astronomical Union (Volume 3, 
1928). The wave-lengths of the hydrogen lines on strip (C) were 
calculated by means of interferometer wave-lengths* of certain 
standard hydrogen lines obtained by Reginald G. Lacount and 
Robert E. Hodgdon in the Rumford Room of the Massachusetts 
Institute of Technology. These thirty-four wave-lengths are marked 
L in Table II. Their accuracy is attested by the fact that the 


3 Unpublished. 
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maximum probable error of the mean of the determinations as ob- 
tained from three plates is 0.0006 A and the average probable error 
0.0003 A. 

The plates were measured by a Société Génevoise comparator, 
mounted in a thermostatically controlled case. The average varia- 
tion of temperature during the time between settings on standards 
(generally 1~3 cm apart) was under o.3° C. 

In measuring a plate the following procedure was employed: All 
the settings (from two to five, or even ten on some standards) of the 
single-thread cross-hair were made from one side of the spectrum 


TABLE I 


WEIGHTS GIVEN TO VARIOUS DETERMINATIONS 


| Second-Order | First-Order 
Plates Plates 
Intensities o3—6 inclusive 4 2 
Intensities o4 and from 7 to 25 inclusive 2 I 


line. Two distinct runs were made with the plate set red left 
(RL, and RL,), then two with the plate set red right (RR, and RR,). 
The mean of RL, and RR, was calculated, then the mean of RL, 
and RR,. The mean of these two means constituted one determina- 
tion of the wave-length of a line on that strip of the plate. The means 
of the determinations made on either part (B) of the various plates 
referred to the iron standards, or part (C) of the various plates re- 
ferred to the Lacount and Hodgdon hydrogen interferometer stand- 
ards, were then used to obtain the weighted mean value of the wave- 
length, given in column M,, of Table IT. 

As the intensities of a given line on different plates were different, 
the values of the various determinations were weighted as shown in 
Table I. Double weight was given to second-order plates. o4 is the 
intensity assigned to a line just visible. Few lines on the writer’s 
scale were stronger than intensity ro. The illumination used in 
measuring was adjusted to fit the intensity of the line being measured. 

The intensities J recorded in Table II are those given by Gale. 
With these the writer’s are in general agreement; in a few cases, the 
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TABLE II—Continued 
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writer’s intensity is substituted, and where the line was not meas- 
ured by Gale, a value is assigned which agrees with that scale. Fur- 
ther, in Table II, \ is the integer value of the wave-length in ang- 
stroms; F, Finkelnburg’s fractional value; G, that obtained by Gale, 
Monk, and Lee; L, the Lacount interferometer value; M,,, the writ- 
er’s weighted mean determination; NV, the number of ‘‘sets”’ of meas- 
urements entering into the calculation of M,,; PE,‘ the probable 
error of the weighted mean, M,,; and v, the “vacuum” wave-numbers 
obtained from M,, or from L, where the latter is present. 

With the interferometrically determined wave-lengths, L, the 
grating determinations, M,, (shown in parentheses and based on the 
iron standards), are listed for comparison. The maximum deviation 
of M,, from L is 0.006, the average 0.002 A. 

The average probable error of the mean for the values of M,, for 
the 512 lines listed is less than 0.002 A. 

In three cases Lacount and Hodgdon measured lines determined 
interferometrically by Gale, Monk and Lee. These are \ 4171.308, 
L = .307; 4 4195.674, L = .670; \ 4212.498, L = .500. They are 
marked by an asterisk in column G of Table II. 


DISCUSSION OF RESULTS 


The writer’s values agree, in general, more closely with those of 
Gale. The average difference between G and M,, without regard to 
sign is o.o11 A; that between F and M,, 0.022 A. 

A comparison of the 34 interferometer values, Z, with the mean 
M.,, reveals that the sum of the positive values of Z — M,, is .034 A 
and the sum of the negative is also .034, showing an average devia- 
tion of .cor in each case, which indicates a random distribution of 
differences and attests the accuracy of the entire set of measure- 
ments. For L — G the sums are + .127 and — .o51; for L — F, 
+ .176 and — .047. 

It is the writer’s purpose to publish the wave-lengths of lines in 
other regions as soon as possible. 

4 It should be borne in mind that the probable error of a mean is of small significance 
in cases where few determinations are involved. 


> Calculated from Kayser’s table. 
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PHOTOMETRY OF THE DIFFUSE NEBULA NGC 7023 


PHILIP C. KEENAN 


ABSTRACT 

The reflection nebula NGC 7023 was photographed in red and in blue light with the 
Perkins reflector. The photographic surface brightness is inversely proportional to the 
distance from the central star for radii of the order of 1’, but the curve is steeper for 
smaller radii. The observations do not fit any of the predicted curves computed on 
the assumption of uniform density, but indicate that there is a concentration of mat- 
ter about the star, which probably lies a relatively short distance inside the nebula. 

The color of some of the brighter parts of the nebula was compared with that of the 
star. For the larger part of the nebula the reflecting power is only slightly greater in 
red light than in blue, as would be expected if the particles are fairly large; but there is 
a well-defined portion in which a red excess of the order of 0.3 mag. was found. 

NGC 7023 in Cepheus is one of the brightest of those diffuse nebu- 
lae which are known to have a predominantly continuous spectrum,’ 
presumably due to light from the central star reflected by the nebu- 
lar matter. Because of its comparative symmetry, in contrast to 
several other objects of its class, as well as its brightness, it was se- 


lected for a reconnaissance of the intensity gradient and the color of 
the nebulosity. 
METHOD OF OBSERVATION AND REDUCTION 

The observations were carried out at the 25-foot Newtonian focus 
of the Perkins reflector, which was diaphragmed down to an aper- 
ture of 60 inches in order to eliminate vignetting arising from the 
smallness of the secondary mirror. Exposures were made on Eclipse 
850 Soft plates to give photographic magnitudes and on Eastman 
1 C plates through a ciné-red filter as used by C. Payne Gaposchkin 
and S. Gaposchkin? to give red magnitudes with an effective wave- 
length of 6300 A. All plates were backed and were provided with in- 
tensity-calibration exposures made through a tube sensitometer. 

Measurements of image density were made for the most part with 
a Kipp registering microphotometer in the Mendenhall Laboratory, 
but some of the supplementary plates were measured on a Stetson 
photometer at the Perkins Observatory. 


*In NGC 7023 the hydrogen lines have been observed to be bright in both nebula 
and star. Cf. V. M. Slipher, Pub. A.S.P., 30, 63, 1918. 


2 Harvard Ann., 89 (No. 5), 93, 1935- 
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The central star of the nebula, BD+67° 1283, of photographic 
magnitude 7.51 and spectrum B2p, was chosen as the standard for 
photographic brightness and relative color. On each plate to be used 
for color determination the star was photographed 12 mm out of 
focus both before and after the exposure on the nebula. Between 
exposures the plate was shifted so as to place the annular star image 
outside the measurable limits of the nebula but on the same back- 
ground of sky fog, in order to eliminate, as far as possible, any effect 
of the sky background on the measures. 

By using an exposure ratio, nebula to star, of 30:1, the density of 
the inner nebula was matched satisfactorily with that of the stellar 
image. This procedure required that the value of Schwarzschild’s 
constant be known for both types of emulsion and for the given 
exposure times. On the basis of a laboratory determination checked 
by exposures on several Selected Areas, the values of p adopted were 
0.83 for the Eclipse plates and 0.74 for the 1 C plates with the red 
filter, but the uncertainty amounts to several per cent. 

Comparison of the densities of the nebular and stellar images in 
each color gave directly the relative reflecting power of the nebulos- 
ity in red and blue light. The color determined in this way is ex- 
pressed in terms of that of the central star. The absolute photo- 
graphic surface brightness of the nebula was then found by measur- 
ing the area of the annular image of the star, which was always pho- 
tographed on the optical axis of the telescope in order to secure a 
circular ring of sufficiently uniform density. 

The observations were greatly facilitated by a rack-and-pinion 
focusing device and a removable filter-holder, which were construct- 
ed in the Observatory shops by Mr. Decker and mounted on the 
double-slide plate-carrier. 


OBSERVATIONAL RESULTS 


Although NGC 7023 possesses so many irregularities of structure 
(Pl. XVIII) that measures along any one radius have little signifi- 
cance, it was found that the means from several diameters followed a 
comparatively smooth curve which gave a fair representation of the 
gradient of photographic brightness in passing outward from the cen- 
tral star. The values obtained by combining data from eight radii on 
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each of three of the best plates taken here are plotted in Figure 1. 
The abscissae in the diagram give the log,, of the distance from the 


Am 











l l l I Lil 


—0.4 —0.2 ° 0.2 0.4 log s 











—o0.8 0.6 

Fic. 1.—Relative intensity gradients in reflection nebulae. The photographic ob- 
servations on NGC 7023 are represented by the black dots. 

central star in minutes of arc, The ordinates represent relative 

brightness in magnitudes, the zero-point being placed arbitrarily at 
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a distance of 1 minute. The significance of the curves in the figure 
will be discussed in the next section. 

From comparison with the central star on three additional plates 
the actual surface brightness at eleven points in the inner nebula 
was determined. These data led to a value of 11.61 + 0.03 mag. as 
the factor to be added to the relative magnitudes of Figure 1 to ob- 
tain photographic surface brightness in magnitudes per square min- 
ute of arc. Thus the luminosity of the nebula decreases from 10.6 
mag./’ at a radius of 0/5 to 12.8 mag./ 0’ at a radius of 3/2. A 
few rough comparisons with stars of Selected Area 18 confirmed the 
absence of any large systematic error in these values. 


TABLE I 


+o.19 mag. +o.23 mag. 
— .19 + .26 
+ .24 + .06 
+0.02 —0.23 


With the inclusion of data from three red-sensitive plates it was 
possible to derive colors for the brighter parts of the nebula along 
three radii. Expressed as the difference between the brightness of 
the nebula relative to the central star in red light and in blue light, 
the color differences for eight points of measurement are given in 
Table I, where positive values indicate red excess. 

The large scatter is due in part to uncertainties in the measures 
and in part to variations of color within the nebula. The average of 
the measures is +0.09 mag. The nebula is thus nearly white, but 
appears to have a slight increase in reflecting power toward the red 
end of the spectrum. The general similarity in the structural details 
visible in the two colors is apparent in Plate XVIII, but it may be no- 
ticed that there is a wedge-shaped dark region about 40” northwest of 
the center that is crossed by a curved bright filament which is much 
more conspicuous on the red plate than in blue light. Three addi- 
tional points measured within this region show a strong red excess, 
the values being +0.30, +0.39, and +0.55 mag. It is evident that 
either there is selective absorption amounting to about 0.3 mag. in 
this region or there is an excess of erispsion in red wave-lengths, pos- 
sibly the result of local excitation giving rise to emission lines. The 
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former explanation appears the more plausible, but observations 
with a fast spectrograph are desirable to settle the point. 

By increasing the red exposures to more than three hours with the 
fastest 1 C Special plates, made available through the co-operation 
of Dr. C. J. Staud of the Eastman Kodak Company, it was possible 
to extend the measures far enough from the center to indicate that 
the mean intensity gradient in red light is similar to that found for 
the blue; but an accurate color chart of the whole nebula cannot be 
made until plates of still greater sensitivity in the red are developed. 

Nebulae such as NGC 7023, which are nearly as blue as the early- 
type stars which they surround, are particularly difficult to photo- 
graph through a red filter. On the other hand, the emission-line 
nebulae have a greater proportion of their light in the visual region. 

A good example is the Swan nebula NGC 6618 (M17), which is 
reproduced in Plate XIX from a one-hour exposure through the red 
filter. If this picture is compared with that taken by J. C. Duncan: 
on an ordinary blue-sensitive plate with the 60-inch reflector at 
Mount Wilson, it will be seen that practically every detail is repro- 
duced with the same relative contrast on the two plates, from which 
it seems reasonable to conclude that there are no very great differ- 
ences in excitation within the nebula. The conspicuous rectangular 
dark marking near the center is as black in reality as it appears to be 
in contrast to the brilliant surrounding regions, for measures on our 
plate showed that it had the same density as the sky background 
near the edge of the field. 

The star indicated by the white arrow in the right-hand image of 
Plate XVIII is of interest as having an unusually low color index for 
an object of its luminosity. From our plates its color can be estimated 
as somewhere between those of an A5 and a Gz star. Since it lies 
within the area of the dark nebula which covers the whole field of 
the plate and is classified as very opaque by Struve and Miss Story,‘ 
it is almost certainly situated between the sun and the nebula. The 
distance of the nebula is uncertain; but if, for the moment, we ac- 
cept Lundmark’s estimate of about 200 parsecs,’ the apparent photo- 
graphic magnitude of 15.5 for the star leads to a value of +8 or +9 


3 Ap. J., §1, 4, 1920, Pl. IV. 
4 Ibid., 84, 203, 1936. 5 Pub. A.S.P., 34, 40, 1922. 








PLATE XIX 








THE SWAN NEBULA, M. 17 (NGC 6618) 


E;xposure one hour through red filter 
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as its maximum absolute magnitude. It would thus appear to be 
intermediate between the main-sequence stars and the white dwarfs, 
and it deserves further attention. 
THEORETICAL INTENSITY GRADIENTS 

The theory of the illumination of small particles by a star was de- 
veloped by Seeliger® and elaborated by Schoenberg’ and by Struve.® 
Its adaptation to the conditions of our present problem is simple. 
We assume that the density throughout the nebula is constant. 

















Te Earth 
a b 


Fic. 2.—Position of star with respect to nebula 


Consider first the case of a star imbedded in a very large nebula 
(Fig. 2a). Then, in comparing the brightness of the nebula with that 
of the star, we can neglect the absorption of the light after reflec- 
tion, for, on the average, this will be very nearly the same as the 
absorption of the direct beam from the star to the earth. At any 
point P the illumination of a particle will be 


, pakr 
E, = = ; (1) 
where £, is the flux over the same cross-sectional area at unit dis- 
tance from the star and & is the linear absorption coefficient. 
The radiation scattered in the direction of the earth is given by 
6 Sitz. d. Akad. d. Wiss. (Miinchen), Heft III, p. 272, 1901. 
7 Handb. d. Astroph., 2, Part I, 163, 1928; Mitt. Universitdtssternwarte Breslau, 3, 


53, 1932. 


§ Struve and Elvey, Ap. J., 83, 162, 1936; Struve and Story, op. cit. 
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KO®(a)E,, where K is a constant proportional to the albedo of the 
particle and O(a) is the factor giving the distribution of the scat- 
tered light as a function of the phase angle a. For a line of sight 
passing through the nebula at apparent distance s from,the star, the 
total intensity of the scattered light is then given by 


*CO eke 

I, = Kn a f ie O(a)dx (r = Vs? + 2°), (2) 
—oO 

n being the number of particles per cm. 

Since the factors outside the integral are constant throughout the 
nebula by our assumptions, the apparent intensity gradient can be 
found by evaluating the integral as a function of s. 

An alternative form of equation (2) is 


n e ks csc a 
I, = Kn Bf S O(a)da, (3) 
° $ 


from which it may be seen that, if there is no appreciable absorption 
and the phase-angle factor can be neglected, the brightness of the 
nebula is inversely proportional to the apparent distance from the 
star. This case is represented by the straight line A in Figure r. 

For the more general case equation (2) was integrated graphically 
for 
sina + (mr —a)cosa 

T 





O(a) ~ 


(large particles subject to Lambert’s law?) and for O(a) ~ 1 + cos?a 
(Rayleigh or Thomson scattering) as well as for O(a) = 1. As ex- 
pected, the phase function was found to affect the curves so slightly 
that it can be neglected in comparison with other sources of error. 
Consequently, only the curve for O(a) = 1 has been reproduced 
(curve B). For small radii the slope of the curve approaches that of 
the line A. 

If the star is at the near edge of the nebula, equation (2) becomes 


~” ek(r+z2) 
I, = KnE, i O(a)dx , (4) 
° 


9 Schoenberg, Handb. d. Astroph., 2, Part I, 64, 1928. 
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since it is now necessary to take account of absorption after reflec- 
tion. This is essentially Seeliger’s expression.’° As shown by curve 
C, the gradient is now somewhat steeper at any given radius. There 
will be a gradual transition from C to B as the star is taken inside 
the nebula. 

If the star lies in front of the nebula, we have (Fig. 2b): 


-k(r’+2—n) 
I, = KnB f° = O(a)dx , (5) 


where 


, “(xe — n) 
r= ——.. 
% 


As shown by the curve D, computed for m = 1, the gradient flat- 
tens out very rapidly near the star if there is any considerable dis- 
tance between star and nebula. 

In all of the plotted curves the zero-point for the ordinates was 
chosen to make the curves pass through the point (0, 0). As a matter 
of convenience, k was chosen as 0.6. Actually the abscissae are 
fixed by the opacity ks; and when the actual distance of any nebula 
is known, the absorption coefficient may be determined approximate- 
ly by sliding the curves along the log s-axis until one of them fits the 
slope of the observed curve at corresponding abscissae. 

In the case of NGC 7023 the observations cannot be fitted well to 
any of the theoretical curves. The observed gradient is nearly pro- 
portional to the inverse radius over most of the curve but becomes 
steeper near the star. This increase in slope is probably real, for the 
observations have been corrected for atmospheric and photographic 
scattering of light from the central star. Probably the discrepancy 
is due to the fact that our assumption of uniform density fails to 
hold for this nebula; the appearance of the photographs suggests 
strongly that there is an abrupt condensation of matter in the imme- 
diate vicinity of the star. Whether this is accidental or is due to 
gravitational influence of the star can best be settled by similar ob- 
servations on other reflection nebulae. Visual examination of plates 
taken here suggests that NGC 2023 and NGC 2071 do not show a 


10 [bid., p. 165, eq. 77. 
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correspondingly rapid increase in brightness toward the star, but 
NGC 2247 appears quite similar to NGC 7023. 

From the observations presented here, it can be concluded tenta- 
tively that NGC 7023 consists of particles having diameters too 
great to scatter according to Rayleigh’s \~‘ law, the concentration 
being so low as to cause very little absorption in the light from the 
star at distances of the order of 10,000 to 30,000 astronomical units 
but increasing rapidly at smaller distances. The star probably lies 
within the nebula but not far from the nearer edge. 


I am greatly indebted to Dr. Struve, director of the Yerkes Ob- 
servatory, for his suggestion of the importance of this observational 
problem and for valuable discussions in the course of the work. I 
should like also to acknowledge the help of Mrs. J. A. Hynek in 
carrying out the computations. 


PERKINS OBSERVATORY 
September 1936 








ON THE POLARIZATION OF LIGHT IN 
REFLECTION NEBULAE 


L. G. HENYEY 


ABSTRACT 

Observations for the purpose of detecting the presence of polarization in the light 
of the reflection nebula NGC 7023, the nebulae in the Pleiades and the nebulous patches 
associated with y Cygni have been made through a Polaroid screen used as an analyzer. 
Each set of observations consisted of two plates having the principal plane of the Po- 
laroid of one plate perpendicular to.that of the other. The relative weakening or 
strengthening of the image along radii drawn from the star parallel to the positions of 
the principal plane of the Polaroid would indicate polarization. 

Measurements on the plates of NGC 7023 indicate the presence of polarization up to 
12 per cent. The plates of the Pleiades and of y Cygni indicate little or no polarization. 

These results are shown to be compatible with the reflection theory of diffuse 
nebulae, provided the reflecting particles are fairly large. In the case of the nebulae 
around y Cygni, the lack of evidence of polarization does not favor the hypothesis of 
the presence of Rayleigh scattering in these objects. This hypothesis was suggested by 
the blue color of these nebulosities. 

1. The investigation of polarization in reflection nebulae which is 
reported in the following article was undertaken to provide informa- 
tion concerning the nature of the nebular particles which scatter the 
light of the associated stars. Certain definite deductions are possible 
from a knowledge of the state of polarization, while certain other 
facts can be elucidated by combining the results of observations of 
polarization with results of other methods of attack. 

The existence of polarization in any amount indicates that the 
scattering particles are not free atoms, for processes of emission all 
give rise to unpolarized light, except in the case of the scattering of 
light in a resonance line. But the observations of Meyer’ on emission 
nebulae show that these objects exhibit no appreciable polarization. 

If we know beforehand that a nebula is not of the emission type, 
observations of its polarization enable us to go a step farther. In 
general, if a particle scatters light in such a way that it does not hold 
the energy for any length of time but simply deflects it without 
change of wave-length, then the polarization of the deflected light 
has the following characteristics: (a) its plane of polarization is 
usually perpendicular to or, more rarely, parallel to the plane formed 

t Lick Obs. Bull., 10, 68, 1920. 
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by the incident and the scattered rays; and (0) the amount of 
polarization is inversely correlated with the size of the particles. 
When the scattering particles are of the order of size of a few hundred 
molecules or smaller we have what is usually known as Rayleigh 
scattering. Here the polarization reaches very large values, as is 
evidenced by our own blue atmosphere, which gives values up to 70 
per cent. For larger particles the polarization is much smaller and 
is quite negligible, except for the case of specular reflection near the 
polarizing angle on the surfaces of very large bodies. 

The orientation of the plane of polarization should give an inde- 
pendent confirmation of the fact that reflection nebulae are illumi- 
nated by stars. For if the nebula is sufficiently spread out, the plane 
of polarization will have different orientations at different points, 
depending on the direction of the line joining the star and the point 
in question. Consequently, if we observe polarization in a nebula 
such that its plane is symmetrical with respect to a star, we can safe- 
ly say that this star is the source of illumination. It must be remem- 
bered that these relationships, strictly speaking, hold true only if we 
consider light which has suffered but one deflection after leaving the 
star; for light which has been deflected more than once the case is 
more complicated. However, in almost any case the planes of 
polarization of multiply scattered light will be so thoroughly mixed 
that all trace of polarization is lost. It therefore seems that radial 
polarization should be present. 

The distribution of polarization over the apparent disk of a 
nebula could, under favorable circumstances, give additional infor- 
mation concerning the structure of the nebula. Here, unfortunately, 
the situation is much more difficult to analyze, since irregularities 
in the space density of the particles in the nebula can influence the 
result tremendously. But it may be well to keep this source of at- 
tack in mind, since undoubtedly in the future it will provide valuable 
information. 

While several attempts have previously been made to detect 
traces of radial polarization in nebulae, no quantitative measures 
have been published. Most nebulae are faint and visual observa- 
tions are therefore out of the question. Previous attempts to make 
observations, together with the difficulties encountered in the con- 
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struction of suitable analyzers, have been described by Dr. W. F. 
Meyer.’ He found that the best results were obtained with a Nicol 
prism. The nebula was photographed in two perpendicular positions 
of the Nicol and the photographic densities of the two images were 
measured with a microphotometer. A comparison of the two sets of 
measures was then made to detect traces of polarization. Owing to 
the large optical thickness of the Nicol, it was necessary to collimate 
the light of the nebula before it passed through the prism. A camera 
lens on the other side of the Nicol focused the light on the photo- 
graphic plate. With this arrangement Meyer obtained a series of ob- 
servations of emission nebulae, of a number of spirals, and of one 
true reflection nebula, NGC 2261 (Hubble’s Variable Nebula), in 
which we are here primarily interested. He found that in this object 
the polarization is small, probably less than 1o per cent. More re- 
cently Struve, Elvey, and Roach? found that the nebula near 
p Ophiuchi exhibits a definite but rather small amount of polariza- 
tion, and they pointed out that this is compatible with the reflection 
theory. Since the effect is radial with respect to p Ophiuchi (unless 
the principal plane of the Polaroid was taken by chance in just the 
right direction), we have a confirmation of the association of the 
star and the nebula. 

2. The detection of a state of radially symmetrical polarization can 
be very simply accomplished by the use of the Polaroid screen. Since 
the Polaroid is thin, it is unnecessary to collimate light which passes 
through it. The Polaroid is primarily a film of some celluloid-like 
material into which are imbedded numerous minute crystals of a 
double-refracting material. The optic axes of these crystals have 
been oriented into parallel directions, so that the film will polarize 
light which passes through it. This polarizing property is not com- 
plete and depends on the wave-length.* It is sufficient to state here 
that the region of transmission is toward the longer wave-lengths 
from 4200 A, and that the region of good polarization (better than 
80 per cent) is between this limit and about 7200 A. The average 
polarization between these limits will depend on the color character- 

2 Thid. 3 Ap. J., 84, 219, 1936. 


4 A description of the transmission and polarization of the Polaroid can be found in a 
paper by Ingersoll, Winans, and Krause, J. Opt. Soc. Amer., 26, 233, 1936. 
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istics of the source and of the photographic plate, but it will usually 
be over go per cent. In the discussion which follows it will be as- 
sumed that the polarization of the Polaroid is 100 per cent. 

On long exposures of bright stars, peculiar streaks were noticed 
passing through the images. Each streak is along a line parallel to 
the principal plane of the Polaroid, and decreases in intensity with 
increasing distance from the star, the total intensity depending on 
the density of the stellar image. The source of the effect is undoubt- 
edly in the Polaroid and may be due to the scattering of the light 
in the polarized component removed from the principal beam. It 
was only on the plates of the Pleiades that the effect was noticeable, 
but even here little difficulty was experienced. It should be remarked 
that the effect is one which is in the opposite sense from the effect 
resulting from radial polarization. 

The use of the Polaroid in detecting radial polarization is illus- 
trated in Plate XXa. This plate shows two photographs, each taken 
through a Polaroid film, of a rectangular glass tank full of a colloidal 
dispersion of rosin in water which is illuminated by a light placed in 
the center. The center is greatly overexposed, so that the lamp is 
not visible. The lines represent the directions of the principal planes 
of the Polaroid in the two cases. Now the light scattered by such a 
medium is polarized in a plane perpendicular to the radius vector, 
so that the Polaroid will weaken the light to the greatest extent along 
the radii parallel to the principal planes. This is readily seen by com- 
paring the two cases. A pair of such photographs of a nebula will 
then provide the means of detecting the presence of polarization. 

3. In Figure 1 let O be a star illuminating a nebula surrounding it, 
and choose two points, P,; and P., lying on radii which are perpen- 
dicular to each other. Let the surface brightness at the two points 
be J, and J,, respectively, and the fraction of the light polarized be 
p, and p, in a plane perpendicular to the radius. Then the light in- 
tensity at each point can be resolved into two components, one 
along the radius and one along the tangent. The two components 
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a) POLARIZATION IN A COLLOIDAL DISPERSION OF ROSIN IN WATER 


b) POLARIZATION IN NGC 7023 
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and for P, corresponding values. Now if we take two photographs 
of such a nebula through a Polaroid disk such that the principal 
plane in one is along OP, and in the other along OP., we will have 





FIG. 1 


but one of the components passing through. The transmitted light 
in each case is given below: 
Case A—Principal plane of Polaroid along OP;: 


ize Siew 


2 


Case B—Principal plane of Polaroid along OP.: 


I - Pr im p2 
JB —_ Sy ’ J.B _ i 
2 2 
It is convenient to express surface brightnesses in terms of the 


sky as a unit. Now on different plates the sky will be somewhat 
different, so the intensities on different plates are not directly com- 
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parable. However, if we take the ratios of intensities on the two 


plates for both points, we have 


Jia Fea 
Jayskky _U©— pr JB,sky Ja,sky _I+ f2 JB,sky ; 
Jie=SCUI +, Jajsky’ Joa I-f, ee (1) 
JB, sky JB, sky 


Taking the ratio of these two quantities, we have 


JA S24 

Jaysky , JaAysky _ (1 — p1)(1 — pz) 

Jip JapsC(t EH P(A + 2)” (2) 
Jz, sky Jp, sky 


The left-hand side consists of known quantities, so we can compute 
the right-hand side. If we assume that the polarization is the same 
at the two points, we determine the quantity 
few] 
i+ a : 
from which p can be computed. In any case, although p, and p, are 
different, p represents an average value. 

4. The objects which have been examined for polarization are 
NGC 7023, the nebulae in the Pleiades and the nebulae around y 
Cygni. 

The plates of NGC 7023 have been obtained with the 24-inch 
reflector of the Yerkes Observatory, on Eastman Hyper-press plates. 
Each plate has been standardized with a tube sensitometer. ‘The 
Polaroid used was a piece mounted in a bakelite rim having a clear 
aperture of 4 cm. It was mounted on the plateholder carrier of the 
reflector, a few millimeters in front of the plate, with the principal 
plane always parallel to the short dimension of the plate. Since the 
carrier may be rotated about the optical axis, the adjustment of the 
orientation could be made by rotating the carrier any amount de- 
sired, according to the divided circle on it. However, only two posi- 
tions were used: the principal plane parallel and perpendicular to 
the hour circle through the nebula. 

Two photographs are reproduced in Plate XXb. The straight 
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lines give the direction of the principal plane of the Polaroid. It is 
seen that an effect is present which is somewhat similar to the effect 
shown by the colloidal solution. The relative weakening of points 
along the two radii is, for points which show maximum effects, 
about 0.3 mag. 
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Fic. 2.—Microphotometer tracings of NGC 7023 


The plates have been measured in two independent ways which 
give the same results. Microphotometer tracings have been made 
along the two radii. These are reproduced in Figure 2, where the 
process of reduction is indicated. The average density along two 
regions on the perpendicular radii was determined by drawing a 
straight line through the region and taking the density of the middle 
point on the line. The difference in magnitude between this point 
and the sky is shown in the figure. Measurements were also made 
at selected points along these radii by means of a Hartmann micro- 
photometer. These measures are collected in Table I. The upper 
part of the table gives the results for one plate, the lower those for 
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a second plate. The principal plane of the Polaroid is indicated in 
the first column. The measures on the left are along the north- 
south radii, while those on the right are along the east-west radii. 
The table gives the measured difference in magnitude between the 
points on the nebula and on the sky, the corresponding intensity 
(plus that of the sky) in terms of the intensity of the sky, and the 
intensities corrected for the sky. The next to the bottom row gives 
the quantity involved in equation (1)—namely, the ratio of inten- 
sities of points on plate A to those on plate B. In the bottom row 
is the mean for the three sets of measures along each radius. Taking 




















TABLE I 
Plate North-South Radii East-West Radii 
Pe Oe —m}) 1.54 | 1.41 | 1.62 | 0.68 | 0.58 70 
E+? | 4:23.) 3:66) 4.48 | TiS | £2.90 || 2:07 
Lt gata | 2.66: | 3.48 1-06-87 0.9 | 1.07 
B. Polaroid EW. pone —m|1.70 | 1.47 | 1.80 | 0.55 | 0.42 | 0.65 
I+1 | 4.79 | 3.87 | 5.25 | 1.66 | 1.47 | 1.82 
I | 3.79 | 2.87 | 4 0.66 | 0.47 | 0.82 
LS eee care a ees emert 0.85 | 0.93 | 0.81 | 1.32 | 1.51 | 1.30 
EE a ee ee eee me es, Server ie) eee et aaa Mean OOO occ ccalnns ce 1.38 

















the ratio of these values, we obtain for these points the mean of the 


quantity | 
a 


= 12 percent. 
I 


from which we find 


It should be repeated that this result is an average for those points 
which show large amounts of differences from plate to plate. Conse- 
quently, the average for the whole nebula is probably somewhat 
smaller. 

The plates of the Pleiades were taken in the same manner as those 
of NGC 7023, except that Eastman Ortho-press plates were used. 
The exposures were approximately 33 hours. The length of the ex- 
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posures and the brightness of some of the stars caused the appear- 
ance of the anomalous streaks already mentioned. However, the 
principal features of several nebulae were well marked, and com- 
parisons could be made satisfactorily. The reader is referred to 
Plate 4 of Barnard’s Atlas of Selected Regions of the Milky Way, on 
which the features described in the following discussion can be easily 
identified. 

Surrounding Maia, there are a number of faint wisps of nebulosity 
and one very well-marked tonguelike patch curving from west of 
the star around the northwest and gradually fading away toward 
the north. Considerable weight can be given to an examination of 
this feature because of its large spread in position angle. No effect 
of polarization was found in it. 

The results for Merope are, similarly, negative. Here too, there 
is a well-defined curved patch, running from the east to the north. 
No effect was noticed along the points of this nebulosity. A compari- 
son of its eastern portion and of other patches nearer the star with 
the extended cloud south of the star also failed to reveal a relative 
change. 

A small trace of polarization was found in the approximately 
straight filament running eastward from Electra. This filament was 
compared with another narrow and straight filament southeast of 
Electra and apparently running through a faint star (BD + 23°510) 
from east to west. It was found that, when the Polaroid is parallel 
to the first filament, it is weakened with respect to the second. A 
comparison with the plate as a whole showed the same effect. It 
seems, therefore, that Electra provides the illumination for this neb- 
ula. The magnitude of the effect, even in this nebulous filament, is 
definitely less than that found in NGC 7023. 

The plates of y Cygni were taken with a Schmidt camera having a 
focal ratio of one to two. The camera has been described in the paper 
by Struve, Elvey, and Roach already referred to. The observing pro- 
cedure is the same as that used by these observers. A visual exam- 
ination of these plates reveals that the polarization, if present, is 
exceedingly weak. Owing to the small scale of the plates and the 
closeness of the star images, it did not seem advisable to measure 


these plates. 
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5. The results given in the foregoing section indicate that the 
polarization in these nebulae is quite weak. It is considerably 
weaker than the polarization of the corona,’ for which the maximum 
is about 27 per cent. The polarization in NGC 7023 is about of the 
same order of magnitude as that found in the zodiacal light, namely, 
15 per cent.° This small amount is in accordance with the reflection 
theory. Scattering by very small particles, such as Rayleigh scatter- 
ing, would, for a probable distribution of the particles of the nebula 
in space, give values of the polarization much larger than this 
perhaps between 30 and 70 per cent. The small amount of polariza- 
tion in NGC 7023, in the nebulae of the Pleiades and around p Ophi- 
uchi justifies the conclusion that the scattering particles are not of 
molecular size, or of a size comparable with, or smaller than, the 
wave-length of light. Since color observations of these nebulae have 
already shown’ that they are not appreciably bluer than the illumi- 
nating stars, it is permissible to assume that the nebular particles 
are fairly large in size and that they scatter the light of the neighbor- 
ing stars in accordance with the theory of Seeliger. A discussion of 
this theory will be given in another paper. 

If Rayleigh scattering is indeed ruled out for the nebulosity as- 
sociated with y Cygni, then we should not expect selective scatter- 
ing to take place. This point seems to be at variance with the color 
index for this object.’ An explanation of this discrepancy may be 
found in the spectrum of the nebulosity. If there are any emission 
lines, they may give rise to such a spurious result. On the other 
hand, the presence of emission lines would not be expected from the 
spectrum of the exciting star, which is F8p, unless there are peculiar 
conditions of excitation. Further observations are needed before 
this point can be cleared up. 

YERKES OBSERVATORY 
October 1936 

5 Dufay and Grouiller (C.R., 196, 1574, 1933) find 26 per cent 10’ from the solar 
limb, while Johnson (Pub. A.S.P., 46, 226, 1934) finds 28 per cent 8’.5 from the limb. 

6 Dufay, C.R., 181, 399, 1925. 


7 Struve, Elvey, and Keenan, Ap. J., 77, 274, 1933; Struve, Elvey, and Roach, ibid., 
84, 219, 1936; and Keenan, ibid., 84, 600, 1936. 


§ Struve, Elvey, and Roach, op. cit., p. 226. 
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Parallaxes with the 60- and 1oo-Inch Reflectors. Fifteenth Series 

WuitTForD, A. E., and Joet Stepsrins. Absorption and Space Reddening 
in the Galaxy from the Colors of Globular Clusters 

WiLpt, Rupert. Low-Dispersion Spectra of Red Stars 

Wiis, Emma T. R. The Change of the Solar Velocity with henteadon 
Distances of the Stars from the Galactic Plane 

WItson, O. C. Variations in the Spectrum of P Cygni 
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A Wealth of Illustrative Material 
for Teachers, Lecturers, and Students of ASTRONOMY 


800 LANTERN SLIDES, prints, and transparencies are now available 
from Yerkes Observatory (29 slides are new this year). This valuable ma- 
terial is the work of such men as Ross, Ellerman, Barrett, Slocum, Frost, 
Adams, Pettit, Fox, Lee, Barnard, Hale, Ritchey, and others. A complete, 
40-page illustrated catalogue will be sent on request. 


A “100 LIST” of these slides has been assembled especially for schools 


wishing a representative presentation of astronomical material. This excel- 
lent set is $62.50, a saving of #12.50 since individual slides are 75 cents. Card 
descriptions are $2.50 a set. 


4. MOTION PICTURES from McMath-Hulbert 


(1 Real, 16 mm., silent) 


A Motion Picture Journey to the Moon 
Sunrise and sunset phenomena in 
Lunar Craters are shown and an occul- 
tation of Delta Capricorni by the 
Moon. 

The Solar Eclipse of August 31, 1932 
Shows the actual progress of the Moon 
across the Sun. 

Jupiter 
Shows the rotation of Jupiter on its 
axis....nebulae....star clusters 
.... and star fields. 

Solar Phenomena 


Analyzes the light of the sun by means 
of the sun spectrum motion writer and 
the phenomena are pictured in the 
light of ionized calcium only. 


Sales price for film . $30.00 each 


Rental for Fo. csci 0%. 
Pave $2.00 each plus carriage 


4, NEW TALKING MOTION PICTURES 


Produced under the supervision of Walter Bartky. Ready fall, 1936. A dis- 
tinctive contribution to educational methods, the films are designed for use 
in the University’s Physical Science Division. 


THE EARTHIN MOTION THESOLARFAMILY THEMOON THE GALAXY 






Sotar Ecuipse or Ava. $1, 1932 


For complete information and descriptive catalogues write: 


THE UNIVERSITY OF CHICAGO PRESS 


5750 Ellis Avenue 
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THE OBSERVATORY 


A Monthly Review of Astronomy, Founded 1877 





Contains full Reports (including the speakers’ accounts of their work 
and the discussions which follow) of the Meetings of the Royal Astro- 
nomical Society and of the Meetings for Geophysical Discussion: 
Articles: Reviews of important astronomical books: Correspondence 
on topics of interest: Notes on current discoveries and research, etc. 





Annual Subscription: 20/— for 12 Numbers. 





The Observatory, founded in 1877 by Srr Wa. Curistie, Astronomer Royal, has been 
edited in the past by W. H. M. Curistiz, E. W. Maunper, A. M. W. Downie, T. 
Lewis, H. H. Turner, A. A. Common, H. P. Horus, S. Caapman, A. S. Epprncron, 
H. Spencer Jones, F. J. M. Srratron, J. Jackson, W. M. H. Greaves, J. A. CARROL, 
W. H. STEAVENSON and R. O. REpMAN. 

Address of Editors: Royal Observatory, Greenwich, England. 
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Annual Tables of Constants AT C 


and Numerical Data 


A NEW VOLUME 
of numerical documentation of SPECTROSCOPY 


To the Readers of the Astrophysical Journal.— 

This new Volume is a reprint of Volume X (1930) of the ATC, All explanations of the Tables 
are given both in English and in French. It contains, as do the preceding ones, all the numerical 
documentation on the Emission-Spectra (Dr. Bruninghaus)—Electromagnetodptics (Prof. Wolfers) 
—Zeeman Effect (Dr. P. Auger) and Absorption Spectra (Prof. V. Henri). 

This latter part is an original work of considerable value. 

310 pp. Bound Copy $10. 


Orders ought to be sent directly to the Secretary of the International Committee of ATC. 


11, Rue Pierre-Curie, Institut de Chimie, Lab. d’ Electrochimie, Paris 5° 
accompanied by the amount in a check on a Paris bank, or a money-order in the name of M. C. 
Marie. 

Reprints of Spectroscopy from previous Volumes will be sent to you at a 50 per cent reduction. 
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